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Summary. — We have been registering the intensity of the penetrating 
component of the cosmic radiation since the 20-th February 1958. The 
absorption and decay coefficients as well as the total barometric coefficient 
were determined from the data of the first year of the observation. We 
determined further the amplitude of the fluctuation of the primary radi- 
ation for each month and the amplitude of the solar daily variation. 


1. — Introduction. 


Within the framework of the I.G. Year and the I.G. Co-operation 1959, 
we have been registering the penetrating component of the cosmic radiation 
at a depth of about 18 m underground (i.e. about 40 m w.e.). The primary 
aim of this work was to determine the atmospheric coefficients of the pene- 
trating component and to investigate the solar diurnal variation of the in- 
tensity (*). 

The measuring apparatus consisted of two identical, semicubical meson 
telescopes operating independently of each other. A detailed description of 


the telescopes may be found in (1). 


(*) Part of the results published in this paper have been presented at the IUPAP 


Conference for Cosmic Rays held in Moscow, July 1959. 
(1) T. SAnvor, A. Somoeyr and F. TeLBisz: Suppl. Nuovo Cimento, in press. 


1 - Il Nuovo Cimento. 
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The experimental station is situated in Budapest, at a geographical lati- 
tude 47.5° N, longitude 18.9° E and at 410ma.s.l. It has been operating 
since the 20-th February 1958. 


2. — The atmospheric coefficients. 
We solved the regression equation of the following type: 


(1) a = AAp+DAh, 


0 


denoting by 1, the average sum of the two fluxes measured by the two tele- 
scopes, by AT the deviation of this sum from I,, by Ap the deviation of the 
barometric pressure from its mean, by Ah the analogous quantity for the height 
of the 200 mb isobaric level and by A and D the relative absorption and decay 
coefficients, respectively. 

We determined further the total barometric coefficient (B), defined by the 
equation 

AI 


ee Sag: 
I, a 


TABLE I. 


Ween gs n A D B ll 
% perem Hg, % per km | % per cm Hg 7 
0 70 P Hoo 


1958 March| 178 | 0.77-40.02 | — 0.54+-0.05| — 0.92+0.27 | — 0.69 +0.03 | 2.1.0.2 
April | 156 | 0.84-+0.01 0.84 +0.04 | — 1.87-L0.20 | — 0.84+0.04 | 2.80.3 
May 234 | 0.54-L0.03 0.77 £0.05 | — 0.45+0.10| — 0.73-L0.05 | 2.4+0.2. 
June | 148 | 0.66+0.03 | —0.50+0.04 | — 1.43-£0.23 | — 0.52--0.04 | 2.3-L0.3 


July | 96 | 0.86--0.03 0.96 +0.07 | — 0.55+0.25 | — 1.04-+0.06 | 1.2-+0.2 
Aug. 181 | 0.78+0.02 1.66 +0.07 | — 0.83-40.18 | —1.72+0.05 | 5.40.4 
Sept. | 108 | 0.37+0.04 | —0.32+0.12 | — 1.02+0.18| — 0.45+0,12 | 4:7-+0.4 
Oct. | 50 | 0.7040.05 | — 0.35 +0.09 | — 1.51-+0.21 | — 0.45 -+.0.08 | 1.9-+0.4 
Nov. 98 | 0.79+0.02 | —0.72+0.04 | + 0.36+0.39 | — 0.71-40.04 | 2.4-L0.3 
Dec. 66 | 0.89+0.02 | — 0.71+0.05 | — 2.04+0.32 | — 0.72 40.04 | 1.4-40.3 
1959 Jan. 119 | 0.88+0.01 | — 0.99 40.04 | — 0.67 -£0.16 | — 1.11-+0.03 | 2.9-+0.3 


Febr. | 188 | 0.88+0.01 | — 0.90 +0.04 | — 0.51-40.18 | — 0.97+0.02 | 2.9+0.3 


1-3-1958- | 1622 | 0.830.004 | — 0.78 40.01 | — 1.20+0.02 | — 0.92+0.01 | 4.64+0.1 
28-2-1959 È 
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Fig. 1. — Absorption coefficients. 
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SOLAR DAILY VARIATION ETC. 


1958 1959 
00 A IV VO VI vil vil IX X XI XW 1 


Yearly mean:-1.20+0.02 


RES 


1% per km 


Fig. 2. — Decay coefficients. 


A, D and B have been determined separately for each month, as well as 


for the one-year interval beginning with the 1-st March 1958. 


may be seen in Table I and in 
Figs. 1, 2 and 3. 

Data concerning the _ isobaric 
layers were taken from the Daily 
Report of the Hungarian National 
Institute of Meteorology. At the 
radio sounding station situated 
20km south-east from our labora- 
tory, balloons are flown up twice 
daily. The barometric pressure data 
were obtained from the readings of 
a barometer at the place of obser- 
vation and from microbarograph 
data placed at our disposal by the 
National Institute of Meteorology. 

In calculating the coefficients 
given in Table I only those bi-hourly 
intervals were used during which 
both telescopes functioned properly 
and the appropriate radio sounding 


a 


The results 
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Fig. 3. — Total barometric coefficients. 


4 T. SANDOR, A. SOMOGYI and F. TELBISZ 


data were also available (*). The total number (k) of such intervals is also shown 


in Table I. 

The errors shown in Table I were calculated on the basis of the method 
outlined in papers (1%). We modified the method so as to take into consi- 
deration also the measuring errors of the atmospheric factors. The following 
formulas were obtained and used in the present calculations: 


a i i [r2(SI,)? + A2IG( Spo)? + D?Lo(dhy)? |? , 

(e) i (2%) 
Or1 — \Opp 

FE). à 3 

Vire LEO: EL (= , RC) | (ones L= (Sh 

Je Opp Veh y O11 1 — thn Opp er Onn i 

Sp) = 1 One len (29) ! | 14 rt ee | oe di 

| Io Om VI — re, | Ou TT NER, 1 — ria \om i 


where r denotes the total correlation coefficient of the equation of regression (1); 

Tps Tint “pn denote the simple correlation coefficients between the quantities I 

and p, I and h, p and h, resp.; 0,,, Gr, and o,, are the second central moments 
h Ta ny } & ] Ag i, rh Ing D pe 5 ray 

of the variables I, p and h, resp., i.e. for instance of, = (1/(k —1)) > (Ap). 


Furthermore, 


4 


dlo= | Spo = —= ola 


where dp and Sh denote the measuring error of one bi-hourly value p and 
h, resp. 

The errors in the atmospheric coefficients due to the uncertainty in meas- 
uring the barometric pressure and the isobaric heights turned out to be in- 
significant (*) as compared to the errors due to statistical fluctuations in the 
measured fluxes, although errors as high as 0.3 mm Hg and 40 m, resp. were 
assumed in the barometric pressures and isobaric heights. 


(*) The appropriate isobaric heights were determined by linear interpolation be- 
tween the data obtained by the two radio Soundings preceding and following th 
bi-hourly time intervals. Sac 

(2) L. JANOSSY and G. D. RocHESTER: Proc. Roy. Soc., A188, 786 (1944) 

(3) L. JANossy, T. SANDor and A. Somocyr: Suppl. Nwovo Cimino 8, 701 (1958) 

(eS) The errors of the decay coefficient turned out to be larger by 15% when i | 
errors in the measurement of the atmospheric data were taken into Sastre the err | 
of the absorption and total barometric coefficient remained nearly a et 


2212 


METEOROLOGICAL COEFFICIENTS AND SOLAR DAILY VARIATION ETC. 5 


The coefficients obtained show a much greater spread than would be ex- 
pected on the basis of their statistical errors. A 72 test yields a probability 
less than 10 * for these spreads to be of purely statistical origin. We have 
inclined to ascribe these large fluctuations in the meteorological coefficients, 
at least partly, to correlations between meteorological factors and other 
(e.g. geomagnetic) quantities influencing the cosmic ray intensity. An example 
clearly demonstrating such a case is given in paper (4). 


3. — The amplitude of fluctuation in the primary radiation. 


In paper (!) we have shown a method suitable for the determination of the 
mean amplitude of fluctuations due to instrumental effects as well as for the 
determination of the mean amplitude of the real fluctuations (¢,), including 
atmospheric effects. 

If we carry out the same calculations using the fluxes corrected for atmos- 
pheric effects instead of the uncorrected ones, we get an average amplitude (0,,), 
which—as far as the procedure given in paper (1) is justified—is neither of 
instrumental nor of atmospheric origin, and which amplitude, consequently 
may be regarded as that of the primary radiation. 

There are several ways for obtaining o,,, according to the different modes 
in which the correction for atmospheric effects can be carried out. 


1) If we reduce our measured fluxes by the coefficients defined by equa- 
tion (1) we have simply 


(8) TAI, 


2) If we calculate the meteorological coefficients separately from the 
fluxes measured by the two telescopes we are able to reduce the data of each 
telescope with its own coefficients. 

In this case 


where a; and d,, resp. are the (not relative) absorption and decay coefficients 
belonging to the data of telescope No. à (î = 1, 2); oi, and 0, resp. are the 
second mixed central moments formed by the (uncorrected) rates of telescope 


(4) T. SANDOR, A. Sowogyi and F. Tergisz: Acta Phys. Hung., 11, 205 (1960). 
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No. i and the meteorological data p and h, resp. For instance 


— 1 LE 


O3p = rer ALND 


where I® denotes the uncorrected flux registered by telescope No. 2. 


3) If we use meteorological coefficients other than those calculated by 
means of the appropriate equations of regression, we get 


2 zi 2 2 Desa € DI 2220 
o, = 0° — a(0î, + 0%) dot, + 05, + ao, F2ado,, + doi, 


where a and d are the given meteorological coefficients. 
The values of o,, shown in Table I were calculated on the basis of the for- 
mula (3). Their errors are given by the formula 


Gece 2 
Nos = VIE, Aig? Ol, È 


which was obtained in the same way as the formulas (2), neglecting errors due 
to uncertainties in the atmospheric data. o,, has the following meaning: 


cee es ALOT 
je 


4. — Solar daily variation. 


After having corrected the bi-hourly data for the absorption and the decay 
effects we determined the amplitude and the phase of the solar daily variation 


TABLE II. 


Month k' | A (%) pe URES 
1958 March 295 1140.2 1240 +0040 
April 300 | 1.6.0.2 1020 +0030 

May | 355 1.0-+0.2 0920 +0040 

June | 291 23.0.2 1230 +0020 

July 282 ey 1540 +0110 
August 293 | 0.60.2 1550 +0130 
September | 347 0.805 1830 40120 
October | 314 | 1.0+0.4 2150-+0120 
November | DTT 14202 1150 40050 
December 318 0.5 +0.3 1550 +0220 

1959 January 308 | 1.1 +0.3 1610 +0100 
February | 304 0.8 -+0.2 1000 +0110 
1-3-1958 - 28-2-1959 3 684 0.68 +0.08 1700 +0030 
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= 


for each monthly period as well as for the whole year of observation. The 
values obtained are shown in Table II and Fig. 4. Every such bi-hourly 
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ee 
Fig. 4. - Harmonie dial. 


interval was used during which at least one of the telescopes functioned pro- 
perly. The total number (k’) of these intervals is also shown in Table II. 


RIASSUNTO (°) 


Sin dal 20 Febbraio 1958 siamo andati registrando l’intensità della componente 
penetrante della radiazione cosmica. I coefficienti di assorbimento e decadimento, 
come anche il coefficiente barometrico totale, sono stati determinati in base ai dati 
del primo anno di osservazioni. Abbiamo inoltre determinato l’ampiezza della fluttua- 
zione della radiazione primaria per ogni mese e l’ampiezza della variazione solare 
giornaliera. 


4 (*) Traduzione a cura della Redazione. 
a 
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Interactions of 6.2 GeV Protons in Emulsions. 


I. — General Results. 


H. WinzeLER, B. KLATBER, W. Kocx, M. Nrxorié and M. SCHNEEBERGER (*) 


Physikalisches Institut der Universität - Bern 


(ricevuto 111 Marzo 1960) 


Summary. — About 700 interactions were found by following 6.2 GeV 
protons « along the track ». An interaction mean free path of (38.2 +1.5) cm 
was obtained. 42 possible proton-free proton collisions were separated 
and the cross section for meson production was found to be 23%{ mb. 
We found an average multiplicity for charged secondary tracks in inelastic 
p-p collisions of 2.8--0.3 and the branching ratio of 2 prong to 4 prong 
events was about 2. These values were compared with the predictions 
of the statistical theory and with a modification of it. The results of 
this investigation would be consistent with a model which assumes 
that at these energies about 50% of all collisions are « peripheral» ones 
which lead to (nucleon) isobar (N.N* or N*N°*) formation only. The 
angular distribution of all thin secondary tracks in the L system have 
been measured and the formula of Castagnoli et al. has been applied 
to a reasonably selected sample of stars. This gave an overestimation 
by a factor of about 2 for the primary energy. Other investigators 
have obtained a similar value at 101? eV. 


Introduction. 


A large amount of experimental work has been performed in the field on .N-N 
collisions at different energies. Several experimental data in the (1-10) GeV 


(*) Now at CERN, Genève. 
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region have been published (+), at 101 eV the available statistics is still con- 
siderable (115) and single events up to 10! eV have been described (1115). 
The appearance of singly charged particles with even higher energies is prob- 
able (15). 

In ref. (5°) and in the present work the proton energy was 6.2 GeV. The 
Dubna machine reaches 10 GeV and the CERN proton synchrotron provide 
us with 28 GeV protons. All these energies can be considered as passage to 
cosmic ray-jet energies (*). Pion, strange-particle and anti-particle production 
increases and, in the framework of general program, the study of its energy 
dependence is of high interest. Unfortunately the average number of particles 
produced is still small even for the highest machine energies and this in some 
Sense complicates the situation. However, with machines one can have 
practically any particle flux and the primary energy is precisely known. Emul- 
sion has the disadvantage of being a complex nuclear mixture, but it seems 
to us, that this difficulty can be overcome and that one can separate, for 
example, the collisions with free protons in a quite satisfactory manner. Be- 
sides this, even collisions with bound nucleons can give much information about: 
the nature of p-N collisions. 


() W. B. FowLER, R. P. SHUTT, A. M. THORNDIKE and W. L. WHITTEMORE: 
Phys. Rev., 103, 1479 (1956). 
(MMM Brock, EE NT ARTE VI I COCCONL BE. Harry W. Bs dowLHR, 
R. P. SHuTT, A. M. THORNDIKE and W. L. WHITTEMORE: Phys. Rev., 103, 1484 (1956). 
(3) W. B. FowLER, R. P. SHuTT, A. M. THORNDIKE, W. L. WHITTEMORE, V. T. Coc- 
coni, E. HART, M. M. Brook, E. M. Harru, E. C. FowLER, J. D. GARRISON and 
T. W. MORRIS: Phys. Rev., 103, 1489 (1956). 
) R. CESTER, T. F. Hoane and A. Kerman: Phys. Rev., 103, 1443 (1956). 
) R. E. CAvanoUGH, D. M. HASKIN and M. ScHEIN: Phys. Rev., 100, 1263 (1955). 
6) U. HaBeR-ScHarm: Nuovo Cimento, 4, 669 (1956). 
) F. N. Hormquisr: UCRL 8559 (1958), p. 7. 
) R. M. KaLBac4H, J. J. Lorp and C. H. Tsao: Phys. Rev., 113, 325, 330 (1959). 
(°) R. R. DANIEL, N. KAMESWARA Rao, P. K. MALHOTRA and Y. TSUZUKI: preprint 
(Bombay, 1959). 
(10) E. M. FRIEDLANDER: Nuovo Cimento, 14, 796 (1959). 
(11) N. P. BogacHEv, S. A. Bunyarov, I. P. MEREKOV and V. M. Srprrow: Dokl. 
Akad. Nauk, 121, 615 (1958). 
(2) M. TeUcHER: Conference on High Energy Accelerators (CERN, 1959); pal: 
(13) Proc. of the 1958 Conference on Elementary Particles (CERN, 1958). Report 
of PICCIONI. 
(4) W. B. FowLER: Proc. of the VII Annual Rochester Conference (1957), p. x1-18. 
(15) M. W. TrucHER, E. LOHRMANN, D. M. HASKIN and M. ScHEIN: Phys. Rev. 


Lett., 2, 313 (1959). co 
(16) G. CLARK, J. EARL, W. KRAUSHAAR, J. LinsLy, B. Rossi and F. ScHERB: 


Nuovo Cimento, 10, 630 (1957). | , + 
(*) In the following called also «ultra high» energies. Energies above 10°? eV. 
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Theoretical predictions can be made according to various theories (RIO), 
Numerical calculations for the above mentioned machine energies, following 
the Fermi statistical theory, have been performed (225). In its application to 
different primary energies the statistical model is practically limited, since 
with increasing primary energy the number of final states to be taken into 
account increases rapidly. Therefore at cosmic ray energies Fermi applies 
the «thermodynamical » (*) treatment. 

For. 6.2 GeV protons, involving certain assumptions, reaction rates for 
pion, strange-particle, anti-baryon production and their energy spectra can 
be calculated with the low energy theory (2), and it is interesting to com- 
pare these with the experiment. Furthermore angular distributions can be 
used to get indications about the geometrical side of the .N-.N interaction 
models, and in this connection one has probably to distinguish between head-on 
and peripheral collisions. Correlated to this and of a certain practical impor- 
tance is also the average fraction of primary energy spent in the creation of 
particles, the «inelasticity ». 

For cosmie ray energies, on the other hand, the situation is kinematically 
(all fs & 1) and perhaps theoretically (thermodynamical treatment) less com- 
plicated than in the low energy cases, but much more undetermined on the 
experimental side, since the primary energies are not known a priori. In ad- 
dition, in most cases, scattering measurements are not sufficient to determine 
the energies of all charged secondary particles and there exist many jets where 
no scattering measurements on the secondary tracks have been possible at 
all. In those cases the main way to determine the primary energy is via the 
angular distribution of the secondary tracks. The formula used by CASTA- 
GNOLI et al. (26), for instance, gives relatively reliable results for the average 
energy of a larger sample of similar interactions, but in a single case it can 


(7) W. HEISENBERG: Zeits. f. Phys., 126, 569 (1949); 133, 65 (1952). 
(18) E. FERMI: Progr. Theor. Phys., 5, 570. (1950). 
(?) L. D. LANDAU: Dokl. Akad. Nauk, 17, 51 (1953); S. Z. BELEN'K1IJ and L. D. 
LANDAU: Suppl. Nuovo Cimento, 8, 15 (1956). 
(2°) R. HAGEDORN: Nuovo Cimento, 15, 246 (1960). 
(2!) J. von BEHR and F. CeRULUS: private communication (1960). 
(22) R. HAGEDORN: CERN 59-25. 
(2°) F. CERULUS and R. HAGEDORN: CERN 59-3. 
) R. HAGEDORN: preprint CERN (1959). A new derivation of the statistical 
theory of particle production with numeral results for p-p collisions at 25 GeV. 
(CSS MM BEL ENTRIT, Vi. M. MAKSIMENKO, A. J. Nrkisoy and J. L. ROZENTAL: 
Fortschr. Phys., 6, 524 (1958). 
(26): C. CasraGnoLI, G. CORTINI, C. FRANZINETTI, A. MANFREDINI and D. Moreno: 
Nuovo Cimento, 10, 1539 (1953). 
(*) Since the Stephan Boltzmann law is applied. 


24) 
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fail by an order of magnitude (27). This can partially be ascribed to statistical 
fluctuations in the charged particles angulary symmetry, (forward-backward 
Symmetry is assumed in the derivation of the «Castagnoli » formula, it is 
evident in the C.M. system for p-p but not for p-n collisions), and to the fact 
that the complete formula contains corrections which depend on the energies 
of the produced particles. The large fluctuations are also due to a lack of 
absolute selection criteria for p-free p collisions in emulsions. Such topics will 
be stressed in the present paper. Especially the « white » stars (stars with no 
gray and black prongs) and stars with only one or two gray and black prongs 
will be discussed. 

With respect to the present work two different ways of investigation are 
available. 


a) Scattering measurements can be performed on the secondary tracks. 
Particle identities are assumed and the kinematics are transformed into the 
C.M. system. This gives the maximum information and naturally questions 
like that of the inelasticity can be answered in this manner. 


b) If one does not perform scattering measurements one can still make 
use of the lab-angular distribution, but then one has to assume certain distinct 
models for meson production and to compare its predictions with the exper- 
imental curves, which in this case exist as angular distributions for different 
shower particle multiplicities. On the average for lower multiplicities the pro- 
duced particles are much more strongly collimated than for larger multiplicities. 
This might give means to approach the geometrical aspect of the interaction, 
namely to distinguish between peripheral and head-on collisions. 


In the present paper we have proceeded according to method 5). Scat- 
tering measurements are in progress and the results will be published in a 
forthcoming paper. 


1. — Experimental procedure. 


6 plates Ilford G-5, (10x15) cm?, 600 um thick, exposed to the 6.2 GeV 
Berkeley internal proton beam, were scanned by the usual « pick-up and fol 
lowing along the track » method. The average proton density in the plan per- 
pendicular to the beam entering the emulsions was 1.1-10° protons per em?. 
The vertical projection of the angular distribution of the beam where it entered 
the plates and after having penetrated 7.6 cm is shown in Fig. 1. The tracks 
were extremely flat. Their average potential path per plate was more than 
6cm. The tracks, being slightly inclined from glass to surface, were picked 
up near the glass, noted individually and followed until they either left the 


(27) M. W. TEUCHER: private communication (1959). Derived for ultra high energies. 
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plate or led to an interaction. The speed of following through was so fast 
that the scanner observed only the motion picture of the track. Every «scat- 
tering » with more than 1° in projection was noted. 


% — At the entrance 
2 (o- 9.89) 
——After 76 cm path 
22 (0-14.78) 
18 
14 
10 
6 
2 


+----- + dee 


EU) ZIO 0 20 - 0 (0) O 2 O0 nr 


Fig. 1. — Projected angular distribution of the beam tracks. 


Sometimes the scanner jumped to an already followed track while fol- 
lowing through which resulted in 28 of the 699 stars being found twice. 
This fact has no further influence on the considerations if taken into account 
in the right manner. The 28 events are not included twice in the diagrams. 


2. — Results. 


Four approximately equivalent samples obtained from 5 scanners contrib- 
_uted to the statistics and a total path of 261m has been followed. Table I 


TaBLe I. — Raw observed values. 
Nica bent Average Possible elastic 
Scanner d À Path (m) interaction p-free p 
interactions hes 
length (em) collisions 
B.A. 199 7205 36.4 2 
(CI, 150 56.80 37.9 1 
BD: Il 
EB. f 201 74.96 37.3 3 
H.W. 149 56.97 38.2 3 
Total 699 261.24 Bical (13) 9 
(*) Uncorrected, see sub-sections 2°2.1 and 2°2.3. 
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contains the raw observed values. As can be seen, the individual fluctuations 
are small. This indicates that the influence of «human failure » is probably 
also small. 


2°1. The p-free p collisions. — A first cut was made in order to separate the 
p-free-p collisions. We observed: 
a) 315 stars with more than 10 prongs; 
b) 347 events with less than 11 prongs and 
¢) 9 possible elastics p-free p collisions. 
21.1. Elastic p free p collisions. — Table IT shows the angles of the 
9 events of class c). They were selected according to the following criteria: 
1) Coplanarity ; 
2) te dy te Oy = 1 — fo; 
3) Energy of the slower proton. 


Tage II. — Angles in the CM system of 9 possible p-free p elastic collisions. 


Event Angle in the CM system 
H 84 4.5° 

B 144 5° 

Ch 19 5.5° 

Br 71 10.5° 

H 137 10.5° 

H 10 ye 

B 29 18° 

Br 28 18° 

Li 34 22° 


Unfortunately there is a strong collimation in the forward-backward di- 
rection and this smoothes out the selectivity of the 3 checks. If most de- 
flections have small angles, test 1 is very rough for geometrical reasons. Con- 
cerning test 2, a small angle deflection of the fast particle corresponds to an 
80°-85° deflection (in the L system) of the slow proton. Small linear distor- 
tions in the emulsion (which, for instance, deform a square into a rhombus 
in the horizontal plane) can cause a larger error in the tangent of the larger 
angle, even though the angle itself may be measured quite precisely (+ 0.1"). 
The same is valid for test 3. Relation 3) is the correlation between the mo- 
mentum of the (slow) proton and its scattering angle: 


pi : cos 97 
FE. ee eee 1 — ((Yprim. proton” 1)/2) sin? ye 


N 
ral 
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In 5 of the 9 cases the proton came to rest in the emulsions. Possibly some of 
the 9 events are only quasielastic collisions, but no statistical separation can 
be made here. 


2:1.2. Inelastic p-free p collisions. — Table III demonstrates how 
the events in class b) have been handled in order to separate the inelastic 
p-free p collisions. An inelastic p-free p collision is, by definition, a collision 
in which at least one particle has been created. 


TABLE ITI. — Separation of inelastic p-free p collisions. 
| N, 
Number of 2 | N; | NT UP 
Sera ao | ar with | cp-free p ela- | Stars with2 | N; È 
ee oy peo cae han 3) Stio» criterion slow p°s Ellipsoid test, 
> grano (ee | exe | si | 0l1|2|3|4|5|6|7 
| À | : 
0 1 
il 40 
| 29 | clean 
2 26 26 24 24 21 — = 
| | 33 | © dirty 
| &) @ clean 
3 55 54 38 35 : Re | o n 
1193) i dirty | 
413 clean 
4 45 32 16 1 11 - - - | 
By 2 dirty 
| || I clean 
5 33 19 7 4 = = 
| 1 dirty 
clean 
6 37 14 4 3 1 = = 
: 1 dirty 
7 9 clean 
5 7 3 1 = = 
? IN dirty 
8 32 1 1 0 È E RIE Sg 
9 28 2 1 0 
10 25 0 0 0 
5 E | 
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N, is the number of stars with less than 3 certain nucleonic tracks. 


N, Comparison with the p-free p elastic kinematics has been made and all 
events with a proton faster than the same-angle elastic proton have been 
eliminated. 


N, remains if those cases in which each of the two protons have less than 
195 MeV kinetic energy are excluded. One can easily show that this is 
a limiting value for inelastic p-free p collisions. 


Up to this column even- and odd-prong stars have been treated equally. 
The next criterion is only applied to the even-prong events. 


N; Many of the remaining stars still show a slow proton. N; remains if all 
events have been excluded in which the proton momentum vector did 
not fulfill a certain angular condition. This means: if an event did not 
belong to an elastic collision and if it occurred on a free proton, at least 
one pion must have been produced. Consequently for such cases a validity 
volume for the endpoint of the proton momentum vector can be calcu- 
lated, which gives a stronger condition than that obtained for N,. The 
L momentum vector of a proton with a given C.M. energy E°" describes. 
an ellipsoid defined by: 


pr yt Hy Ben 0080" + V (PE) — (Bouyémmn)? sin? 9° 
SEE dar (Béuyou)” sin? a i 


Upper indices denote the frame of reference, Æ is the total energy, P the 
Momentum, f and y are the velocity parameters and 2° is the angle between 
I momentum vector and primary direction. If E°" > ES, the ellipsoid cor- 
responding to E° lies entirely within the one for ET“. So all protons with 
0< EM EM must have L momentum vectors which lie within the ellipsoid 
for EM. For a given pion multiplicity n,, POM is determined by: 


(Eu)? ae mi = (My + Nr) 


EX 
max —_ 9 FCM 
2E Gta 


The pion multiplicity n, can be replaced by the minimum number #, of the 
charged pions. The presence of neutral pions causes a decrease of E and 
the criterion is fulfilled a fortiori. Thus two-prong stars which were not in 
agreement with the elasticity test were subjected to the ellipsoid test with 
n, =1. For the 4-prong stars one can insert n, = 2, etc. 

33 stars are left which can belong to inelastic p-free p collisions. We dis- 
tinguished between «clean » and «dirty » events. In our terminology a clean 
event is neither accompanied by a blob nor by an electron. As one sees from 
Table III there still exists a high fraction b = 15/40 of clean events among 
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the odd prong stars (which certainely belong to collisions with bound nucleons). 
This indicates that a certain number of clean events among the even-prong 
stars can belong to complex nuclei collisions. On the other hand one must 
algo assume a certain number of spurious blobs and electrons among the 
inelastic p-free p collisions, (about (0--10)% (*), see fi. ref. (?8)). Thus the 
fraction f of clean events among the p-free p collisions is smaller than 1. With 
the knowledge of f and b one can calculate in a statistical manner the real 
number x, of inelastic p-free p collisions using the equations: 


æ, + a, = Number of even-prong stars = possible inelastic 
p-free p collisions. 


fv; + ba, = Number of clean even prong stars. 


With b = 0.37 and f = 0.95 (**) we obtained: x,=19. With this the inelastic 
cross-section for p-free p collisions becomes: 


= +7 
ON ZO SIDE 


The method of separation applied here is a statistical one and does not 
enable us to decide in detail which of the 33 possible p-free p inelastic colli- 
sions are true p-free p collisions. So far the branching ratio (see Table III) 
of charged particle multiplicities—f.i. the ratio 2-prong stars/4-prong stars = 
= 21/11 — necessarily contains a contamination from p-bound nucleon col- 
Jisions. 

Thin track multiplicities and angular distributions for the 33 possible in- 
elastic p-free p collisions are given in the next sub-section. 


2:2. The collisions with bound nucleons. — Fig. 2 shows the frequency distri- 
bution of 662 stars with n, (grey+black) prongs. 319 stars show more than 
6 grey +black prongs. With the rough assumption that most of the grey +black 
prongs belong to nucleonic particles (protons, deuterons, tritons, alphas, etc.) 
one can say that most of these 319 stars are due to the heavy emulsion nuclei 
(Ag and Br), which contribute about 70% to the total geometrical cross- 
section. We differentiate roughly in the following between heavy and light 


(8) Y. EISENBERG, W. KocH, E. LOHRMANN, M. NikoLié, M. SCHNEEBERGER and 
H. WINZELER: Nuovo Cimento, 8, 664 (1958). 

(*) In principle this percentage can be obtained from the elastic p-free-p collisions. 
Among our 9 cases only one event had an associated blob. 

(*) Naturally f depends on the kind of stack, i.e. on exposure and development, 
but the results are insensitive to its precise value since it is large anyway. 
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nuclei by attributing the 
Stars with n,>6 to p-bound oo 
nucleon collisions in heavy 

nuclei. (Concerning Fig. 2 80 
it must be said that usually 
such distributions are plot- 
ted on a semilogarithmic 
scale (2° 39), but our available 
Statistics is not large enough 50, 
to expect additional infor- 
mation from this.) Fig. 2 40. 
also contains the distribution 

of n, for stars with one 3 
thin track (nm; =1). 


39'Stars' with one thin track only 


Change of DE interval 


] - 1Star 


WA Stars with N= 1 


' 
resterà 


Fig. 2. — Frequency of stars 
with », gray+black prongs. 


0 ON COTE NET"; a 


22.1. Thin track multiplicities. — Table IV shows the thin track 
multiplicity distribution for the 33 possible inelastic p-free p collisions and 
for 94 stars with n,= 0, 1 and 2 (see discussion). The corresponding average 
n, are 2.06 and 2.09. In the latter 94 stars p-free p collisions are not included. 


TaBze IV. — Thin track multiplicities for the 33 possible inelastic p-free p collisions and 
for the stars with n,=0, 1, 2. 


n 0 1 2 EN PE NAT 6 7 


Inelastic p-free p collisions 2 12 | 6 6 0 0 0 


stars with n,=0, I, 2, 
excluding p-free p collisions 


Or 
VO) 
© 
i 


0 20 30 35 


A summary of the raw observed number of thin tracks emerging from 
662 stars is given in Table V. An increase in thin track multiplicity with in- 
creasing number of gray +black prongs is indicated. The effect, which must 


(29) M. W. Trucner: Zeits. f. Naturfor., 8a, 131 (1953). 
(50) W. HEISENBERG: Vorträge über kosmische Strahlung (1953), p. 79. 


2 - Il Nuovo Cimento. 
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be due to secondary collisions inside the hit nuclei, becomes more obvious 
from Fig. 3 which shows the frequency distribution of the stars as a function 
of the thin track number n, for heavy (n,> 6) and light (7, < 7) nuclei. 
For n,<7 and > 6, ©, becomes 2.27 and 3.04 respectively. 


Tage V. — Thin track statistics as function of gray+black prong number (9 possible 
p-free p not included). 
SR Number Thin tracks, Average number 
È of stars total of thin tracks 
0 70 131 1.86 
1 53 107 2.02 
2 45 95 2.11 
3/4 86 192 2.23 
5/6 89 257 2.89 
7/8 47 123 2.62 
9/10 46 147 3.20 
11/12 38 128 23) 
13/14 33 101 3.06 
15/16 30 89 2.97 
17/18 30 93 3.10 
19/20 27 80 2.96 
21/22 20 50 2.50 
23/24 23 74 3.22 
25/30 25 89 3.56 
Total 662 1753 2.65 


It seems to us necessary to make some remarks about the notions « stars » 
and «interaction ». In the case of the p-free p collisions the terminology is 
clear. We have distinguished between elastic and inelastic collisions as dis- 
cussed in subsection 21.2. For proton collisions with complex nuclei, how- 
ever, there is some arbitrariness. We call «star» every event which has 
at least one gray +black or at least 2 thin secondary tracks. We have also 
noted 39 «scatterings », the projected scattering angles of which are > 1°. 
In a certain fraction of these events the only difference between the « stars » 
with n,—1 and n,> 0 is that neutrons have replaced the charged evaporation 
products. For uniqueness of definition of an «interaction» we perform a 
suitable extrapolation by using Table V. As can be seen from this table the 
average number of thin tracks for events with », — 0 is lower than the neigh- 
bouring values for stars with n, = 1,2 and more. We perform here an extra- 
polation by postulating the average thin track multiplicity for the events 
with n, = 0 to be 2.1, which is an average value for the stars with n, =1, 2 


2226 


INTERACTIONS OF 6.2 GeV PROTONS IN EMULSIONS - I 19 


and 3. This extrapolation results in only 55 instead of 70 events with n, — 0 
having to be taken as «stars». Thus this procedure can be used for defining 
the notion of «star ». 


| 
Number of slars with ns thin tracks 
ol 
Ha 
mot | 
I fa 
\ 15 
i (oO 
10 ey & 
ES 
È I 
S 
CR o 
no CI 
Ù ' | 
80! ‘: SS} === Stars with less than 
| Il 
A | 7(gray+black) prongs 
n mel 
60! 319 Stars with more than 


e 6 (gray + black) prongs 


507 


40; 


10 


+—_—_ + 

O e ee 

Fig. 3. — Frequency distribution of 662 stars with n, thin tracks (excluded 9 possible 
p-free p elastie collisions). 


22.2. Thin track angular distributions. — Fig. 44 show the 
lab-angular distributions for different thin track multiplicities separately for 
heavy and light nuclei. Fig. 4-b correspond to the thin tracks of the 33 pos- 
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Sible inelastic p-free p collisions. Fig. 4-c consist of the angular distributions 
of the thin secondary tracks of the stars with n, — 0, 1 and 2. Such stars are 


Ng = 6 
aS el O 
; [L Se Dr DE Pr PORRI 
(0) 20 40 60 80 100 120 140 160 180 
a) 
np? 6 


+ + fe sas =} 4 
0 20 40 60 80 100 720 140 160 180 


Fig. 4. — L-angular distributions of thin tracks. 


often in the analysis of cosmic ray induced jets as representative for « pseudo- 
nucleon » collisions. The distributions will be discussed in detail in Section 3°2. 


22.3 Interaction mean free path. — We observed 699 interactions 
in a total path of 261m. After subtraction of 15 events with one thin track 
only (see sub-section 2°2.1), we get an interaction mean free path of 


A == (36.2 = 1.5),.cm.,. 


3. — Discussion. 


In the following we compare the results obtained from the p-free p colli- 
sions with numerical results of calculations according to the statistical model. 
In addition a rough « velocity » test on the secondary particles is performed. 

The collisions with bound nuclei, their thin track angular distributions and 
number of gray +black tracks are discussed in sub-section 3°2. 


31. The p-free-p collisions. — The statistics obtained for the p-free p col- 
lisions is small but within its limits the elastic and the inelastic cross-section 
of 23 mb agrees with a corresponding value obtained by KALBACH et al. (8) 
who applied a completely different scanning method. For the relative fre- 
quencies of 2 and 4 prong events however we obtained approximately the 
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reciprocal of their ratio; namely 15 mb and 7 mb, respectively. This point 
has to be clarified with better statistics. 

On the other hand we can compare our values with the results of extensive 
numerical calculations for 6.2 GeV protons performed by CERULUS and HAGE- 
DORN (?°). PAR 

These authors used the Fermi theory in the low energy form (statisticai 
model) with the expression: 


T(28j + 1)% 
P,=1 SL) E i | fae CE, Ma, Mn, P = 0), 


P, is the probability for a distinct final state b. W takes account of the iso- 
topic spin T conservation. It is the probability for » particles of 7-3, s par- 
ticles of T —1 and ¢ particles of the T= 3 to form the state of total isotopic 
spin 7. The next factor stands for the spin multiplicity, Q for the interaction 
volume and o for the phase space weight. 

The calculations available (2!) yield different results if the pion-pion isobar 
is included or not. The corresponding average multiplicities for charged sec- 
ondary particles become 3.71 and 3.21. The observed average was 2.8, which 
might indicate that the Fermi theory overestimates the multiplicities at these 
energies. Furthermore the Fermi theory predicts reaction rates of 2, 4 and 
6 prong events which do not agree with the experimentally observed ones. 
The reaction rates predicted with the Fermi theory (we shall call the Fermi 
theory model A) for charged particle multiplicities as calculated by HAGE- 
DORN, V. BEHR and CERULUS are shown in Table VI. 


TABLE VI. — Branching ratios for the emission of charged secondary particles according 
to the statistical theory. 


| Predicted percentages | 
| with nucleon isobar, with nucleon isobar 
| without pion-pion isobar and pion-pion isobar 
2-prong stars 43.4 25.5 | 
| 4-prong stars | 52.6 63.5 | 
| 6-prong stars | 4.0 11.0 


The bad agreement between theory A and the experiment indicates that 
an interaction of the Fermi type alone is probably not able to reproduce the 
experimental observation satisfactorily. 

A more refined model (model B) was considered by v. BEHR and Cr- 
RULUS (21). In that case the validity of the Fermi theory was presumed for 
« head-on collisions » which in this model contribute, by definition, à certain 


a 
N 
nN 
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fraction a to the total inelastic cross-section. Peripheral collisions were treated 
differently. They were only allowed to form one or two nucleon isobars (3-3). 
The relative frequencies of the 3 types of interactions (Fermi, two isobars: 
N*N*, one isobar: N.N*) are called a, b and c. 

There exist different possibilities to evaluate a, b and ¢ from the experiment. 
With sufficient statistics this model will easily prove or disprove itself. 

A first method to get an idea-about a, b and ¢ uses the reaction rates of 
the charged secondary prongs. According to v. BEHR and CERULUS a, b and e 
could be completely determined by the following equation system. 

Without pion-pion isobar, but with nucleon isobar: 


0.434a + 0.8b + ¢ = fraction of 2-prong stars = 21/33 — 0.64 


0.526 a + 0.2 = fraction of 4-prong stars = 11/33 = 0.33 
0.049 a = fraction of 6-prong stars = 1/33 = 0.03 


With pion-pion isobar: 


0255a + 0.8b +e =0.64 
0263 DIM ZE)EZIO = (LS 
0.119 a 40509 


i 4 i 
Contrib. from Fermi-part N*N* NN* 


The six prong events apparently are so rare that one cannot make use of 
the third equation. Since our main interest is directed to the magnitude of a, 


With nucleon 


/ isobar 


/ with nucleon- 


ione) c 
0.2 and pion-pion isobar Th 
+— +— ce 
0 1 2 6 
ELIS 


which stands for the abundance of the Fermi part, we plotted a as a function 
of cJb in Fig. 5. The figure shows (according to our small statistics) that a 
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can vary from 0.30 to 0.63. With the reasonable assumption that ¢/b is about 
unity and presuming the validity of model B one can say that head-on col- 
lisions contribute only about 50%. So much for the analysis of the branching 
ratios. 

There exist other possibilities for getting a, b and ¢ from the experiment. 
For discussing one of them we consider the situation in the C.M. system 
according to model B. Fig. 6 shows a diagram with 8°" and J™ on the axes. 


og! Protons Protons 


Gray track region 


0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 


Fig. 6. 


The shadowed area on the right contains those particles which lead to gray 
or black prongs in the L system (meaning that p’ < 0.7). For head-on col- 
lisions isotropy in the CM system is assumed. For NN* and N*N* formation 
isotropic decay in the system of the isobars is assumed and further that the 
transverse momentum of these isobars in the CM system is small, i.e. that 
their line of flight is practically identical with the primary direction. With 
this crude picture one gets the CM angular kinematics for peripheral collisions 
shown in Fig. 6. According to this the nucleon of the backward isobar (or also 
the backward nucleon in the case of N.N* formation) always produces a gray 
or black track in the L system. For head-on collisions, on the other hand 
with isotropy in the CM system, 98 % of all tracks SAI be thin in the L 
system. 

We now make use of the real observed percentage of gray +black prongs 
from the p-free p collisions in estimating a, b and pai V. BEHR and CERULUS 
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give the following reaction rates for the decay of N'*N°* and NN°*: 


TABLE VII. — Branching ratios for the isobar decay. 


{p fptretn: 0.18 

Roe | Deir 02.0) 
19) 

e tinte: 0.58 


[nn m'nt: 0.04 


vs fp pr: ON 
È lptn+a: 0.83 


One can easily see that with the above considerations one gets the two 
equations: 


0.98a'+ 0.72b' + 0.71¢' = 0.74 = fraction of observed thin tracks , 


0.98a'+ 0.330'+ 0.426’ = 0.37 = fraction of observed white stars , 


We distinguished between a, b, c and a’, b’, c'. If really, as it was assumed, 
all reactions from head-on collisions had lead to isotropic particle emission 
the primed and unprimed values should be equal. If one obtained a much 
smaller value for a’ than for a this would indicate that not all interactions of 
the « Fermi-part » could show isotropic particle emission. The first equation 
is not very sensitive to the real value of a’. The second one looks better but 
there the number of white stars on the right is only 12 in our case. It is thus 
meaningless to calculate a’ from these equations. We can only say that the 
observed fractions on the right hand of the equations seem mostly to be deter- 
mined by d' and ¢’ and this might indicate a small value for a’. 

The method applied above is some Kind of integrate velocity test. It was 
performed visually and can thus only be considered as an estimate. While 
for scattering measurements only flat tracks can be taken into consideration 
one could here use all tracks. (In our experience a practiced scanner can judge 
relatively reliably whether a track is gray or thin, and moreover there are 
only a few «border line » tracks.) 

According to v. BEHR and CERULUS 20% of the fraction b decay into 
p--p-+x++7-, and practically all these events should yield stars of the type 
nn, =1 and n,=2. We observed 5 events of this type among 33 possible 
p-free p collisions. If the Fermi-part shows throughout isotropic emission, we 
can estimate b to be ~ 25/33. These cases will be investigated in more detail 
in a forthcoming paper. The events n, —1 and n,=3 will serve especially 
to study N*N* formation. 
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Finally in Fig. 7 we show the L angular distribution of the 3 thin tracks 
from the reaction: pipoN*+N*>p+p+r'+= as ‘alculated with the 
model discussed above. The strong peak in the forward direction arises from 
the protons of the forward isobars, while, as already mentioned, the Dos 
of the backward isobars, lead to gray tracks. We shall use this figure in the 
discussion of the angular distribution of the thin tracks from all stars in the 
next section. 


100! Mainly from the decay protons 
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of the foreward isobar 


50} from decay pions of the foreward isobar 


From decay pions of the backward isobar 
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0 + + + + + + + + + = À = 
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Fig. 7. - Calculated L-angular distribution of thin tracks for N°*.N°* formation according 
to model B. 


Concerning this sub-section, we cannot decide definitely whether the ob- 
served facts indicate a high abundance of N°N* and N"N* from peripheral 
collisions or whether they indicate a strong deviation from isotropy for protons 
in head-on collisions. If in the Fermi-part the protons tend to keep their pri- 
mary directions, one could obtain the observed percentage of gray tracks 
without the assumption of a large fraction of nucleon isobars. However with 
better statistics it should then be possible to decide by means of the observed 


average multiplicities of charged secondary particles which of the effects is 
predominant. 
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32. p-bound-nucleon collisions. 


82.1. Mean free path for interaction. — We begin with a consider- 
ation that makes use of the mean free path in the emulsions for a derivation of 
the cross-section for «star production » of a single nucleon. With the assump- 
tion of uniform nucleon density, of Hofstadter’s nuclear radii and standard 
emulsion composition one obtains a mean free path in nuclear matter and thus 
a cross-section of 38 mb for the single nucleon. (Interaction m.f.p. 38.2 em). 

This value has to be compared with the p-free p data. We obtained 9 elastic 
and 19 inelastic p-free p collisions (see sub-section 2°1.2) and a reasonable as- 
sumption is that every elastic and inelastic p-bound nucleon collision has lead 
to a star. One would therefore say with this approximation, that the 38 mb 
obtained indirectly must be compared with the 34 mb total p-free p cross- 
section. The only assumption is that p-p collisions and p-n collisions are equi- 
valent at these energies. 


32.2. The angular distribution of thin tracks and the ap- 
plication of the « Castagnoli formula». — In Figs. 4 the L angular 
distributions for different thin track number are shown. Obviously we have 
stronger collimation for lower multiplicities and for small »,. This effect, on 
the one side, must be ascribed to the nature of the p-nucleon collision and 
on the other hand, to secondary collisions inside the struck nucleus. Since 
the difference between the stars with n,< 7 and n,>6 is so evident, one 
can separate the whole sample of observed stars into two groups. The one con- 
sists mainly (if »,< 5) of such events which, in first approximation, can be 
considered as répresentative of the angular distribution of the p-free nucleon 
collisions while the other group consists mainly of stars in which secondary 
interactions happened. At the moment we are not much interested in the 
latter sample. If we compare the angular distribution of the stars with n,< 7 
with that from the p-free p collisions, we see that both show a similar sensitive 
dependence on »,, while the corresponding angular distributions for the stars 
with n,>6 are highly insensitive to ny. It therefore seems to be justified 
for the purpose of this discussion to accept the angular distribution of the thin 
tracks of the stars with n,< 7 as a p-free nucleon like angular distribution. 

For the stars with n,< 7 the formula used by CASTAGNOLI et al. (28) has 
been applied. The general expression. is 


> log tg 97 + log C, 


CN À 


1 
log vom = — — 
W 
and we took the first order approximation with log? —0. The « derived » 
energies for the primary proton are given in Table VIII. In principle this 
formula is «valid» when f%,, = fe: In our case BZ, is 0.88 while foto 
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can vary within relatively wide limits. In addition a necessary condition is 
forward/backward symmetry in the CM system. But here we are considering 

a mixture resulting from p-p and p-n collisions. For the latter type of inter- 
action symmetry is a poorly fulfilled assumption (3), not a logical necessity 

as for the p-p collisions. Several other objections can be made but now we 
are interested only in qualitative results. 


TABLE VIII. — Primary energies as calculated with the formula of Castagnoli et al. 


Ns (nr < 7) Eee proton (GeV) 


470 
26 
15.7 


H> wh = 


Sì © 
| 
D A © 


Looking at Table IV one can conclude that for the p-free p collisions at 
these energies more than 4 thin prongs are rather infrequent. One can therefore 
see that the stars with n,>5 are no longer typical ones. Probably the 
« Castagnoli » formula gives the best approach for stars with ns = 3 and-n, = 4. 
In spite of this the primary energy for n, — 3 is overestimated by a factor 
of 2, and for n, =4 the overestimation is still about 30%. The deviation 
between the energy obtained in the first order approximation and the real 
energy can be expressed by the magnitude of C. The value of © for the stars 
with n,— 3 is 0.7 and for ns= 4 it becomes 0.9. At energies of 1012 eV 
SCHEIN et al. obtained C = 0.7 (82). 

The energy value obtained for n,—1 is wrong by orders of magnitude. 
This can be understood if one compares the angular distribution of the thin 
tracks from the stars with n, =1 with the calculated distribution as shown 
in Fig. 7. The two distributions look similar. This is possibly a chance co- 
incidence and does not necessarily mean that these stars mostly belong to 
events in which nucleon isobars only have been produced. Nevertheless this 
comparison serves well to illustrate the real situation. It shows that the 
stron peak in the forward direction belongs to « penetrating » protons. One aint 
expect that, on the average, the «other» proton was slow and produced a 


(#1) W. F. Fowrer, R. P. Suurr, A. M. THORNDIKE and 
Dre i W. L. W MORE : 
Phys. Rev., 95, 1026 (1954). WHITTEMORE: 


(2°) M. Scnery, D. M. Haskin, E. Lourmann and M. W. TEUCHER: Multiple Meson. 
Production in N-N Collision at 10!2 eV, preprint (Chicago, 1959) 
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gray track (see discussion in Section 3°1). These protons can not be distin- 
guished from the evaporation tracks and, using the « Castagnoli » formula, one 
is thus biased in favor of an overestimation of the energy. 


82.3. Frequency distribution of thin tracks. — In Fig. 2 we 
plotted separately the frequency distribution of stars with %,=1 a8 a function 
of their gray+black prong number. One can observe a shift to lower Ny? 
(This effect has already been found in early works (2%%)). For n,=1 the aver- 
age n, is %, — 7.2 while for the rest of the stars it is 10.5. This could be 
explained by means of two different effects. Once, we are biased if we take 
only the stars with n,=1 since, for stars with larger n, the probability of 
secondary interactions is larger and therefore more evaporation tracks are 
expected. Two, this shift can arise from the effect that, as stated above, the 
nucleon which was emitted backwards in the CM system had a low kinetic 
energy in the L system (especially in quasielastic collisions) and, if it has made 
an interaction, that should result in a lower excitation energy for stars with 
", = |: 

Concerning secondary interactions we want to make more remarks. As 
can be seen from Fig. 3 the stars with larger n, show, on the average, a 
larger n,. This must be due to secondary interactions for a twofold reason. 
First, if 7, is larger the chance for secondary interactions increases and, second, 
fast particles undergoing secondary interactions still have some chance to 
lead to additional thin tracks by producing other fast particles. That the 
latter effect really happens can be deduced from a comparison of the average 
ns of the inelastic p-free p collisions (n, = 2.06) with the total average number 
of thin tracks (which is 2.65). 


82.4. Frequency distribution of heavy tracks, the star size. 
— We want to say a few words about the distribution of gray+-black prongs. 
In a paper of CAMERINI et al. (3°) the star size distribution (represented in our 
notion by ”, and »,) has been studied as a function of the primary energy. 
These authors find no indication for a saturation in the average n, with in- 
creasing primary energy. Our average n, is 9.7 as compared with their 18 
(for 6 GeV). This is a remarkable difference which is probably due to an es- 
sential loss of small stars in their area scan. According to measurements of 
BIRNBAUM et al. (34), n, at higher energies is no longer a good measure for 
the primary energy. Our statistics are not sufficient to decide if there also 


(33) U. CAMERINI, J. H. DAVIES, P. H. FowLER, C. FRANZINETTI, H. MUIRHEAD, 
W. O. Lock, D. H. Perkins and G. YekutIsLI: Phil. Mag. 42, 1254 (1951). 
(34) M. Brrnpaum, M. M. SHapIRO, B. SmLLER and F. W. O’DELL: Phys. Rev., 


86, 86 (1952). 
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exists some kind of «knee» for our case of constant primary energy (see. fi. 
refs. (2) and (*)). In ref. (3°) a detailed discussion of the physical reason for 
this «knee » is given. Of course n, does not only depend on the primary energy 
but also on the atomic number, i.e. on the number of charged nucleons per 
nucleus available for evaporation and is therefore influenced by the compo- 
sition of the emulsion. 


82.5. Selection criteria, stars with small visible excitation 
energy. — In Table IIT we have demonstrated how the different separation 
steps, in order to separate the p-free p collisions, have influenced the number 
in the columns. Events with odd and even prong numbers were treated 
equally and use was made of the appearance of slow electrons or blobs: This 
procedure finally gave us the necessary additional information to deduce the 
real number of p-free p collisions among the possible ones. 

About 2 of the 33 possible inelastic p-free p collisions had to be attributed 
to p-bound nucleon collisions (Table ITI column N; and sub-Section 2°1.2). 
In this connection it is of interest to ask whether or not the fact that, after, 
strong selection, the sample left still has an indistinguishable contamination 
of p-bound n collisions can influence results which were derived for p-free p 
collisions only. This is of special interest for jets produced by ultra high 
enerey particles. There one generally does not know the primary energy and 
mainly the angular distributions of white stars serve to estimate the primary 
energy. In the following we shall discuss preferentially the white star or stars 
with small visible excitation energy. Certainly it will not be possible to com- 
pare our 6.2 GeV stars throughout with ultra high energy jets, but at least 
some points can be discussed qualitatively. 

The total number of white stars (excluding «scatterings » but including 
p-free p collisions) is 30, which makes about 5% of all observed stars. This is 
rather low as compared with the 10--15% for stars at ultra high energies (32). 
One probable explanation is that (in p-free p collisions at 6.2 GeV) some of 
the meson and one nucleon produce gray tracks, while at ultra high energies 
most tracks are thin. (Backward protons from N*N* formation can there 
still produce gray tracks, but .V*.V* formation becomes much more infrequent 
with increasing energy). In the following we shall briefly extrapolate from our 
Stars to cosmic ray jets. 

The fraction of white stars among the 33 possible p-free p inelastic colli- 
sions is about 0.4. If one believes this fraction to be about 1 at energies of 
10” eV one can extrapolate roughly (*) to higher energies with our 30 white 
stars. Assuming that there are enough cases, where the proton collided with 


(35) M. W. TEUCHER: Zeits. f. Naturfor., Ta, 61 (1952). 
(*) The influence of secondary interactions is neglected. here. 
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only one bound nucleon without giving visible excitation energy to the rest 
nucleus. The extrapolated number of white stars at 1012 eV is 30/0.4 = 75, 
which in fact makes about 11% and which is in agreement with the (10 —15)% 
obtained directly at that energy. Only about } of these (10--15)% belong to 
real inelastic p-free p collisions. The rest must belong to complex nuclei leading 
to odd and even prong stars, and one question is whether the angular distri- 
bution of the (thin) secondaries from these stars looks very different from the 
angular distribution obtained from p-free p collisions alone. The same may 
be asked of the average multiplicities and the reaction rates of the (thin) 
secondary tracks. Often the appearance of a white star as an even-prong 
event has been attributed to p-free p collisions, and for odd-prong events 
white stars are believed to show angular distributions and multiplicities still 
representative for p-free nucleon collisions, since one has concluded that the 
average excitation energy given to the rest-nucleus in white stars is small, 

The average thin track multiplicity of 33 possible p-free p collisions is 2.06, 
for the stars with small visible excitation energy (n, = 0, 1, 2) it is 2.09; and 
the total average is 2.65. Thus at first sight it seems that the stars with small 
visible exeitation energy do not show a considerable contamination with stars 
in which secondary interactions happened, but we believe that this, in reality, 
is not true. For the calculation of n, = 2.09 the 24 stars with 1 thin track 
only have been included (see sub-section 2°2.1). Therefore 2.09 is certainly 
not the right value to compare wity 2.06 from p-free p collisions. If we leave 
out the 24 critical stars we get 2.50. This indicates that already the stars 
with small visible excitation energy have a considerable contamination with 
stars in which secondary interactions happened. 

From multiplicity considerations alone we believe that even the white stars 
contain a non-negligible fraction of stars in which secondary interactions 
happened. On the other hand one can see from Figs. 4 that the angular distri 
butions of thin tracks for p-free p collisions, for the stars with small visible 
excitation energy and for all stars with n,<7 do not differ significantly. 
Therefore we can conclude that if one performs a strong selection (as in sub- 
Section 2°1.2) of possible p-free p collisions the selected sample is represent- 
ative of p-free p collisions as regards multiplicity and angular distribution. 
At higher energies it can happen that occasionally white stars occur in which 
secondary interactions were made. Such single events naturally are no longer 
representative for the correspokding pfree nucleon collision. 


We are highly indebted to Prof. SCHEIN for making the plates available 
to us and to Drs. HEARD, LOFGREN, TEUCHER, LOHRMANN and to the Be- 
vatron crew for providing us with this excellent exposure. 


3 - Il Nuoro Cimento. 
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RIASSUNTO (*) 


Si sono trovate 700 interazioni seguendo protoni di 6.2 GeV «lungo la traccia ». 
Si ottenne un percorso libero medio di interazione di (38.2+1.5) cm. Si separarono 
42 possibili collisioni protone-protone libero trovando che la sezione trasversale per 
la produzione di mesoni è 23*7 mb. Abbiamo trovato una molteplicità media per tracce 
di secondari carichi in collisione anelastiche p-p di 2.8 +0.9 ed il rapporto di branching 
fra eventi a 2 e 4 rami era circa 2. Questi valori vennero paragonati a quelli pre- 
detti dalla teoria statistica e da una sua modificazione. I risultati di queste ricerche 
sarebbero in accordo con un modello che suppone che a queste energie circa il 50% 
di tutte le collisioni sono «periferiche » portano solo alla formazione di isobare 
(di nucleoni) (NN* o N*N*). Si sono misurate le distribuzioni angolari di tutte le 
tracce secondarie sottili nel sistema L e si è applicata la formula di Castagnoli ef al. 
ad un campione di stelle ragionevolmente selezionato. Questo porta ad una stima per 


eccesso dell’energia primaria per un fattore di circa 2. Altri ricercatori hanno ottenuto 
a 10 eV un valore simile. 


(*) Traduzione a cura della Redazione, 
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On the Two Fluid Model in Multiple Production of Mesons. 


M. HAMAGUCHI 


Department of Physics, Kyoto University - Kyoto 


(ricevuto il 12 Marzo 1960) 


Summary. — The mechanism of multiple particle production in high 
energy nucleon-nucleon collisions is investigated using the two fluid 
model. The theoretical inelasticity is numerically calculated. 


1. — Introduction. 


The present note describes a picture of the dissipating system using the two 
fluid model (1), in which the core part of the incident nucleon behaves as the 
perfect fluid, while the cloud part surrounding the core behaves as the viscous 
fluid. We may imagine that the latter part is composed of 7-mesons, K-mesons 
and nucleon-antinucleon pairs, ete. At the first stage of the expansion, various 
fluids produced by high energy nucleon-nucleon collisions are enclosed in the 
flat-circular volume, where each constituent of these fluids is near thermal 
equilibrium at extremely high temperature. These fluids expand into the outer 
space with the development of time, keeping up the local uniformity of the 
energy. Since the core part does not contain the viscous effect, it may be 
easy to find its conspicuous distribution at the fronts of fluid motion. Although 
the cloud parts are lagging a little behind the core fluid in the expansion pro- 
cess because of the viscosity, they will also show the distribution concentrated 
intensely in the vicinity of the fronts of the expanding cloud in a one dimensional 
fluid motion. 

Here, we shal] remember the empirical model (?) proposed by Coccont, 
but not further discuss the connection between our model and Cocconi’s one. 


(1) R. P. Feynman: Phys. Rev., 91, 1291, 1301 (1953); 94, 262 (1954). 
G. Cocconi: Phys. Rev., 111, 1699 (1958). 
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To express the covariant energy-momentum tensor (5), we have used 


(1) TICO, rene re 


ke = i(ccre), (cloud) 


— ¥ 0,0” + 8, 0% + (org + v'O")} — Z Lg of}, 


j(cloud) j (cloud) 


where è4= 0,-0,0% g(a =1, 2,3)=1, g*=—1 and C=1. Y represents 
k=i,j 

the sum of each component of the core and the cloud specified by the indi- 

Cesu. 

In Section 2, we have solved the equation of motion 0,7,,— 0 in hydro- 
dynamical expansion process, where Reynolds number is very large #> 1. Ana- 
logously to Heisenberg-Landau’s (*) idea, we think that the use of a turbulent 
flow is most favorable for the hydrodynamic description of multiple particle 
production. In the domain considered, there is only a small fluctuation due 
to the viscosity. Therefore, an asymptotic expansion with respect to viscous 
coefficients will be used. 

On the other hand, the dependence on the temperature of the viscous coefti- 
cient should be directly taken into account, but we must keep in sight the all-over 
description of the phenomena in excess of our stress on this point only. So, 
we should turn our attention to the large Reynolds numbers mentioned above, 
even if there is a slight fluctuation of viscosity with temperature. 

In Section 3, we have obtained an expression of inelasticity using the 
(0, 0)-components of energy-momentum tensor at the critical time, and made 
out a few tables. As a consequence, the trend of the deviation of inelasticity 
in the two fluid model will be investigated. 


2. — Equation of motion. 


From the fundamental equation in relativistic hydrodynamics 8,7. = 0, 
. . . v pv 
we have the following simultaneous equations (*) of a one dimensional motion 


(5) G. A. KLUITENBERG: Relativistic Thermodynamics of Irreversible Processes (Leiden); 
E. C. G. STUECKELBER& and G. WANDERS: Helv. Phys. Acta, 26, 307 (1953); C. ECKART: 
Phys. Rev., 58, 919 (1940). 


(4) W. HEISENBERG: Zeits. f. Phys., 133, 65 (1952); L. D. Lanpau: ev. Akad. 
Nauk SSSR, 17, 51 (1953). 

() In the derivation of eq. (2), we have omitted the derivative of viscous coef- 
ficients with respect to space and time. The volume of the System at a given time is 
Va xa*x (a, linear dimension of system), so the energy doni e — Efra?x ~ T4. Here 
the T-dependence of viscosity being assumed as an example 4 toc T3, we can REST. 
automatically an x-dependence of viscosity n— D/xt, (D, a ra const). There- 


2244 


ON THE TWO FLUID MODEL IN MULTIPLE PRODUCTION OF MESONS 37 


in the direction of positive x-axis, 


40 1 de 0 ov ev 0 [ov dv]| 
= (A E 59= 120 + 02) — = I De 
siete Ia matsaltal * ae: 


+ y a fèfuè 120) 10/20 1 20] 
sen || Of 20) 20e 2 DEI? 
(2) 
ine aq oa . ; do 
a) ae. =(5) Se RA Li A te | 
3 0t 3 0& \v? j (ctoud) |a t OË| 20& |v? Ot è AI 
nf : (Ki 0 [dv 1 ov 1 0 [1 ov TL | 
scisso ' |2 Ctl Gt do de) 4 velo di — 208 GEIL? 


introducing the new variable co-ordinate É—=t—x. In the derivation of the 
above expression, we have imposed the ultra-relativistie condition upon the 
four velocity, 


(e) VU =0>1, Vo — 1 & 1/(20) , 15ESAÀ, 
where A is a Lorentz contracted width of the system. To solve the above 


simultaneous equations under the condition (¢), we make use of the method 
of asymptotic expansion with respect to n, È, 


o = (0) + FY mo) — Dd Lu}, 


j(cloud) j(cloud) 

e-0+ Y ne 3 Le 

Eri | li j cy ero) 
j(cloud) j(cloud) 


where higher order terms of 7, and ¢; were omitted due to the assumption of 
small viscosity fluctuation. At first sight, the propriety of application of asym- 
ptotic expansion is not always obvious for higher order derivatives of eq. (2). 
But, we will not meet with any difficulty, if only convergence of this series 


fore, the derivatives of momentum transfer in the equation of motion are estimated 


On [Ov n dv à [ov 1) 
D a VEN NS e 
ag a (4 — E) \0& 0& \ 0& E \ 0g 


as follows: 


In the domain of t > £> A, the former is negligible as compared with the latter. 
It will be also proved that in a general way this conclusion does not depend on the 


above special example. 
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and non-singularity are satisfied. This kind of method had been used by Heisen- 
berg in past years (5). 

The solution of the above equations can be carried out quite similarly after 
the method of Sections 1, 2, in I (9), in which at first a power function of 
(t/£) was assumed as the velocity, and after the simultaneous equations of € 
with respect to t, £ were solved. Loue the lengthy calculation gives à set of 
solutions, 


t UNE 
v2 a PE D rer (e) 


J (cloud) 


* o* oy £ LO a, el 
(3) (9) E RE exp | (ing ting ml in 5) i 


j(eloud) 


UNE Er 
an DI (183% — Ce") (=) ste > 5 Of, E À) + DI C;°O2(, cs À) 
S 1 J. 


4 * 0* 
SME E 
k =7 cloud),j(cloud) fed 


Both of the last two terms 0,, 0, of the energy density are negligible compared 
with the second term, therefore these additional terms will be omitted in the fol- 
lowing calculation. e°* is an initial value of the solution in the non-viscous equa- 
tion, and e!*, «?* are also those in the 1-st order equations with respect to 
viscosity. 7°, f; and f? are given by f° = 4V In (t/4)/ln (E/4), fio — 0.83e1*/e0*, 
ete. As is clearly seen from the region in which our solution is applicable, if we 
expediently assume f° & 1 and ¢;* & &* =— l280*, though a suitable method to 
determine the latter quantity does not exist, (J means the dimension of the length 


of order 1), fj & fj = 0.830 ef*/e* =y,<1 will be obtained. Then, the 
expression (3) is also written 


7 ic 
vt > pit (à) , 
iS; s 


3 (cloud) 
a È L In yi ! a È I 
1° Aa Va i ti 


i dEi eles 
“2% [rl (ni ems, fin Sn 1) tt) (E). 


(3) ~ Tel exp 


(5) W. HEISENBERG: Ann. d. Phys., 74, 577 (1924). 

(6) M. HamaGUucHI: Nuovo Cimento | 
quoted as I, II respectively. 

(*) «=4.54 and B=2.21 are the numerical values given in I. 


, 4, 1242 (1956); 5, 1622 (1957), hereafter 
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From this, we see that the energy densities of the cloud grow gradually during 
the expansion process due to the viscosity effect (n;< 0, È;> 0). 

Applying now the same method to this case as in the derivation of the 
limit of a one dimensional motion in I, we can also obtain the following limit 
using the velocity mentioned above, | 


2 g\2(a-1 
(4) nat 14 > vm ia) (5) |. 


la : 
Ss J (cloud) S 


Then, the expression of critical energies in the conical flight is also given by 


Di E tool ia 
= (in In, 7 int ‘In i) P(nx— 0x) (4) | , 


where k — à (core), j (cloud), and 7,;=¢;=0. Applying Vin (t,/A)/In (£/A) wy 2 
derived from f° ~ 1 to the above expression, we can rewrite it in a favorable 
form for the following discussion, 


ES -13 p\ B+4 
Ke a er (* D HE Cy) (:) È 


By inserting the limit (4) of a one dimensional motion into (5), we obtain the 
respective critical energy densities for the cloud particles (viscous) and the core 


(30) Ba (7) lex 


particles (non-viscous), 


as \4 9 A a\za=D È . E a\2B+8 
(6,) Epetbeay = ae (E) i I =F 3 , | VilN; PE G3) (2) | 1 ln; [æ)) (2) | ’ 
‘ ~ j(cloud) \ 


a A 8 al a 2(a~—1) 
(6) CHE a es" (, | L | 5 > AVE 6) (3 | i 


j(cloud) 


where #, is a critical time in free particle scattering. Here, we pay attention 
to 2(a—1) & 28 +8 and y;=7, 6; =. 

We can generally expect that the existence of viscous effect of the cloud 
results in an increase of the total entropy which is covered by the decrease 
of internal energy of the system. Here, assuming local equilibrium in the 
elements of mixing fluids, we may assume that this decrease of internal energy 
is covered also by the decrease of the energy density of non-viscous core parts. 
This fact differs from our concrete concept in treating of the energy density 
of the usual perfect fluid. A matter of primary concern for us is that the critical 
energy of the cloud shows an increase larger than that of the core in the 
course of the development of the system, due to the viscous effect of the last 


+ 
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term on the right hand side of (6,). Then, we have 


By & È = a sia 
(6") vie. eslet* À i a 


As we do not know the suitable quantitative estimation of 7, ¢, we cannot 
definitely settle the above inequality, but discuss about it to a certain extent 
in due consideration of the different critical temperatures of the core and the 
cloud (see Section 3). The inequality is also a mathematical expression which 
characterizes the properties of mixing fluids viewed from the two fluid model. 


3. — Hydrodynamic expression of the inelasticity. 


In this section, we calculate the mean energy densities of the core particles 
and of the cloud particles at the end stage of the collision process in the centre 
of mass system, and then express theoretically the inelasticity with these 
quantities. 

The (0, 0)-components of the energy-momentum tensor follow from the 
relation (1), 


7 oo _ 4 ae : Ov 3 Ov DI , OÙ 1 da 
(1°) k mol = AYA Et Te di To Abe Sk o ci . 


The velocity of the fluid element in the conical flight is also given by 
(7) ve t/a 


which is obviously valid for the mixing fluid (see eq. (47) in I). Substituting 
relation (7) into (1'), we obtain 


" 4 t\2 il t2\ 
a") ne GI 
e 3 fe > Me |= ee il aN 


using the volume viscosity €? = ¢,+2n,. 


Here, the mean values of 7° at the critical point t=? satisfying the re- 
lations (6,), (6,), are well defined by 


E fr fus si 


(*) See eqs. (41), (48) in I as to the integral domain. 
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Then, we have for the core and cloud particles, 


4 \2(a-1)-4 
Tiene) = 4(e7€ oft Tele een 
Le 9[3—2(a—1)] ict 6) À = Heiei )°, 
(8) Abd con & 4 ee")? Il 3 I IF a 2(a—1)—4 
(cloud) > (& OM oe ir PZ 6) 1 = 


Cei eee 

2[ [28 DEN le | | + Os(n, Le) = 4(e 65")? , 
respectively (*). As is easily verified, the viscous term 03(7, €) induced from 
last two terms of (1") should be neglected according to the relation 


Di a He e od = ET 
ln A ej > C-, (n wl). 


7, C” represent the mean values of viscous coefficient in the domain considered. 
Now, we shall define the theoretical inelasticity by a quantity Æ, which 
means the ratio between the energy of cloud particles and the energy of the total 
system, 
Ty, fel")! 
af eet ia Di \ Ss > (e € ) 


(9) 5 t(core) PIE i (core) 
K IE o CR 


j(clnud) i (loud) 


(*) The quantity e_;; defined by (8) is closely connected with ef. By straightforward 
calculation and omission of higher order terms of 7, £, we have 


c'y CA 1?(7 — ¢) Ve 
(ef 185) = (ee /¢4) È (2B +4) (i) | | 


If we fix together this expression and the inequality (6”), a relation (e{'/e%*) > (e/e?*) 
is also obtained. Remembering that the viscous contribution cannot affect the non- 
dissipative term, that is, Reynolds number #> 1, we can easily find 


Assuming the relation between energy-density and temperature e“ oc T4 which holds. 
in a way approximately analogous to black body radiation, we have 


1 
oe bP 


where 7° and T, are respective critical temperatures of the core and the cloud, and 
the initial stage of the total system is regarded as about 7' ~ co. 
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In our model, the cloud particles surrounding the core nucleon form the 
mixing fluids consisting of the 7-mesons and the K-mesons. The ratios of initial 
energy densities of these particle are given as follows: 


(eo lee) = (Gn PO) /9. 9 G(O)) & 0.857, 


(10) ig 
(ep lex) = (Fo P()/9.9 FO) & 1-148 . 


Here, we have used the relation 


Male 
ee a 


together with the statistical weights g, = 3, g,=—4 andg,»=gp=4. Let us 
now replace the energy in eq. (9) by the particle density through the medium 
of the relation (7) 


Sec 0 po.) 

i(core) — !li c~ 70 DI 
(10) as 
cf Jes A (7) av; 

€5 (cloud) — ny ui i P(x,) ;) 


0 


plc Mie NT Em cyt an dE TETE TN) SITA 
which 7°, 7°, are the critical temperatures of the core- and the cloud-particles, 
respectively. 

By inserting the relation (10), (10') into (9) we finally obtain an expression 


which after a simple, but lengthy, calculation may be written in the following 
form, 


Te (Poli 
(2 ee) i 
F (x9) 


(2n 5)? T° ai] 


Le ke 
ie 


(11) 


C\I o ) Cn) 3 È 
(Ma)? (0420 Ta) + (né)? (0571 


Here, we shall compute the inelasticity in the special case in which the cloud 
e erate pate lim ce and the core 
To Ti T.jV?2. The stage of Ta T #3 
in which the viscous effect 
critical energy e° 


temperature shows a change of 
corresponds to the extreme case 
of the cloud has a maximum contribution to 
; In that case our perturbational calculation may be destroyed. 


(7) Functions (9, F) and (9, #) were gi A 
, é 9, F) were given by Belenkii-L: -Milehi 
tati EA g y nkij-Landau-Milehin for the 


S. Z. BELEN’Kis and Jy. 9D. MU ) 
Nauk, 56, 309 (1955); S. Z. BeLEN’Krr : I Si PAR nd 
a EN KIJ and S. A. MILEHIN: Zurn. Esp. Teor. Fiz., 
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For the numerical values of F(a) and (x), we have used the Table given 
in reference (7). 

The following table gives the results of our numerical computation for 
the inelasticity. 


TABLE I. — K = inelasticity. 


| Composition of various particles | 
AE K ven + ee | Model 
np ni | nee | 
TOT, 0.73 0.04 0.35 Il | one fluid 
ZIE) 0.84 OOo 0.35 1 two fluid 


From Table I, we can understand the tendency to fluctuation of the inelasti 
city in the two fluid model. Æ value for she < F, is generally larger than 
CLONES 

In the case of T,<1.1m,c, the inelasticity increases gradually with the 
decrease of the critical temperature 7, (and also T,). We shall expressly point 
out the fact that this value attains the order of magnitude of KA ~ 0.9 at 
T, = 0.8m,¢?. 

If we apply the high critical temperature T,> 1.1m,c? to our model, it 
may be possible to make a forecast of the low value of inelasticity. In fact, 
at the critical temperature T° = T, =1.3m,c?, we could obtain K & 0.64 from 
eq. (11), after we made use of the appropriate quantities required for the 
calculation. In conclusion, we can expect a fluctuating value of inelasticity 
K =0.6~0.9 for the critical temperature 7’, = (0.8 ~ 1.3)m_c? considered. 


4. — Discussion. 


We have here discussed the mechanism of multiple particle production based 
upon the two fluid model and sought for the plausible expression of inelasti- 
city of this model. This expression is of the form containing the energy den- 
sities at the initial and the critical times. As is clearly known from (6°), if 


it were possible to suppose e5*/et = e;*/ef (i.e. n=C=0) we have a simple 


: . a OA i NET 
expression for the one fluid model, A = Siege > ele’, (es; critical 
j(eloud) j(cloud) 


energy of the total system). We cannot, however, approve the above assumption 
owing to the essence of two fluid model, in which the critical temperature of 
the core is not equal to that of the cloud. Then, we have had to use the rela- 
tions (9), (11) for calculation of inelasticity. That is to say, such a treatment 
gives enough ground for intervention of our two fluid model. 
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RIASSUNTO (*) 


Si analizza il meccanismo di produzione multipla di particelle nelle collisioni 
nucleone-nucleone di alta energia, usando il modello a due fluidi. Si calcola numerica- 


mente l’anelasticità teorica. 
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The Excited Nucleon’s Model of Meson Production i 
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(ricevuto il 12 Marzo 1960) 


Summary. — A model of meson production from excited nucleons is 
described for both nucleon-nucleon and nucleon-nucleus collisions induced 
by cosmic ray particles. The results are compared with observations 
in photographic emulsions, and reasonable agreement is found. 


1. — Introduction. 


A few years ago TOKAGI (') and KRAUSHAR and MARKS (?) suggested that 
at high energy nuclear collisions, mesons are produced from two centers. 
According to KRAUSHAR and MARKS (2), in nucleon-nucleon collisions two ex- 
cited nucleons are formed that live longer than nuclear reaction time; later 
on they decay independently into mesons and nucleons. Observations of LIn- 
DENBAUM and YUAN (*) at lower energies ((0.5--3.0) GeV) suggest a similar 
process of meson production in nucleon-nucleon collisions. At this range of 
energy one or two nucleons are excited to the first ($3) isobar level found in 
xt+-p scattering experiments. The isobars have a resonance-like mass distri- 
bution with a peak at 1.230 GeV and a life time longer than 10-**s. Sub- 
sequently each excited nucleon decays spontaneously into a nucleon and a 


pion. These observations were successfully extended by LinDENBAUM and 


(*) This work has been sponsored in part by the Geophysics Research Directorate 
of the Air Force Cambridge Research Center, Air Research and Development Command, 
United States Air Force under contract AF61(052)-58 through the European Office ARDC, 

(1) S. TAKAGI: Progr. Theor. Phys., 7, 123 (1952). 


(2) W. L. KrausHaaR and L. J. MARKS: Phys. Rev., 93, 326 (1954). 
(3) S. J. Linpenpaum and L. C. L. Yuan: Phys. Rev. 93, 917 and 1431 (1954); 


103, 404 (1956). 
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STERNHEIMER (4) and by HOLMQUIST (5) up to energies of few GeV, with the 
conclusion that mesons are produced predominantly from excited nucleons. 
The recent discovery of second and third resonance levels in photomeson- 
production (5) and x-p scattering experiments (7°) indicate the possibility of 
extending this theory to even higher energies. Further analysis of nuclear 
events produced by cosmic ray particles by Cocconi (°) and CIoK et ats (8) 
gave new evidence for meson emission from two centers. At the same time se- 
veral authors (1:11) observe that in these collisions mesons are produced with 
a low transverse momentum with a peak at about P,= 0.4 GeV/c. Both of 
these observations are compatible with a model of meson production from ex- 
cited nucleons. These recent developments make it worthwhile to explore the 
possibility of extending the excited nucleon’s model of meson production to 
cosmic ray energies. 

The purpose of this paper is an attempt to describe meson production by 
cosmic ray particles in terms of an excited nucleon’s model. We shall assume 
that in a nucleon-nucleon collision two excited nucleons are produced with 
masses M, and M, and momenta P, and P, given by an excitation function 
W(M, P, M, P,|E.) where E, is the total energy in C.M.S. of the collisions. The 
two excited nucleons have a life time longer than the collision time, and leave 
the reaction range as two free particles. Subsequently each excited nucleon 
decays independently into a number of low energy pions, and a single nucleon. 

A model of meson production in nucleon-nucleus collision will be described 
as follows: The energetic incident nucleon sees the nucleus as a gas of A in- 
dependent nucleons. The first phase of the reaction is a nucleon-nucleon col- 
lision between the incident nucleon and one nucleon of the nucleus. From this 


*) S. J. LINDENBAUM and R. M. STERNHEIMER: Phys. Rev., 105, 1874 (1957). 


) 
) F. N. Hozmquisr: Ph. D. Thesis (University of California, 1958). 
6) R. R. WiLson: Phys. Rev., 110, 1212 (1958). 

7) H.C. Burrowss, D. 0. Carpwezz, D. H. Friscu, D. A. Hitt and D. M. Ritson: 
Phys. Rev. Lett., 2, 119 (1959). 


(8) R. R. CRITTENDEN, J. H. SCANDRETT, W. D. SHEPARD and W. D. WALKER: 
Phys. Rev. Lett., 2, 121 (1959). 
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(E 
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( 


(9) Ss Cocconi: Phys. Rev., 111, 1699 (1958). 

(1°) P. CIok, T. CocHEN, J. GreruLA, R. HoxyxgsKI, A. JURAK, M. MIESOWICZ, 
T. SANIEWSKA and J. PERNEGR: Nuovo Cimento, 10, 741 (1958). 

(7) S. HasEGAWA, J. NISHIMURA and Y. NISHIMURA: Nuovo Cimento, 6, 979 (1957). 

(12) O. ae eee Y. NISHIMURA, M. Tsuzuxi, H. YAMANOUCHI, H. ra H. HASE- 
CAE NES ET AS MIO RUN AG A, Ye Fusimoro, § . HOSEGAWA, J. NIisHIMURA, K. Niv, 


K. NISHIKAWA, K. IMAEDA and M. Kazuno Curl Nuovo Cimento, 11, 125 (1959). 


(13) B. Epwarps, J.-Losry, D. H. PERKINS, K. PINKAU and J. REyNoLDS: Phil. 
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primary collision two excited nucleons are emitted in a narrow cone in the 
direction of the incident particle. We shall assume that the excited nucleons 
have the same interaction mean free path in nuclear matter as an ordinary 
nucleon. This implies that the mean collision time of an excited nucleon in 
nuclear matter is of the order of the mean reaction time, and is smaller than 
the mean life time of an excited nucleon. Thus most of the energetic excited 
nucleons, will make another collision before they will have a chance to decay 
into mesons. We further assume that one more excited nucleon is produced 
in the next excited nucleon-nucleon collision. In this manner a cascade of 
excited nucleons is developed in the nucleus. When the cascade of excited 
nucleons leaves the nucleus the excited nucleons will decay independently 
into mesons and nucleons in a narrow cone. This narrow cone of energetic 
particles is the jet of shower particles associated with high energy cosmie ray 
eventsin photographic emulsion. Some of the slow excited nucleons will decay 
into pions and nucleons in the nucleus. These particles in further collisions 
transfer energy to the entire nucleus and produce the large star with grey 
and heavy prongs observed in nuclear emulsion. 


2. — Meson production in nucleon-nucleon collisions. 


According to the excited nucleon model, meson production is a two phase 
process: 1) creation of excited nucleons in nuclear reaction; and 2) decay 
of excited nucleons into mesons and nucleons. Both of these reactions are 
strong, and conserve among other things, isospin, strangeness and parity. In the 
most general case the excited nucleons and subsequently the mesons, may carry 
different quantum numbers of charge, angular momentum, isospin and strange- 
ness. For the moment we shall ignore this specification and deal only with 
the simple neutral nucleons and pions. Later on, the theory may be modified 
to include the different charge and strangeness states. 

Following LINDENBAUM and STERNHEIMER (4) the relative probability for 
creating two excited nucleons of masses M, and M, in a nucleon nucleon col- 
lision is: IH |*(|P;| HB. /(£,+ E.)), where H is the matrix element for this 
reaction, E,P,; E,P, are the energies and momenta of the two excited nucleons 
in C.MS. In this system E,+E,= E. the total energy of the system and 
|P.|=|P.|. At low energies | H(M,M,) le has a resonance found in 7*-p scat- 
tering experiments. At higher energies there are indications for second and 
third resonances corresponding to isobars of masses M —1.750 and M —2.071 GeV 
respectively (55). 

The behaviour of |H | at even higher energies is not yet known, and in 
order to apply the isobar model to cosmic ray energies we shall have to 
guess the general character of |H FE Later on by comparing the results of 
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that model with experiment we can find out how well the average features 
of |H|° were pictured, and in this way learn something about its behaviour 
at higher energies. 

The most simple assumption is that at high energies there are many isobar 
levels of the same strength and thus |H | averaged over the different reso- 
nance levels is a constant. This assumption yields a meson multiplicity that 
rises like H, with energy (15). Experimentally we know that meson multiplicity 
increases very slowly with energy, and thus ‘|H°|) must be a decreasing 
function of the mass of the isobar. As a second attempt let us assume that 
in the energy range (1--1600) GeV |H(M,M,) |* averaged over all resonance levels 
could be approximated by exp [—(M,+ M,)/x], where « is a free parameter 
to be adjusted so that the calculated multiplicity in the energy range 
(25-1600) GeV will be as near as possible to the observed one. Furthermore 
we shall assume that the excited nucleons are produced isotropically in the 
C.M.S. of the collision. 

For the decay of excited nucleons, we shall assume that mesons are emitted 
from excited nucleons, one after another isotropicaliy and with a constant 
low energy Æ,. A choice of H, between 0.3 and 0.5 GeV is suggested by 
several authors (141617), As we shall see later, this type of decay yields in 
the ©.M.S. a narrow distribution for the transverse momenta of pions with 
a peak at PIN Em. 


3. — Monte-Carlo calculation of meson production in nucleon-nucleon collisions. 


In order to study the excited nucleon’s model of meson production in 
nucleon-nucleon collisions, a Monte-Carlo calculation was programmed for the 
WEIZAC electronic computer. The calculation was carried out for a given pri- 
mary E. At this energy a number of collisions were followed. For each col- 
lision two excited nucleons of masses M, and M, were chosen according to the 
weighting function W(M;M,) = exp [—(M,+ M,)/«](| P|’ E,E,)/(E,+ E,)), where 
E; = P\+M; and E;= P?+M? are the energies of the excited nucleons in 
their C.M.S.; E\+ E.= E, is the total energy of the collision in that system 
and |P,|=|P,|. Both M, and M, are heavier than 1.230 GeV, the mass of 
the lowest excited nucleon that decays into a meson, the (3 $) isobar. « is a 
free parameter to be found by comparison with experiment. The direction of 
each excited nucleon P, (—— P,) 


in the C.M.S. was chosen isotropically. Sub- 


co) WE Maor and G. YEKUTIELI: T.S.R. no. 1 (unpublished, 1959). 
(18) S. KANEKo and M. Oxazaxr: Nuovo Cimento, 8, 521 (1958) 
(7) E. LOHRMANN: Nuovo Cimento, 5, 1074 (1957) 
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sequently each excited nucleon emits one after another n, mesons. Each meson 
is emitted isotropically and with a constant energy E_ from the excited nucleon. 
Whenever a meson is emitted the excited nucleon recoils to conserve momen- 
tum. The last meson is emitted with a momentum chosen to insure that 
the last excited nucleon will go into a nucleon. Calculation of each collision 
gives the momenta in the Laboratory System of » created mesons and two 
outgoing nucleons. This output is used to study the following features of 
nucleon-nucleon collisions: 1) Meson multiplicity; 2) Energy distributions of 
mesons and nucleons, and 3) Transverse momentum (P,) distribution of me- 
sons and nucleons. 

In order to find the values of the two free parameters of this model, « 
and Æ,, several small samples with different values of x and E, were tried 
out at seven different incident energies: E, = 25; 50; 100; 200; 400; 800; and 
1600 GeV. The average numbers of mesons obtained in this way are pre- 
sented in Table I. The P, distributions of the different samples are all very 


TapLe I. — Meson production (n_) in nucleon-nucleon collision for different values 
of E_ and a. 


Free parameters Incident energies in GeV 
a Br ope DM a MS whe Aco 800 1600 
(GeV) (GeV) ris 
3.0 .375 9 11.2 SQUAT 
9.0 SSIS) 8.6 10.6 18.7 Mies 18.6 
3.0 .400 (*) TEL 10.7 13.4 14.9 | 16.5 Ted ee 
4.0 Sells | | 28.3 
4.0 .400 | 18.3 21.0 DDS 20.8 17.6 
All samples N = 10. () Ir = 250 


similar. The values of the transverse momenta are mostly below 1 GeV/c with 
a pronounced peak near 0.350 GeV/c. The peak of the P, distributions in the 
different cases varies slowly between (320-380) MeV/c and it is in reasonable 
agreement with the experimental observations. Therefore the average meson 
multiplicity in this model depends mainly on the choice of «a, i.e. Da the 
average mass of the excited nucleons. The only experimental information on 
meson production in this range of energy is from observations in nuclear emul- 
sions. In using the nuclear emulsion data we made the assumption that the 
group of stars with NA <2 represent mainly nucleon-nucleon ame The 
energy of the incident particle is measured by the median angle and the CASTA- 
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GNOLI (18) methods. The average numbers of shower particles ns of nuclear 
events with Nh <2 as observed by CASTAGNOLI et al. (18), EDWARDS et al. (15), 
GUREVICH et al. (19), LOHRMANN (17) and BRISBOUT et al. (2°) are presented in 
Table IT. 


Tage Il. — Meson multiplicity nz as function of incident energy E in stars with N,< 2. 
is average number of shower particles and n is the number of observed events. 


n 


s 


È Nucleon-nucleon collisions 
Emulsion data for x=3.5 GeV and E_=.400 GeV 
Energy range | (E,> | : 
(GeV) (GeV) Ns Na n ns re 
25— 150 70 T8 Os | 12 9.9 13.4 
150— 300 250 12.5 17.3 2 10.9 14.9 
300— 600 560 10.0 156 1 1220 16.5 
6001 200 905 15.5 DIL | 2 12.5 Wes 
1200—2 400 | 1050 11 16.1 | 7 1227 ITS 
Mx 000 1277 8.6 11.4 185) 02 13.8 
6002400 1790 12.5 le? | 9 12 LES 


The average number of mesons per event is given by n_ =1.5(n,—1) on 
the assumption that all mesons are pions. By comparing Table I to II it is 
evident that the best choice of the free parameters is x=3.5 GeV, E_=0.4 GeV. 

Using these values an extensive calculation of 250 collisions at each of 
the seven incident energies was carried out. The results of these nucleon- 


nucleon collisions are preseted with those of the nucleon-nucleus collisions 
in Section 6, 7 and 8. 


4. — Nucleon-nucleus interaction. 


Igt 
In Section 2 a model for nucleon-nucleon collision is described, we shall 
extend this model to nucleon-nucleus interactions, 


18 È NO 
(#5) Sì CASTAGNOLI, G. CORTINI, C..FRANZINETTI, A. MANFREDINI and D. MORENO: 
Nuovo Cimento, 10, 1539 (1953). 


(19) I. .I. GuREWICH, A. P. Micxacova, B. 
Journ. Exp. Theor. Phys., 34 (7), 185 (1958). 


.. (°°) F. A. BRrISBOUT, C. DAHANAYAKE, A. ENGLER, Y. Fusimoro and D. H. PERKINS: 
Phil. Mag., 1, 605 (1956). 
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The nucleus is pictured as a gas of free nucleons with a spherical sym- 
metric distribution o(r). Observations by HOFSTADTER (21) and collaborators 
on electron nucleus scattering, show that heavy and medium nuclei can be 
«described by a core of constant density 0, and a periphery around it in which 
‘the density decreases gradually. This type of nucleus can be approximated 
by a trapezoidal density distribution: o(r < R,) =- 0, = constant, where R, — 
=— 1.56+1.08A*is the radius of the core measured in units of 10-13 em; 
and A is the number of nucleons in the nucleus, and 0(R, <r< R,) = 
= 0,(r — R;)/(R;— R,) in the periphery. The width of the periphery is constant 
R,— KR, = 3.12-10-% em. For instance a nucleus of an AgBr « Atom » in photo- 
graphic emulsion will have A = 94; R, = 3.45-10- em and o = 1.74:10%5 em-3 

The incident nucleon is of very high energy (above 20 GeV) and we shall 
assume that it interacts singly and independently with nucleons of the nucleus, 
with the same cross-section for nuclear interaction as in nucleon-nucleon col- 
lisions. According to BRENNER and WILLIAMS (22) the cross-section for nucleon- 
nucleon interaction between 20 and 50 GeV is 21 mb. We shall use this value 
for nucleon-nucleon collisions up to 1600 GeV. This cross-section corresponds 
to a mean free path for nuclear collision of À = 4.7-10-13 cm in the core of 
an AgBr nucleus. 

The first collision of the incident nucleon in the nucleus is a nucleon-nucleon 
interaction and is treated according to the model described in Section 2. In 
this collision two excited nucleons are produced. We shall assume that each 
excited nucleon has the same interaction mean free path in nuclear matter 
as ordinary nucleons. An excited nucleon has a mean life time longer than 
a nucleon-nucleon reaction time, and thus the mean time for nuclear collision 
in a nucleus is about 1.5 times longer than its reaction mean time (it is pro- 
portional to the ratio: «mean free path for collision in nuclear matter » over 
«range for nuclear collision »). 

Now, the mean life time in the system of the nucleus of a very energetic 
excited nucleon is increased due to its Jarge Lorentz factor y, and it is cer- 
tainly longer than its mean time for collision in the nucleus. Thus an ener- 
getic excited nucleon with a y factor larger than a critical value y, on the 
average makes a collision in the nucleus before it has a chance to decay. On 
the other hand a low energy excited nucleon (with y < 7) will on the average 
decay inside the nucleus. The critical value of y, depends on the ratio of the 
mean life time of an excited nucleon to its mean time for nuclear collision in 
the nucleus. Considering the assumptions made before on the nuclear inter- 
action of excited nucleons we shall take y, = 1.4 for the critical Lorentz factor 


(21) R. HOFSTADTER: fev. Mod. Phys., 28, 214 (1956). 
(22) A. E. BRENNER and R. W. WiLLiaMms: Phys. Rev., 106, 1020 (1957). 


2259 


52 U. MAOR and G. YEKUTIELI 

of excited nucleons. This is the first condition for the development of à cascade 
of excited nucleons in the nucleus. Only excited nucleons with y < yo will 
participate in this cascade. We shall further assume that in excited-nucleon 
nucleon. collisions, as in nucleon-nucleon collisions two excited nucleons M, 
and M, are produced with relative frequency W(M,M,) exp [— (M,+ M,)/«]- 
-(| P, | £,2,/(£,+ E,)) with the condition that the masses of each secondary 
excited nucleon is not smaller than that of the primary excited nucleon. To 
complete the conditions for the development of a cascade of excited nucleons 
in the nucleus, we assume that a target nucleon in the nucleus may interact 
only with one excited nucleon and that there is no mutual interaction between 
two excited nucleons. According to this model the incident nucleon will ini- 
tiate a cascade of excited nucleons in the nucleus. Most of the excited nucleons 
(with y> 1.4) will leave the nucleus and decay outside the nucleus into pions 
and ordinary nucleons. The decay of these particles will follow the same pro- 
cess as described for excited nucleons produced in nucleon-nucleon collisions. 
(See Section 2). These particles will constitute the jets observed in nuclear 
events produced by cosmic ray particles in photographic emulsions. Some 
of the excited nucleons (with y < 1.4) will decay in the nucleus and give up 
their energy to excite the nucleus and to produce the cosmic ray star asso- 
ciated with high energy jets of cosmic ray events in nuclear emulsion. 


5. — Monte-Carlo calculation of nucleon-nucleus interaction. 


A Monte-Carlo calculation of the nucleon nucleus interaction was pro- 
grammed for the weErzac electronic computer. The block diagram of this 
program is shown in Fig, 1. The actual calculations were carried out for two 
nuclei: a heavy nucleus with A = 94 representing the nuclei of Ag and Br 
atoms in the emulsion and a light nucleus with À —14 representing the nuclei 
of the atoms ©, N and O in the emulsion. For each nucleus the reaction was 
studied at five different incident energies: EH = 100; 200; 400; 800 and 1600 GeV, 
and for every incident energy 100 cascades were followed. Tiron a the calcu- 
lation as in the nucleon-nucleon collisions the parameters x-3.5 GeV and E_ 
=0.4GeV were used. The output of this numerical calculation is used to 
Be ice features of nucleon-nucleus interaction. 1) Meson mul- 
Iphicity,; 2) Energy distribution of mesons an 3; vers 
momenta distribution of mesons and D The ed 
way are presented in the next Section. 
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1. Input dato: energy of incident nu- 12. If all cascades were followed, energy 
cleon. Size of target nucleus. Num- spectra, transverse momenta distri- 
ber of cascades to be followed. | bution and multiplicity of mesons 
SE nes : ae | and nucleons are calculated and 
| printed out. End. 
2. Number of completed cascades is | 
checked. | | 
3. If not all cascades were followed, 
a new cascade is started. A point 
of entry to nucleus is chosen. | 
LS ae | Sas SÈ i 
: ERA | = È 
The Loc n of nucleons passed | | 11. If cascade is completed, numbers, 
before collision takes place is chosen. | transverse momenta and energies of |. _ 
— | produced nucleons and mesons are 
| | stored. 
5. The nucleon’s density in a tube with | <= . 
a cross section o=21mb along the fi 
particle’s path in the nucleus is ra # a 7 
integrated until mn nucleons are | 6c. If it is an incident nucleon. Event 
counted, or the particle leaves the | is numbered to calculate transpar- 
nucleus before that. Nucleons that | ency of nucleus. 
were knocked out of this tube in SS 3 === = 
previous collisions are not counted. 
ae Fa vs | | 
| $ | È 
6. At the position along the path where 6a. Particle leaves nucleus without col- 
n nucleons are counted a collision lision. Identity of particle is checked. 
takes place. ar — 
| ; af Doe ae | 
7. Masses and momenta of two excited 6b. If en. particle decays into mesons 
nucleons (e.n.’s) are chosen. y of and nucleons. Momenta of decay Fa 
en. is compared to yo. particles calculated and stored. 
| 
} 
8. Fast en. (y > yo) stored for next Lee) 8a. Slow en. (y < yo) decays and gives 
generation of cascade. up energy to excited nucleus. 
Ÿ Jedi 
9. Store of e.n.’s is checked. = a 
1 
| = 
; | 
10. If store is not empty new e.n. is | 10a. If empty, cascade ends. STIA 
taken out of it to continue cascade. 
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Fig. 1. — A block diagram of a cascade of excited nucleons in nucleon-nucleus 
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6. — Meson production in nucleon-nucleon and nucleon-nucleus collisions. 


The average numbers of mesons #, and shower particles , produced in 
nucleon-nucleon and nucleon-nucleus collisions as function of the incident 
energy are given in Tables II and III. Assuming that all mesons are pions, 
the average number of shower particles », is given by: ny = 0.57 y+ 0.66”, 
where # y is the average number of relativistic nucleons produced per inter- 
action. If there are 20% K-particles among mesons n, is reduced by 4%. 

The main features of the cascade of excited nucleons in AgBr and CNO 
«nuclei» are summarized in Table III 


Taste IIT. — Average numbers of energetic nucleons (Ny) and mesons (n,) in nucleon 
nucleus collisions. 


AgBr = 21) mib AgBr o = 40 mb 
E (GeV) Ny Nr |W ye N y Na nr/N yp 
100 4.8 18.9 6.0 Told 38.9 5.5 
200 Dali 38.6 7.6 8.11 54.9 6.8 
400 ee 46.7 9.0 8.17 69.5 8.2 
800 | 5.1 50.4 9.9 8.83 85.2 9.7 
1600 | 5.6 61.6 11.0 9.59 103.0 10.7 
| = 
0.04 [H] = 0.16 [CNO] + 
CNO = PA ida + 0.80 [AgBr] 
(oy == PAL tomo} 
| 
100 | al 19.6 6.2 4.4 26.8 6.1 
200 | 50 25.0 Ta 4.7 35.4 Lod 
400 | 3.4 | 30.5 9.2 4.7 43.0 9.0 
800 | 3.6 35.5 9.8 4.7 46.7 99 
1600 | 059 34.5 10.5 all 55.4 10.9 
|- ele È x 


Results are given for four different types of cascades: a) A cascade of 


excited nucleons in heavy nuclei of the emulsion, the AgBr group, with 


nt) AVA ae È 
o=21 mb for nucleon-nucleon collisions. This value for the cross-section is 


the one measure NNER « 
sured by BRENNER and WILLIAMS (22) at the highest energy meas- 


ured so far, 50 GeV. 6) The same cascade but with o = 40 mb to study the 


effect of the nucleon-nucleon interaction cross-section on the nucleon nucleus 


collisions. ¢) An excited nucleon’s cascade in a light nucleus of the CNO group 


} | En) ' A Says) : i 
with o = 21 mb, and d) The results of a cascade of an average nucleon nucleus 
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interaction in photographie emulsion. The average interaction is given by 
0.04 [H] + 0.16 [CNO]+ 0.80 [AgBr], the relative contributions of the three 
types of collisions in emulsions were calculated as follows: Oy». = 21 mb for 
the nucleon-nucleon collisions and the nucleon-nucleus cross-sections obtained 
from the Monte-Carlo calculation, for nucleon-CNO and nucleon-AgBr col- 
Jisions. 

The average number of fast nucleons produced per interaction N,, is equal 
to the average number of excited nucleons that participated in the cascads. 
It appears from Table III that N y depends very much on the size of the 
nucleus and the cross-section for nucleon-nucleon collision of the cascade. At 
the same time the number of fast nucleons in the same nucleus increases slowly 
with increasing incident energy. The average number of pions produced per 
excited nucleon (#_/N x) behaves the opposite way. It depends only on the 
incident energy £,, and it is practically constant for different types of cascades 
induced by the same primary energy. Thus the total number of mesons pro- 
duced in a given cascade depends on the primary energy through (n IN x) 
(the average number of pions produced per excited nucleon), on the size of the 
nucleus and on the nucleon-nucleon collision cross-section, through (N \;) (the 
average number of nucleons produced per cascade). 


7. — Observable features of nucleon-nucleus interactions. 


Unfortunately we cannot compare directly the results of the excited nu- 
cleon’s model with experimental observations. Most of the experimental ob- 
servations between 100 and 1600 GeV were made in photographic emulsion 
exposed to the cosmic radiation. In these experiments it is impossible to 
measure directly the energy of the primary and the size of the target nucleus 
of each individual event. 

Several methods (215) were used to estimate the primary energy of high 
energy nuclear events in photographic emulsion. They are based on the as- 
sumption that all mesons are produced in a single center from which they are 
emitted with a forward-backward symmetry. According to CASTAGNOLI ¢ 
al. (18) the Lorentz factor of the emitting center y,, is related to the angular 
distribution of the shower particles in the laboratory system by: 


my. = nz? Yn cotf = <ncot 9, 


where y? = (1— f?)-! and B, is the velocity of the emitting center. 


(283) N. M. Durrer and W. D. WALKER: Phys. Rev., 23, 215 (1954). 
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Equivalent results are obtained by the median angle and the #-plot (?*) 
methods. The <Incot 0) parameter was found very useful for estimating 
the primary energy £, in nucleon-nucleon collisions, where Ey/m, = yo = 2y2—1. 
The same method was applied to very high energy nucleon-nucleus colli- 
sions (2425), on the assumption that all the N struck nucleons of the target 
nucleus form a single cluster, moving with a Lorentz factor y,, from which 
mesons are emitted with a forward-backward symmetry. In this case again 
Iny, = <Incot > and the primary energy E, is given by: 


Em, y = N(y —1) +y VN? (y2 — 1) By ; 


Thus the estimation of the primary energy Æ, in the nucleon-nucleus inter- 
actions depends on an additional assumption—the formation of a single cluster 
of struck nucleons, and on an extra parameter, N, the number of struck 
nucleons in the target nucleus. 

According to the excited nucleon’s model mesons are emitted from many 
different centers, thus one may think that the relation In y. — <Incot @) is 
meaningless for this model. However as In cot 9) can be evaluated for each 
event, we may use it to define a fictitious system of emission with Lorentz 
factor y.. The parameter y. so defined is treated as a stochastic variable, and 
an attempt is made to verify statistically the relation 


Voter (Yo de N)(1 + 2yoN == nN?) 
for different values of primary energy 
E=MVo ; 


where N is the number of the excited nucleons per interaction. 

To this end the distributions of Iny, = {In cot0> as calculated by the 
Monte-Carlo pregram are studied for nucleon-nucleon, nucleon-CNO and 
nucleon-AgBr collisions at different incident energies E,, and fixed number 
of struck nucleons N. 

Some examples of <In cot 0) distributions evaluated by the excited nu- 
cleon’s model for H, CNO and AgBr are shown on Fig. 2 and 3. The distri- 
butions of the nucleon-nucleon collisions have pronounced maxima and can 
be approximated by Gaussian curves with an averaged dispersion o = 0.2. 
The mean values of In y. for these distributions are summarized in Table IV. 


(24) G. Cooconr: Phys. Rev., 98, 1107 (1954). 
(25) E. M. FRIEDLANDER: Nuovo Cimento, 14, 976 (1959) 


2264 


THE EXCITED NUCLEON’S MODEL OF MESON PRODUCTION oT 


In the same table the expected y,’s ( ny (yo+1)/2) are compared with 
the corresponding mean values 


TABLE IV. — Comparison between expected and calculated In y y, as function of incident energy. 


Incident Expected The mean In y, of Te tous calculated 
energy E, In y, | the distribution |”, expected 
| 
25 1.32 1.46 1.16 
50 1.65 1.74 1.09 
100 1.99 2.09 1.08 
200 2.34 2.40 LO 
400 2.68 2.72 1.04 
800 3.03 3.05 | 1.02 
1600 3.48 5.40 1.02 
30 T si 


£,= 50 GeV E,=100 GeV 


i Ce 
100 150 200 SCA > 150 200 250 300 


legal 
E=400GeV + 


n cot 8 


SATA] ST 
200 250 300 359, co 
n CISA 


E, =800 GeV 


lee TÈ 
250 300 350 400 a 300 350 400 450 
i i An cot6) 


Fig. 2. — <In cot 6> in nucleon-nucleon collisions. 


It appears that for the excited nucleon’s model the relation In(Cy,) = 
= (In cot 0) holds, where is a small correction factor that varies from 1.16 
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to 1.02 between 25 and 1600 GeV. Thus the <In cot 0) parameter may be 
used to estimate the incident energy in nucleon-nucleon collisions. 


T ml 
E,=100 GeV CNO RESSE GeV 
201 13 2 1 | | 
% an! 
10+- 


T T et T 
100 150 00 : 00 TSO 2 
5 2 2% Wee” (SO 200 “elt 
| T oi 
E= 400 GeV CNO 
5 4 2 
20 
id 4 
% 
10 
Ts 
So, 200, 250 00 150 COCO 20) S00) 
(ln cot 8) 
| T gua yl 
E,=1600 GeV E,=1600 GeV A 
8 
207 i 
% 
10 
200 250 300 350 200 
{ 250 00 
(n cote) ù 
Fig. 3. — <In cot 0) in nucleon-CNO and nucleon-AgBr collisions. 


The Gn cot 0) distributions of nucleon-nucleus collisions are more compli- 
cated Each distribution is a superposition of several distributions of colli- 
sions with one, two or N struck nucleons. The arrows above each distribution 
show the expected In y, for different numbers, N, of struck SEE as cal- 
culated by the relation: y, = (yo + N)(1 + 2y, N+ N°). The corresponding 
maxima suggest that with good statistics the Gn cot 0) distribution can be 
used to estimate both the number of struck nucleons antl fie incident ener 
in nucleon-nucleus collisions (compare FRIEDLANDER (22) 4 

The average values of (In cot 0) for cascades with the same number of 
struck nucleons N are plotted on Fig. 4. On the same figure the et 
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In y, according to relation Ve = (Yot N)(1+ 2y)N+ N°) is shown by full line. 
The good agreement demonstrated by Fig. 4 verifies the relation In Ve 
= {ncotO> for nucleon- 
nucleus collisions as well. 
Finally in order to estimate 
the dispersion of «In cot 0) 


E =200 GeV 


in nucleon-nucleus collisions, SÙ Crete 
several (In cot > histograms 219 Se == 
are plotted on Fig. 5 for Ens 
cascades of fixed N, at 
E, = 400 Gov. 20 

In spite of the small 15 
statistics the dispersion of 1.0 


le 2-3-4 5-5 6=7-8=9 EN Et) 
the nucleon-N nucleons col- 


lisions is of the same order Fig. 4. — {In cot 0) as function of number of struck 
; o nucleons. Targets: eH, o CNO, x AgBr. 
as that obtained for nucleon- 5 ? 8 
nucleon collisions. 
The results obtained with the help of the excited nucleons model for both 
nucleon-nucleon and nucleon-nucleus collisions show that the verification of 


the relation In y, = <In cot 0) does not prove meson emission from a single 


coe 
À 
20. (28 
{In cor 8) 
c | di J | | 
Sea N=4 N=5 
Li je % 
È i { 
T in | 
19 20 ds n 20 259 iS 2.0 2.5 
In cot 8) 
Fig. 5. — <Incot6> distribution for cascades of different N in nucleon-CNO and 


nucleon-AgBr collisions at £,=400 GeV. 


center. On the other hand we believe that for a number of models similar 
to the excited nucleon’s model this kinematical relation will hold. It expresses 
the fact that for these models, mesons are emitted from a UMN of centers 
fluctuating around a fictitious system moving with a Lorentz factor y. = 


= exp (In cot 0). 
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The above discussion shows that the Castagnoli parameter <In cot 0) can 
be used to estimate the incident energy in nucleon-nucleus colisions only when 
the number of struck nucleons N is knewn. According to the excited nucleon’s 
model the number of struck nucleons in a nucleon-nucleus collision depends 
mainly on the nucleon-nucleon interaction cross-section and very little on the 
incident energy. The frequency of occurrence of struck nucleons in nucleon-CNO 
and nucleon-AgBr collisions, the transparencies of these nuclei, and their re- 


lative interaction cross-section in 


emulsion, are given in Table V. 


Tage V. — Frequency of occurrence (%) of struck nucleons in nucleon-CNO 


and nucleon-AgBr 


collisions for © y = 21 mb. 


Relative occurrence (%) of N struck nucleons Nuclear| Relative 
T t trans- cross 
tea ù parency |section in 
1 2 3 4 5 | 6 | 7 8 | 9 (%) | emulsion 
a | ss : 
CNO | 41 24 13 10 5 4 2 1 0 80 0.16 
AgBr| 22 | 14 | op i | ible 7 5 8 46 0.80 
fiat an 
N25 
(30) 
[Relative frequenc 
in emulsion | 
NES AE 
[19] 
<\nco/ e> 
1.5 2.0 25 3.0 35 4.0 
Fig. 6. — Average multiplicity <n.) versi 
g ipheity <n,) versus {In cot 6) for cascades of fixed Nat dif: 


ferent i 


neident energies daro 
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Meson multiplicity in nucleon-nucleus collisions depends strongly on the 
number of nucleons, N, struck in each reaction. In order to study the cor- 
relation between shower multiplicity, n, and number of struck nucleons, N, 
we plot in Fig. 6 average n, versus average (In cot 9) for cascades of the same 
incident energy E, and with the same number of struck nucleons N. Fig. 6 
summarizes the average observable features of the excited nucleon’s model. 
Each point on this plot represents a nuclear reaction given by its multiplicity 
n, and its «effective logarithmic incident energy »: In Yo = {In cot 0). The ver- 
tical curves are lines of constant primary energy £,, and the horizontal curves 
are lines of constant number of struck nucleons per collision, N. For the 
purpose of analysing nuclear events produced in photographic emulsions by 
cosmic radiation, the relative frequencies of collisions in photographic emul- 
sions with N struck nucleons, and the relative primary cosmie ray inten- 
sities are also given in brackets. If the fluctuations around their means of 
both n, and In cot 0) were small, we could use Fig. 6 to find the incident 
energy (H,) and the size of the nuclear cascade (N) of individual nuclear events 
in photographic emulsion. In reality the fluctuations of n, and «In cot 0) are 
not small and do mix high energy large size cascades with low energy small 
size cascades, and thus only the average features of cosmic ray events can 
be analysed. 


8. — Comparison with cosmic ray events in photographic emulsion. 


Several authors observed high energy cosmic ray events in photographic 
emulsion. Each event was analyzed in terms of: 1) its charged multiplicity ,, 
2) the Lorentz factor of the emission center y, evaluated by the median angle 
or the Castagnoli e¢ al. method, and 3) its star size N,. The first two obser- 
vables can be compared directly with the similar parameters of the excited 
nucleon’s model. The relation between the star size N, and the number of 
struck nucleons N is more complicated. Nuclear events observed in photo- 
graphic emulsion are divided into three groups according to their star size N,: 
1) nucleon-nucleon collisions including N = 1 reactions in CNO and AgBr stars 
with N, <2; 2) nucleon-CNO collisions with probably some low excitation 
nucleon-AgBr collisions excluding nucleon-nucleon collisions in CNO and AgBr; 
—stars with 2< N,<7, and 3) all other nucleon-AgBr collisions—stars with 
DERE 

For comparison with the excited nucleon’s model, the N,<2 group is 
equivalent to cascades with N =1 i.e. nucleon-nucleon collisions. All nu- 
cleon-CNO reactions with N>1 are included in the second group of stars 
with 2< N, <7. Most of the nucleon-AgBr reaction with N>1 belong to 
the third group. As we do not know how to distinguish a nucleon-AgBr col- 
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lision among events with 2< N, <7, we shall tentatively include all nucleon- 
AgBr reactions with N>1 in the group of stars with N,> 7. 

On Fig. 7 the observed cosmic ray events are plotted according to their 
multiplicity n,, (in cot 0) and star size N,. The three curves on the same 


25 50 100 200 400 800 1600 £ (n-n) GeV 
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Fig. 7. - n, as function of In y, in photographic emulsion. 


figure represent the expected <n. as function of <In cot 0) for three types 
of reactions as calculated by the excited nucleon’s model. The curves cor- 
respond to nucleon-nucleon collisions i.e. N =1, nucleon-CNO reactions with 
N>1 and nucleon-AgBr reactions with N > 1 respectively. In evaluating these 
curves the incident energies H,’s weighed according to the primary cosmic 
ray spectrum C-E7?? (25), and a oy = 21 mb is used for the nucleon-nucleon 
interaction cross-section. The contributions of large size cascades are mainly 
at low y.’s, thus when the incident energy is cut off at Z,=1600 GeV the 
averages n, of the AgBr group is underestimated above In Ve = 2.D. To over- 
come this difficulty, the incident energies are extended up to 3200 GeV. To- 
gether with Iny, the effective nucleon-nucleon incident energy E(n—-n)= 
= m,(2y:—1) is given also in GeV. 2 

The individual events of the three groups observed in photographie emul- 
sion are plotted on Fig. 7, and are seen to fall around their respective averaged 
curves as calculated by the excited nucleon’s model. The average values of 


(26) V. L. GINZBURG: Progress in Elementary Particles and Cosmic Ray Physics 
vol. 4 (1958), p. 343. “ies 
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the observed multiplicities of stars with N,<2:; 2 Ni aN ei SN onde 
events are compared in Tables IT and VI with the predictions of the excited 
nucleon’s model. 


TABLE VI. — Comparison between n, observed in nuclear emulsion and the results 
of the Monte-Carlo calculations (M.C.). 

M.C. ) | | MC. [ota EC TRE 

Observed; tie Observed: | Observed : M 
FEE CNO T1= N i | AgBr | esa 7 average! 
DT N>1 n | V>1 iano collisioni 

aa | e n ns In y, Ng INT Iny,| n n n, 

RE CA SESSO DI Eri = | PR = 

1202 SES 8 ZOO MSOs 8215 8 28.2 1.94 | 17 15 23.1 
240A I TOA MISI 23:60 02:381] 29.9 5 | 33.0 | 12.4 20.8 “110 25.6 | 
2.69 | 15.6 | 3 | 34.5 | | 
2.800 "28.4 6 | 28.2 | 2.90 |26.2 | 16 2820 
3.45 | 27.0 Il Diep 9 191987 4 | 38.0 3.48 | 10.330) 31:4 | 

(*) The »s at this energy is underestimated, because 11 of the events are with N7< 2. 


It appears from these tables that the multiplicities predicted by the excited 
nucleon’s model in general are in reasonable agreement with the observed 
values. The agreement between theory and experiment is fairly good at high 
and intermediate energies, but at lower energies the model predicts too high 
values. This defect can be corrected by using a smaller value of in the 
excitation function of the excited nucleons. The effect of this change will be 
to reduce the inelasticity of the nuclear collisions. The same effect can be 
achieved by allowing semi-elastic collisions in which only one nucleon is 
excited. 

The excited nucleon’s model predicts a ratio of 1:10 of nucleons to me- 
sons in an average collision in nuclear emulsion around 1000GeV. If we 
allow 10% (or 20%) K-particles among the produced mesons we shall find a. 
value of (1+1)/11= 0.18 (or (1+2)/11 = 0.27) for the fraction of non-pions 
among all produced particles, in good agreement with the value of 0.27 + 0.8 
observed by EDWARDS et al. (7°). 


9. — Transverse momentum distribution of mesons and nucleons. 


The transverse momentum (P,) distributions of mesons produced in nu- 
cleon-nucleon and nucleon-nucleus collisions were calculated by the excited 
nucleon’s model. Few samples are shown in Fig. 8. All these distributions, 
for different types of nuclear collisions and over a wide range of incident ener- 


me 
~ 
n 
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oies are quite similar: most of the mesons have transverse momentum below 
1 GeV/c with a peak around 0.350 GeV/c. However at higher incident energies, 
above 800 GeV, the P,-distribution becomes flatter, with a lower peak and a 


0.2 
| Nucleons Pions Pions 
04 | E=200 GeVinAg Br £,=800 GeVin Ag Br £,=100 GeVin Ag Br 
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0.2 | : 
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O2 


Nucleons Pions Pions 
E=200 GeVin n-n F,=800 GeVin nn £,=100 GeVin n-n 


0.1 


1 2 3, O5 1.0 OS 10 
G eV/ c 


Fig. 8. — Transverse momentum distribution of nucleons and pions. 


longer tail toward higher momenta. We believe that the P,-distribution can 
be corrected by assuming that the excited nucleons are produced with an angu- 
lar distribution collimated along the primary direction instead of the isotropie 
distribution used so far. 

The experimental observations by EDWARDS et al. (1) and others (1:14) on 
the P, distribution of mesons produced in nuclear interaction, agree well with 
the results of the Monte-Carlo calculations. A comparison between exper- 
imental and theoretical results is given in Table VI. 

A few samples of the P, distributions of nucleons in nucleon-nucleon and 
nucleon-nucleus collisions are shown in Fig. 8. It appears that while the nu- 
cleons produced in nucleon-nucleus collisions have a peak in their P,-distri- 
butions, the P,-distributions in nucleon-nucleon collisions are flat. However 
in all these cases most of the nucleons have transverse momentum between 1 


and 2 GeV/c (compare Fig. 8) in good agreement with the experimental obser- 
vations of EDWARDS et al.(18). 
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spectra obtained for the CNO target. The integral energy spectra of nucleons 
and pions produced in .N-CNO collisions are given in Fig. 10. 


11. — Conclusions. 


The excited nucleon’s model offers a consistent theory for meson pro- 
duction in nucleon-nucleon and nucleon-nucleus collisions between 100 and 
1600 GeV. It predicts with reasonable agreement the charge multiplicity and 
P, distribution of secondaries in nuclear events induced by cosmic ray par- 
ticles. Comparison between theory and experiment shows that the excited 
nucleon’s model can be improved in the following ways: 1) Meson production, 
in nucleon-nucleus collisions mainly at low energies, could be reduced. This 
can be achieved either by using a lower value for x in the excitation function 
or, by allowing collisions in which only one excited nucleon is formed. 2) The 
P, distribution of mesons above 800 GeV has to be more concentrated around 
its peak. This can be obtained by assuming that the excited nucleons in the 
O.M.S. of production, are emitted anisotropically, and with a strong collimation 
along the primary direction. 


TABLE VII. — Characteristic of P, distribution of mesons. 
Incident energies | Peak momentum | È i 
(GeV) | (GeV/c) Observed or calculated by 
100-- 500 0.800 —0.400 EDWARDS ef al. (13) 
500 1 000 | 0.300 0.400 EDWARDS et al. (18) 
30 000 0.500 EDWARDS et al. (1° 
| 38 000 | 0.365 4-0.030 MINAKAWA et al. (12) 
84 000 | © 0.407-40.032 MINAKAWA et al. (12) | 
100 | 0.320 This model with H,—400 MeV 


The study of the excited nucleon’s model yields some results applicable 
to other models as well. The Castagnoli relation In, =n cot) will be 
valid for all models of meson production in which particles are emitted from 
one or more centers fluctuating around a system moving with a Lorentz 
factor y.. Therefore the experimental verification of the Castagnoli relation 
does not necessarily prove the emission of mesons from a single center. In a 
similar way the P, distribution of mesons and nucleons at this range of energy 
seems to be rather a kinematical property, connected with the production of 
low energy mesons from fast moving centers. On the other hand, it appears 
that the real test for the validity of a theory of meson production is the 
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experimental verification of its dynamical properties such as: 1) meson mul- 
tiplicity as function of energy in nucleon-nucleon and nucleon-nucleus col- 
lisions; and 2) energy spectra of mesons and nucleons in nuclear reactions. 


We would like to thank Prof. C. L. PEKERIS for making the wEIZAC elec- 
tronic computer available for this work; the crew of the computer for their 
kind collaboration; and to Mr. R. YA’arr who participated in the early stages 
of this work. 


RATS S UN Ons) 


Si descrive un modello di produzione di mesoni da nucleoni eccitati per collisioni 
nucleone-nucleone e nucleone-nucleo prodotte da particelle di raggi cosmici. I risultati 
vengono confrontati con le osservazioni nelle emulsioni fotografiche, e si trova un 
accordo ragionevole. 


(*) Traduzione a cura della Redazione. 


19 


N 


IL NUOVO CIMENTO Vos MOINE, Nc dl 1° Luglio 1960 


Theory of Low Energy Nucleon-Nucleon Scattering - I. 


D. AMATI, E. LEADER (*) and B. VITALE (*) 
CERN - Geneva 


(ricevuto il 21 Marzo 1960) 


Summary. — The nucleon-nucleon interaction at energies below the 
threshold for inelastic processes, is studied using a method based on 
the Cini-Fubini approach to the Double Integral Representation tech- 
nique. The essential idea is to treat all the singularities of the ampli- 
tudes, for values of the variables lying near their physical region, taking 
full advantage of the symmetry of the Mandelstam representation. The 
spectral functions are calculated using unitarity, both in the nucleon- 
nucleon and in the nucleon-antinucleon channels. In the latter the two 
pion contribution is calculated in terms of pion-nucleon scattering and 
contains amplitudes characterising the low waves of the NN +27 process. 
The construction of integral equations for the partial wave amplitudes 
is discussed. The nucleon-nucleon amplitude is developed in terms of 
linear combinations of the Fermi operators chosen to have simple pro- 
perties under crossing, and whose invariant coefficients are shown to 
satisfy a Mandelstam representation in 4th order perturbation theory. 


1. — Introduction. 


We shall develop in this paper a theoretical study of the nucleon-nucleon 
interaction at low energies; i.e. valid for energies below the threshold of pion 
production, and perhaps for energies at which the inelastic scattering is still 
negligible. 

The method we shall use is based on the two dimensional spectral repre- 
sentation proposed by MANDELSTAM (!) and follows a more simplified approach 


(*) Elsie Ballot Fellow on leave of absence from St. John’s College, Cambridge. 
(**) On leave of absence from Istituto di Fisica Teorica, Università di Napoli. 
(1) S. MANDELSTAM: Phys. Rev., 112, 1344 (1958) 
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introduced recently by Cini and FUBINI (2). This approach consists mainly 
in taking correctly into account only the contributions to the causal amplitudes 
due to singularities near the physical region of the variables: the contribution 
of higher singularities being expressed by some frée parameters in the final 
results. In the Cini-Fubini paper the case of the scattering of scalar neutral 
nucleons is discussed; the physical case of interaction of two fermions through 
a charged pseudoscalar meson field is dealt with in the present paper by a 
suitable extension of the Cini-Fubini method. 

Application of dispersion relation techniques to .V-.N scattering was first 
developed by GOLDBERGER, NAMBU and OEHME (*). They made use of dis- 
persion relations, for the coefficients of a particular set of invariant operators 
at fixed momentum transfer. The two pion contribution to such dispersion 
relations was recently worked out by GOLDBERGER and OEHME (‘) in a way 
such that the dependence on one of the momentum transfer variables has 
been explicitly exhibited. However, because their starting point is a one 
dimensional dispersion relation, the dependence on the other momentum transfer 
variable (which dependence contributes to the angular variation of the scat- 
tering amplitude) is hidden in an arbitrary unknown function. Alternatively 
such a method only exhibits the properties of exchange forces (in lowest order 
due to the exchange of a charged pion between neutron and proton) and not 
the properties of the direct ones (neutral pion exchange in lowest order). We 
shall discuss more extensively the differences between our approach and that 
of GOLDBERGER, OEHEME and NAMBU at the end of this paper. 

The two dimensional dispersion relation approach to nucleon-nucleon scat- 
tering has been applied by CHEW (?) to the determination of the pion-nucleon 
coupling constant g from the angular distribution in .N°-.N scattering; by CENT, 
FUBINI and STANGHELLINI (5) in order to obtain fixed angle dispersion relations, 
which are used both for the determination of g and for an analysis of the low 
energy deviation from the shape independent approximation in proton-proton 
scattering; the latter problem has also been treated by Noyes and WONG (7) 

We shall proceed as follows: in Section 2 we shall define the kinematics 
and separate the causal amplitude into a suitable set of invariant operators 


(2) M. Cini and S. FuBINI: Theory of Low Energy Scattering in Field Theory, in 
Ann. Phys. (to be published). 2 

(3) M. L. GoLDpBERGER, Y. NamBu and R. OEHME: Ann. Phys., 2, 226 (1957), 
henceforth referred to as G.N.O. 

(4) M.L. GoLpBERGER and R. Onume: Application of Dispersion Relations to Nucleon- 
Nucleon Scattering: the Two Pion Contribution, EFINS 59-64 (to be published), hence- 
forth referred to as G.O. 

(5) G. F. CHEW: Phys. Rev., 112, 1380 (1958). J 

(6) M. Cini, S. Fuprni and A. STANGHELLINI: Phys. Rev., 114, 1633 (1959). 

(7) H. P. Noyes and D. Y. Wone: Phys. Rev. Lett., 3, 191 (1959). 
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(which we shall choose as linear combination of Fermi invariants, symmetric 
or antisymmetric under crossing). In Section 3 we shall write down the two 
dimensional representation for the scalar functions which multiply those in- 
variant operators, and express them later on by means of dispersion relations 
at fixed energy. This will allow us to use unitarity in the nucleon- antinucleon 
channel and determine the corresponding two pion contributions. The way 
of dealing phenomenologically with the intermediate states with energy higher 
than two pion masses is then analysed. In Section 4 the way leading to 
integral equations for the individual angular momentum amplitudes is discus- 
sed, and the results obtained are briefly summarized and discussed. Appendix it 
contains a general analysis of the problem of the choice of a set of invariant 
operators and Appendix II contains some technical details on the separation 
of the s- and p-wave two pion intermediate states from the states with 
higher angular momentum. 


2. — The scattering matrix. 


21. Kinematics. — We wish to discuss the elastic scattering of two nucleons, 

with initial and final 4-momenta n,, p, and n,, p, respectively (see Fig. 1). 

Our process will be characterized by three inde- 

TN-N Channel pendent 4-vectors (because energy-momentum 

ni Np conservation gives: 2,+ Pi = M: + Pr), which we 
shall choose as: 


> N-N Channel N.= (ns + n); 
(2.1) P= (PPS 
p p 
2 A = (m — M2) = (Pr — Pr) 


Ich Fea 
Wy LG Le AE, 
As the initial and final particles are on the mass shell, only two independent 
cn can be constructed from these three vectors; we shall define the three 
scalars: 


| WET 


9 9 
(2.2) ù ni) 
| Ne 
where e.g. pì = p?— pì = 2 
ere €.g. D; = P’— Pp, =— m*, where m is the nucleon mass. 


These are connected through the relation: 


(2.3) w+t+t=4m?. 
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We shall also need in the following a set of four orthogonal 4-vectors, which 
we shall choose as: 


(2.4) OV Nee Laer N PAL 


à = Equve M,K,4 


(0) 


Definitions (2.1), (2.2), (2.4) give: 
(2.0) N— Pa To i); AA=—t; Mis K3=--}: Ir--wii. 
In the C.M. system of the two nucleons (2.2) reduces to: 


w= A(k?+m?) = 4E? 


(2.6) t = — 2k*(1— così) 


t =— 2k*(1+ così), 


where %, E, and 0 are, respectively, the momentum and the energy of each 
nucleon and the angle of scattering in the C.M. system. — 

The same 4-vectors n,, #,, p, and p, will describe also the process of elastic 
nucleon-antinucleon scattering; for instance, that process in which ,, p, and 
— N», — Pi are the initial and final momenta of the nucleon and antinucleon 
respectively. In the nucleon-antinucleon C.M. system (2.2) now gives: 


w= — 2p?(1 + cos y) 
(2.7) t= 4(p?+ m?) = 4F, 
t = — 2p*(1— cosy), 


where p, E, and y are, respectively, the momentum and the energy of each 
nucleon or antinucleon and the angle of scattering between the two nucleons 


in the C.M. system. 


22. Scattering and causal matrices. — We define a transition amplitude D 
through its connection with the S-matrix: 
(2.8) Sox = Opa + 1045 — Pa) D pa « 


Tn the case of elastic nucleon-nucleon and nucleon-antinucleon scattering the 
transition amplitude is connected with the differential cross-section in the C.M. 


system by the relation: 


(2.9) 0.10, pi |D;,hP, 


where £ is the C.M. energy of each particle. 
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bo 


The scattering amplitude T is such that 


(2.10) 03,(9, ®) | TO |? 
so that 


(2.11) D, = (1/aB)T;,. 
The covariant Feynman amplitude F is given by 


(2.12) F,, = 4n(Ejm*)T,, - 

The causal amplitude M, with which the rest of this paper is concerned, can be 
defined as in G.N.O. (form. (2.14)) and coincides with the Feynman amplitude 
for physical scattering processes. 

We shall deal in the following also with the nucleon-antinucleon absorption 
into two pions (NN — 7x); for this process it will be convenient to introduce 
an absorption amplitude t (5), which is connected with the scattering ampli- 
tude by: 

1 m,/q 


os T, (NN = rx) = TE 
( ) Bol N° > TT) ox FE /3 Tha » 


where p and q are the C.M. momentum of the nucleons and pions respectively. 
(We note here that the relation (2.10) still holds, when o(, y) describes the 
differential cross-section for the process NN 77.) We have therefore, from 
(2.10), (2.13): 

ql M Ton 


(2.14) Oa (WN? > HT) = We 
pal a p | 2H zr | 


The absorption amplitude is directly related to the Feynman amplitude 
for the process (NN 27): 


(2.15) Tg. = — Zu Fe {NN TT), 


where w is the 7-meson mass. 


23. Separation of the causal amplitude M into invariants. — The causal 
amplitude M is a covariant matrix in the spin and isotopic spin space of the 
incoming and outgoing particles. When using matrix elements M,, for the 
elastic nucleon-nucleon and nucleon-antinucleon Scattering, we shall assume 


(5) G. F. CHEW, M. L. GoLDBERGER, F. E. Low and Y. NAMBU: Phys. Rev., 106 
1377 (1957). | 
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in the following that the spin and isotopic spin operators in M are saturated 
as follows: 


My (NN > NN) = (n) Wel po) MY w, (ns) up), 


(2.16) UM, (NN > NN) = 0; (— mn) vp.) Mi un) x(— p,) = 


Ba 


= B,(n") Ux(De) Mii u;(m)0,(0') , 


where M*’, by the general substitution law (?), is the same function for both 
processes. Note that the spinor saturation of M is between », and », and 
between y», and p, independent of whether they are nucleons or antinucleons. 
Since the operator M will be written as a sum of products of spin operators, 
we shall simply indicate the saturation of these operators by an index (n) or (p). 
The quantities (9) and (9) are positive energy spinors satisfying respectively 
the equations: 


(iv:q+m)u(g)= 0, 


(iva —m)vg)=0, 


4 


where g, =-+-(q?+ m?) 

Besides the spin and isotopic spin operators, M, being a covariant matrix, 
will depend only on the scalars that can be formed by means of the three in- 
dependent vectors (2.1); in particular, we shall take M as dependent on the 
scalars w, t and ¢ defined in (2.2). Therefore M can in principle be written 
as a linear combination of a set of invariant operators, 0;, whose coefficients 0; 
will be scalar functions of w, t and È: 


(2.18) My, > o,(wti) Bl O;|& 


(or course 0; will contain, in the general case, not only the spin and isotopic 
spin operators but also the dynamical variables, as, for instance, in GINO} 
There is a large freedom in the choice of the O;; each problem requires 
a choice of invariant operators most suitable to it. We shall see that our cal- 
culations, where use of the two dimensional dispersion relations is made, will 
require in an essential way a choice of invariant operators such that their 
coefficients can be expressed by a Mandelstam representation. Further a 
simplification occurs if these scalar functions behave simply UL crossing 
(Gt), which corresponds in the nucleon-nucleon channel to the interchange 
of the two initial (or final) particles. 
(9) J. M. Jaucx and F. ROHRLICH: The Theory of Photons and Electrons (Cambridge, 
Mass. 1955), pp. 161, 258. 
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Perturbation calculations up to the 4-th order indicate, as will be shown 
later, that the most useful set of invariants is given by what we shall call in 
the following «perturbative invariants ». 

This set of 0, is given by (we disregard here the isotopic spin part of 0,): 


| Pra; P,=(iyt-pl+1"iy? N); 
(2.19) 
| Rici DER Ri 


These invariants, however, present no simple behaviour under crossing 
(which is equivalent to a Fierz transformation on them); we have therefore to 
find a new set of invariants, linear combinations of the P defined by (2.19), 
with suitable coefficients that do not spoil the analytic behaviour of their cor- 
responding scalar functions, and which transform in the simplest way under 
a Fierz transformation. This had led us to the use of linear combinations of 
Fermi invariants. 


We can express the perturbative invariants P,; in terms of the Fermi invariants F, 
defined as usual by: 
| ESS = len FV a= yh ye 
(2.20) | 


FINE TES UNION PTE NT PVI APE CA I NEED ED AE ET al aoe 
Py = 1 => vale): F,= A a (05 Wa) (U5 Yn) F;= EP — VE VS > 
US v 


[see Appendix I] by the relations: 


PS? 
1 | —w 7 = 

RAZZE lg mV 4 pei, i 
m 4 4 

(2.21) 

pid (one 4m?P |, 

ae 

IE = JP 


and we note that rio additional singul 


wa arities are introduced in going from the P, to 
he F;. : 


2 4. ESET symmetry. — It is easy to see that five independent invariant 
operators which are symmetric or antisymmetric under 


crossing (i.e. under a 
Fierz transformation P(2, 4) 


which brings from [u(n.)O™ a(n n D up.) 
to [u(n) 0Mu(p)] [a(p,)O%u(n;)] and : ee Wu) | ECs u(ps) | 
ZIO Di P:)OVu(n,)|] and therefore exchanges # into #) can be con- 


structed as linear combinations of the Fermi invariants combinations which 
Ds ) D 
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we shall denote henceforth by 0;. The choice of C; will eventually depend 
on, arguments of simplicity. 
Let us define the index j in such a way that: 
0; is antisymmetric under P(2, 4) if j is odd (j MESI 
C; is symmetric under P(2, 4) if j is even (j — 2, 4) 
so that: 
99 dai 
(2.23) C; Peay (— iL) C; . 
We shall now go back to (2.18) and separate the spin and isotopie spin 
parts in 0,. Let us introduce the two projection operators in the isotopic 


spin space, P, (T— 0,1 is the total isotopic spin of the NN system), such 
that: 


(2.24) PINNA q(T) 

(in terms of the isotopic spin matrices t” and +’, P, and P, are given by: 
(2.25) Py = 4 — t** t?) ; P, = (28 + 1-7), 

where, as usual, t” is to be saturated between (n,) and w(n,) and 7? between 


U(p.) and u(p,)). 
We can therefore write (2.18) as: 


As a P(2,4) Fierz transformation is equivalent to an exchange mn =, 
the Pauli principle applied to the initial NN state will give: 


(2.27) ‘na pa |M [np = — (Map | MP), 
where position, spin and isotopic spin of the two incoming particles are inter- 


changed. 
By using the crossing symmetry (2.23) and by noticing that: 


(2.28) (ny ps | P,| mp) = — (1) ps | PP) , 


A "6 a om ; = D2 EA 73 AP ni à = 
we obtain, from (2.23), (2.26), (2.27) and (2.28), the following crossing con 


ditions on the cf: 


9,99) ci (wit) = (— 1) +7 6? (wtt) = ¢,,¢;(wit) , 
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which defines the meaning of the symbol &,,, which will be used in the fol- 
lowing to simplify our notations. 

We shall now briefly discuss the choice of the set of C;s which we shall 
use in what follows. As our method will eventually lead to integral equations 
for the transition amplitudes corresponding to definite initial and final spin 
states (helicity states or singlet and triplet states), a simplicity criterion can 
be provided by the relations between the coefficients of the O,’s, ¢;, and the 
above-mentioned transition amplitudes. A set which turns out to be conve- 
nient from this point of view is given by: 


SESTA eA 
MOST AE OP ye 
(2.30) ee (OS = Vas aie 
= 3847 + 3P\2 
DEA SIETE) I 


(we recall that the O;’s with odd index j are odd under a Fierz transformation 
P(2,4) and those with even j are even). 


The usefulness of the set (2.30) can be easily seen by expressing the C; as Fermi 
operators, defined by a saturation between spinors different from the standard satu- 
ration (2.16). We shall here merely sketch the argument which leads to the connection 
of the O,’s, as given by (2.30), with the transition amplitudes for definite helicity states. 
Let us symbolize for simplicity the standard saturation (2.16) of the Fermi invariants 
between spinors by the expression (1, 2; 3, 4). If now we define a new set of Fermi 
invariants, 7;, by sandwiching the Fermi operators between initial-initial (and final- 
final) particles, we shall get, for instance, an expression of the type: (2, 4; 3, 1). It is 
easy to see that this permutation of the spinors can be obtained by the consecutive 
operation of a P(1, 2) and a P(2, 4) (Fierz) transformation to the standard order (1, 2; 3, 4). 
P(1, 2) is diagonal [see ref. (14)] and gives plus sign for S, A and P and minus sign 
for V and 7; if we now apply the product matrix P(1, 2), P(2, 4) to the Fermi inva- 
riants we get the linear combinations given by (2.30). It is therefore demonstrated 
that those combinations reduce each to a single Fermi invariant, saturated between 
the two initial and final particles. With the saturation (2, 4; 3, 1) one finds that: 


(2.31) S=0;; V=0,; T=0;  A=6,, MP 


| To go now from the C's to the scattering amplitudes corresponding to definite 
spin states, we have to saturate the invariant operators with spinors with definite 
helicity or with combinations of them corresponding to singlet and triplet states. It 
is clear that this connection will be much simplified, in the nucleon-nucleon C.M. system 
by our new saturation; our invariant operators will be represented by single Vera 
invariants, even or odd under the exchange tf, and the y o contained in 
them will now connect the two incoming (outgoing) particles, which have opposite 
momenta in the C.M. 


84 


THEORY OF LOW ENERGY NUCLEON-NUCLEON SCATTERING - I Wei 


As an example, we shall give here the connection between the singlet and triplet 
transition amplitudes 7; (as defined, for instance, in G.N.O.) and the coefficients C; 


of the invariant operators O;: 


mm | 
Ti = ——[nl2n — (y —1) sin? 0], — 2(n? — 1) cos Oey + [2 + (1 — 1) sin? 0]e, |; 
47) 


m 
Tim = = (y —1) sin? One, — ag), 
oa) 
2m 
Ig = ven (7 — 1) sin 6[n cos 6c, — (n + l)e; — cos 0G], 
01 ] ] 2 Ù 
470 
(2.32) 
2m 
= v — (n —1) sin 0[(n + 1)6, + cos Oc, — 7 cos OG], 
4a 
m 
To 3 [ (7? — 1) cos 0e, + [7 — (n —1) cos? 0] + n[1 + (7 — 1) cos? OCAE 
2711] 
m 
tel) code, yi], 
277 


where = 4/w/2m and the six 7;; are reduced to five independent amplitudes by the 
relation : 


(2.33) (Ty, — Tia — Too) = V2 ctg (Ti+ To) » 


which follows from the time reversal invariance of the theory. 


By using (2.21) and (2.30) it is easy to find out the relations between the 
invariants P; and the Cÿs: 


P,= 310 + C++ C2 + Ci + 05], 


L à! | 
P, = LL (m2 —w — 24) 6, — 2m 0, + 2m 0, + — (2m* —w) 0. + 
i 4|m 
LL 14 (2m2 + w+ 21) Ci |, 
m 
2. = =: ; 
exe ley = L [(5m? — 2w — 24) C, — (3m? w) 0, + (m2 — w) 034 
4 
+ m?0,— (3m? — 2w — 21) Cs], 
1 ; 5 
JB: = sl 20 | 0; Ca + 205], 
5 Y Y C | C I Y | 
ie ree 0, Si] ARS 
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25. Isotopic spin in the nucleon-antinucleon channel. — We shall find useful 
in the following to label our matrix elements according to the total isotopic 
spin I in the NN channel, while our label T corresponds to the total isotopic 
spin in the NN channel. The connection between these two sets of T- and 
I-labelled amplitudes can be easily obtained in the following way. 

Let us divide M,,(NÒN >), which is a matrix in the isotopic spin space, 
into two parts: 


(2.35) MU (NN mr) = Nd, + ET, GIN, 
where N, and N_ are proportional to the transition amplitudes with definite 
total isotopic spin J= 0, 1: 


(2.36) No) (WAG INS Ne =a Ne 


From (2.35) and the orthogonality of N. and N_, which follows from (2.36), 
we obtain: 


(2.37) U(NN > NN) =3NIN, +2N*N vor = 3 Mt + 2M x2? 


As M = P,M°+ P,M (where the label 0, 1 now refers to the total isotopic 
spin T in the NN channel), using (2.25) it is easy to obtain the desired re- 
lation: 


| me=3ut—oM-, 
(2.38) 


| M=3Mt +2: 


3. - Two dimensional dispersion relations. 


31. — Let us consider the fourth order Feynman graphs of Fig. 2. 


See ears KEANE RE 
i i xX de È I 1 x oi 
I ! x ' I ea 
' il sa > ‘ : BAN 
È 7 a 3 AT è CAR 
I 1 I 2 P, Il Na P, = n, 


Fig. 2 


am nee computat'on of graph I (the results of the perturbative calcus 
lation will turn out to be a particular case of the results of this section) shows 
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that the coefficients of its expansion in terms of the perturbative invariants 
(2.19) have the simple form 


DR ns 0 
el fui @=w)y 0 


(2m)? (2)? 


The graph IT gives just the same contribution to those coefficients (except 
for a plus-minus sign (*)) with w substituted by £. This means, as we have 
discussed before, that the contributions of graphs I and II to the coefficients 
of the invariants C; of (2.30) have the necessary analytic properties in order 
to satisfy a Mandelstam representation (except for one more subtraction due 
to the pole at infinity introduced by the transformation (2.34)). 

Graph ITI can be obtained from graph I by exchanging », into p, and vice- 
versa. This means that the coefficients of the « perturbative » invariants with 
the y matrices saturated between n, and p,, and p, and », (crossed perturbative 
invariants) will have the simple form: 


o(xy) 
(3.3) | di | dy i 


(2m)? (2)? 


Now, what can be said about these crossed perturbative invariants? They 
will be related to the perturbative ones of (2.19) in a very complicated way 
in which extra singularities will appear. So we shall not be able to assert that 
the four graphs of Fig. 2 give coefficients of perturbative invariants satisfying a 
Mandelstam representation. However, the crossed perturbative invariants 
will be related to the crossed 0, invariants just in the same way as the per- 
turbative ones where to the 0,. But the crossed C, are just the C; except 
for a + sign: this means that the contributions of eraphs III and IV to 


(*) This symmetry can be understood easily from (2.21). In fact, the amplitude 
for graph II can be obtained from that of I by the substitution n,<—n,. This sub- 
stitution makes w + 7, changes the sign of S, A and P of (2.20), and leaves V and OF 
unchanged. From (2.21) one sees that P,, P; and P; are antisymmetric for this ex- 
change while P, and P, are symmetric. Besides, the isotopic state (+) introduced in 
the last section is antisymmetric for that exchange while the (—) state is symmetric. 
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the coefficients of the 0, have again the correct singularities so to satisfy a 
Mandelstam representation. 

So generalizing the results of the fourth order perturbative calculation we 
shall suppose the validity of the Mandelstam representation for the functions 
ci(wtt) of (2.26). We note that up to now there is no more indication for the 
validity of Mandelstam representation, even for the case of scalar particles, 
than that given by perturbative arguments. 

-Tnstead of writing the representation in the usual way, we shall put in 
evidence the kinematical coefficients coming from (2.34) by writing: 


SA 6m? —w—2t 5, i 
i(w't) =| pi(wtt) + —— pi(wit) + (5m? — 2w — 2t)- 
i Mm 
ue a a ite i 
- pa(wit) — 4pi(wti) + pico) ( 1)7{t< >t} LL = 1( —, 
rs 


rl 


A al 
co(wit) = fi [pi — 2mpî — (3m? — w) ps + 2pî ni 


t — t) 
oh 
bada 105 100 
= 1 
ci(wit) = n [pt + 2mpî + (m2 — w)p3 — 2piî TA} 
(3.4) Nes 
SI Latina) 
Wy — W 
È il 2m2 — w 
Te È fps + SPE" pi ment + vi] + 
m | 
Dit 
LED 
Eat HR ne 


2m? + w + 21 
(wt) =; Di pela Ra. - po (3m? — 2u 2t) ps + 4pi Ds | 
DI n 
pe Lt >t} — dec di 

Wow 

where 

ee | pt =3pt+ 2p- 

| PES pt Oe 
with 
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D,(i—t) is the residue at the deuteron pole at w, = (2m —B), B being the 
magnitude of the deuteron binding energy. D, is given in terms of B and the 
asymptotic ratio of S and D waves in the deuteron wave function: explicit 
expression for such residues, as coefficients of two dimensional invariants, are 
given in formula (4.43) of GNO. g is the usual renormalized rationalized pion 
nucleon coupling constant, i.e. g?/4a — 14.4. 


Clearly (3.4) to (3.6) are absolutely analogous to a subtracted Mandelstam repre- 
sentation for the c;(wtt) written in the usual form (12). The cut in w begins at (2m)? 
while that of ¢ or £ starts at (2)? which corresponds to the lowest energy of the process 
NN — 27. Of course the range from (24)? to (2m)? corresponds to an unphysical region 
of energy for the NN system. Another feature is automatically taken into account 
in (3.4)-(3.6); this is the fact, true at all orders in perturbation theory, that ¢ and f 
cannot both be in the region (2)? to (2m)?. By simple inspection of a general Feynman 
diagram one immediately realizes that if looking in one of the NN? channels it does 
not have any baryonic line in an intermediate state, it presents at least two baryonic 
lines in the intermediate states of the other NN channel. 


The lowest cut in t in (3.6), de. that going from 4u? to 9u?, is due to the 
intermediate two pion contribution; let us single out this contribution by 
writing (3.6) in the form: 


=(apt! “v= (at i x 

1. HE CAT NÉ i oe 1) [ro > é], 
dd To} C—O 

au? am? 


where o(t’w) has a cut in w starting at (2m) and y(xt) one in ¢ that, however, 
begins only at 9u?. Following our general approach we Shall make the as- 
sumption that in N-N scattering (in which ¢ and ¢ are negative) the depen- 
dence of the last integral in (3.7) on t is sufficiently weak so as to contribute 
only to the lower partial waves. We shall discuss in Section 4 the consequences 


of such a hypothesis. 


32. Evaluation of the two pion contribution. — We shall use now unitarity 
in the nucleon-antinucleon channel in which ¢ is the energy and w and? are 
the two momentum transfers (see (2.7), where the relation bteween vw, è and ¢ 
and the momentum and scattering angle in the NN C.M. system is given). 

The unitarity condition in terms of the transition matrix D defined 


by (2.8) is: 
(3.8) (Di, —D,,) = > O*(P, — Pa) Dp Pye . 
y 
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co 
bo 


Relation (3.8) leads, when written in terms of the causal matrix M and of 
the absorption matrix + for the process NN — 27 (see (2.13), (2.14), (2.15)), 
to the following relation : 


> NN |t+|an> <an|t| NN) + J,, 


(3.9) i<NN|Mt— M|NN) = Got 8 
where the sum is carried out over all 2x intermediate states that conserve 
four momenta. The contributions of the intermediate states in M*M have 
been split into those coming from intermediate two pion states and those co- 
ming from higher energy intermediate states with baryonic number zero (J,). 

Also, one sees easily from (3.7) and the fact that all the invariants we have 
introduced are self-adjoint, that the two pion contribution to 


(3.10) i(M+— M) =2 > (Ime,)0; 


is given by means of (3.4), (3.5) in terms of 


(3.11) Im (p;) = 0;(0t) = (—1)0,(#) . 


Besides, because we have put in evidence in (3.4) the kinematical factors coming 
from the transformation from P; to 0;, Im p; turns out to be given by the 
coefficients of the perturbative invariants in the expression for (7/2)(M+— M). 
Then, due to (3.9): 


| Ne AM 

D 17 (200) Se (2115 Sp 

(3.12) È. [o;(wt) = (—1)o(#)] P; = Dam)? 
where 

(3.13) M = > Dip. \t* |20><20|t|m,7%) , 


(27 


so that by evaluating .# we shall be able to obtain from (3.12) the function o 


A corresponds to the graph of Fig. 3, in which the intermediate pions are 
the mass shell, where n= n,, n'=— n, are the momenta of the incoming 
nucleon and antinucleon, respectively and 0 —=Der 

n n' p' =— p, those of the outgoing nucleon and antinueleon. 
eS 2x | t| m7» is the matrix element for the process 

ae NN > 27 and is an analytic continuation of that 


cata Y Y g-44q Corresponding to pion nucleon scattering. 
2 We write it in the usual way (3): 


16) Cerra] ra) =U) (48,5 bey, T4, ) + 


ie 
D 


ie Be | "a 
+ ty Q (Br: + Ets, 7,1B, )lu(n), 
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where x, f are the isotopic spin indices of the two pions and 
R q_-q' 
(3.15) ta 


A, and B, will be functions of 


t =—(n+n')2, 
(3.16) Sn) — (n* 9g" — 2g cop), 
S$, = — (n — q’)? (n? + q? +2nq cos y), 
(3.17) (€ + Sn +S, = 2m? + 2u?), 


where y, is the WN — 27 scattering angle; n and q the modulus of thenucleon 
and pion momentum, respectively, in the NN C.M. system. 

We see that the ¢ defined here coincides with the ¢ of (2.2), s, and Ss, depend 
on the scattering angle in the process NN — 27. Similar reasoning can be 
applied to 


(Paps |e |200 


We split the sum over intermediate states in (3.13) into a sum over the iso- 
topic indices of the pions (which is trivially done) and an integration over 
their momenta; we obtain then 


(3.18) A =3Mt+4 209 - OL, 


(FLO) 7 = Java u(p)(— A#°+ iy: QBs") v(p'): 
-O(n') (—A® + ty -QBr) u(n): 044 + qg'— p —p'). 


Note that, in the region we are interested in (94° <1<4y"), g is real but p is pure 
imaginary. The unitarity condition has therefore to be extended to the complex 
p-plane by using S*(p*)S(p)=1 (see, for instance, N. VAN KAMPEN: Phys. Rev., 89, 
1072 (1953)). The arguments of A* and B* in (3.19) contain therefore p* and not p, 
which is equivalent to a substitution sz25, as is clear from (3.16). The results 
obtained by using this extension of the unitarity condition are completely equivalent 
to those that can be obtained by working always in the physical region (t> 4m’), 
where p and q are both real (and so are therefore s and s), and going then by ana- 
lytical continuation in t to the region we are interested in (4u?°<t< 9"). 


In order to perform the integration over the two 6 functions it is simpler to 
so into the NN centre of mass system. In such a system 4 has no spatial 


components, so that 


(3.20) eee oe 
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iti i inor space, Wi i ion: have 
Writing .#+ as an operator in spor space, W ith our usual notations, we 


(3.21) Ms == (t — 4u*)) | (—AË + iy QBE*) (— Art ty QBr) d9, 
where the integration is now extended only over the whole solid angle of Q, 
its modulus being given by 


(3.22) o=- 


In order to proceed in the evaluation of @ we must introduce in (3.21) 
some knowledge about A and B. From their Mandelstam representation we 
know the location of their singularities, so we can write for them the following 


expression: 


1 i 1 ol ert ss 
A*(sst) == fares t) (, ets ATTE = | a ae 
IT s G ; 


— 8 s'- 8] n= 
(m + 10) (24)? 
1 1 DA FAR LR i 
B=(sst)= 9° — + eil ds’ + 
(3.23) LEE ms ms) a SEI Seo Pe SOIR 
Cm +) 


where the upper sign corresponds to isotopic index (+) and the lower to the 
index (—). 
We distinguish for A and B three contributions; first the pole that repre- 
sents the graph of Fig. 4; next the cut in s and s. This incorporates the x-.N° 
rescattering corrections to A and B and is that studied by CHEW, 
Be A GOLDBERGER, Low and NAMBU for z-. scattering (5); it is domi- 
SERE x nated by the 33 resonance. And lastly, the cut in # brings into 
the game the x-r interaction; this term has been considered by 
FRAZER and FULCO (!°) in the analysis of the electromagnetic 
structure of the nucleon, and recently by Bowcock, COTTINGHAM 
and LURTÉ (!') in the analysis of r-.N° scattering, under the hypothesis of p-wave 
m-r resonance (*). In these last investigations is discussed the location of the 
singularities of v, and v, in s and 5 as derived from the Mandelstam repre- 
sentation for A and B. These singularities are cuts in s and 3 that, however 
begin only at s=(m-24)? and 5=(m+2u)2. The whole philosophy of our 


(19) W. Frazer and J. FuLco: Phys. Rev. Lett., 2, 365 (1959). 


(11) J. Bowcock, W. N. CortINGHAM and D. LurIé: Nuovo Cimento, 16, 918 (1960). 


() Note that o,(s't) and o;(s't) also have cuts in t beginning at 1642 (see (11)). 
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approach consists in making use of the weak dependence of every quantity 
in those variables for which the singularities are far from the region of interest. 
Or, alternatively, we lump the dependence of these variables into some con- 
stants which come from the lower terms «of a development of the quantity 
in question. 

The development parameter is just 


(3.24) $ — 5 = 4pq cos Q, 


g being the c.m. scattering angle in the reaction NN +27. Then a weak 
dependence on s—sS means that only low waves of the two pions occur in 
the process NN >2x. To fix better our ideas let us suppose that only s 
and p-waves occur in the development of the last terms in (3.23); it is clear 
however, that higher waves of the two pions shall be present by virtue of the 
strong s and 5 dependence of the other terms in (3.23). 

We stress that the restriction made is justified only for the values of # in 
which we are interested, i.e. a total energy of the two pion system of a few 
pion masses. The inclusion of higher waves coming from the integrals over # 
in (23) would be a simple extension of the following arguments. 

We notice first that not only the last term in (3.23) contributes to the s 
and p-waves; the other terms that we shall call ./+Z+, where 


L TA 1 Re 
@*(sst) = — Joris t) | ee. 5) ds’, 
1 ‘ 8 - A 


§—s Sii 
o 5) (m + pu) 
(3.25) eee i eee fant do 
Opel e = — IS - =, NE = 1s 
(ss ) x Jane 3 ) ec 8 ! RIESI 3 ? 
0 
with 
(2.26) o*(s', t) = 2g?6(s’— m2) + ox (s', t) 


also contribute to them. In Appendix IT it is shown that the 27 s, p-wave 


contributions of .7 and Z are given respectively by (omitting the + indices 
when inessential) 
sets 
a) Ot) — Tone 
(3.27) 
| p =F Pitt) ’ 
where 
ao (t) =: A, 
Sii dl MO 
x0=3( A1 + A5). 
(3.28) 2 \pq P 
| Br(f) = — #5), 
ETAT RENE 
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+1 
D =| f Ba loos 9) A 0089 

= il 

+1 


B,\t) =| 4P, (00s p)d cose. 
-1 


Let us denote by 4 and B,, the whole s and p-wave contribution to A 
and B respectively. Then, because we supposed that the last term in (3.23) 
contributes only to such waves 


I (es 
ASSE ARRET 
JUNE TT 
(3.30) Nu A 
yp s—s 
Bo PAESE 
8+3 | C= 
(2u)* 


so we can write 


| A= at An, 
(3.31) 
| = 4—8+8,,. 


n 


In analogy with (3.27) we can write for A4, and B,, the development 


(E) 
(3.32) A = At) + ,(t) È i — 
Bs» = y(t) 9 
with 
(3.33) AS = A7 = ni Sg 0 n 


We are now ready to go further in the evaluation of (3.21). For every A, 
B,, A, or B, we shall use an expression of the form (3.31) with (3.27) i 
(3.32) were S— 8 will be either s,— 3, or Sh —$,, respectively. In Appen- 
dix II it is shown that the iteration of the terms containing A,, and B,, with 
Duo containing %—a and 4—f is zero. This is physically clear beosli 
interference between terms containing only s and p-waves (as As, and B ) 
po a — a ni dat which by construction do not contain mene ae 
certainly vanish. Since « ¢ : 
(3.32)) a PR wi ni on ace i n È: Re SR ee 
S Statement is that the iteration of « 


ine] 


Ae 
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and B with 7 and Z gives the same contribution as their iteration with « 


and B respectively (*). 
2 


Then we can write for (3 
[ic A+ iy QB) (— Ly + ty: QB,) + 
xla 4 iy QB.» —Bln)} dQQ : 


1) 


Il t 
[AS GU: T vy: Q(B, 


( (AE, 
It is not difficult to develop (3.34) into invariants. We obtain readily 


(3.34) 
Bn) ( 


M = (a + 2me + m’e)1® 1° — (6 + md)[1" iy? NH iy" PB LE] 
+ e[iy?-N iy" -P] =, fly?-y"] 5 


(3.35) 
where 
| 1 t 2 Li " 9% * * 
di i a owe LA An + Afp Asp, — Xp An] dQ, 
1 t aie *# 3 n) 
b= 3 ti-=- LÀ (BE SD, + Re (Bry, As) — Po%n]" 
ò -M -K 
il t z Di yk 
Cs ( Cr n) fer An + Re (Boy, Asp) — Po%nl 
QM QK 
Tm — Ke | Se 
pone { (1 we) fire. PE el 
So) 8MK \ \4 
x 2 - K)? 
ier (2K? + MP) 4 oo | (2M? + K+) — Q*(K?-+M*)|dQ¢ , 
1 t 5 dires > 9 
Br Ba ar By Bite by Bal: 
È mel we)) | i a 
(E (QE (je 2 + QUM' — K)] d00, 
M: ie 
IL t 3 ap* * A . 
he sl di n | [BA + Bry, Bes, — BoBal 
À - M}? - K)? 
3 e PE ee di — ce "| dQQa 3 


interference of s and p-waves with the whole 
the s and p part of ./ and 9, 


(*) This reflects the fact that the 
< and Z is the same as their interference with only 
i.e. « and p. 
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where M and K are the spatial components of M and K, as defined in (2.4), 
in the .NN centre of mass system. In such a system the time component of 
M and K vanish, so that 


(3.37) M?=—w, ki] Se 


Let us consider the integration over . and Z in (3.36). We see that the 
dependence on 2, for every A or Z is explicitly exhibited in s and 3 (s,, Fn 
and s,, $,) in the denominators of (3.25). Then the integration over Qo can 


be performed explicitly. We prefer, however, to express the results in para- 


fr if 1 x 
| (at) (=, RE gee | da , 
; o—t ax—-w 


metric form 


where the dependence on w of x is determined. We shall see later the reason 
for this form. In order to integrate the terms in (3.36) containing As,Bsp % 
and 6, we take for them their expressions (3.27) and (3.32) in which their 
dependence on 29 appears explicitly in s,—8, = 4Q:N in one factor and 
s,— 8, =4Q-P in the other. The integration over 2, can then be performed 
easily. By collecting the results we find: 


i tain 1 
bs > furia) = 1 | 7 ua 
de È æ — w 6 
am? 
[Re (228) — ates (eft — ya) 
(ni A1) — Br] a 
I D 1 / 
eas) el faq 14 1 
T, | ) reve, dæ 
4m? 
Ah Se 1 do 
FF = SE = ta oo 
ife 2) PR A En à da, 
4m gl 
FI 1 Ft 1 
ee =— | x, (tx ae | 
| ) she, da, 
4m? 
il it 1 | 
+ — = | y= (ty e A = ae i ‘ 
Lf ars Seana) 
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with 


be (ie vi] [as ds" oF (8't) o¢(s"t) K(s's"ta) , 


i /t UE =e 7 gle! = g'gl' 
ie = [Ja as (stat | = 2 Aa i i He 
me D — dm? + i x 


pee ‘to sn RUE) EE, 
pi x 


} xa (CE) ) [far ds" 07(8't) 07 (s"t)- 
(3.39) 


i_ 4m? EE ZL? Wel. 
È a—A4m +1)? K' |(x — 4m? + 0) 


Vir 2A 
x, (te) = aa OR (s't)Og(s"t): 


| 20—4m?+t 1. 72 y2 | 
| a2 — 4m? +i? K || ; 


— 
SD 
[uen 
= 
= 
teste) | 
t | 
Î 
~ 
— 
19 


] 


Wi AM ET) TE 


SE VAI Bie a MEN ENT: | . - | 
Sii dedita (osi e 
xp (at) Fi | fas ds" a7 (8't) 03(8° E (2 — 4m? Lt) Kf’ 


Ca 


with 


UA 


(3.40) ssa) = (a [v/s' + +/8"]*) (@ — [Vs —48"1?) — 4m(s'— 8")? — 


— dar { 3'8"  [u2— m?]?— (u? + m?)(8'+ 8") + wa } Tae 


NT 


(3.41) Z(s's"t) =s'+ s"— Qu? + m°) +6, Y(s's") = ses 


where the integrations over s’ and s” begin at m? and go up to a curve (in 
the s’—s" plane) given by: 


t , Lo 
(3.42) s's"4 (s' + “| UE n°) Em ru)? 7 (t — Au?) 4 


a Vis? Ls'(t—2u? 2m?) + (m?— u2)][s"? Ls"(t—2u? Im) (m? wy) = 0. 


This limits, for any fixed x and t, the possible intermediate masses 4/s! and 4/s”7 
of Fig. 5 that can contribute to (3.39). Alternatively, the 
same integration region of (3.38) and (3.39) can be obtained. Vs" 
by limiting the integration over « for every value of s', s” and 


tin the form 
2 Ve 
for fas far Ss 
SE Fig. 5. 
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where x, is just the value given by solving (3.42) for æ. This means that for every pair 
of intermediate masses +/s’ and 45" the integration over x is limited to lie above 
a curve in the x — t plane given by (3.42) (see ref. (12)). 

If we take of =o#= 0, so that only the 6 functions of 5 in (3.26) remain in (3.39) 
the expressions for the x are just real analytic expressions for t > 4u?, x > 4m?. This 
is, on the other hand, just the perturbative fourth order (renormalized) contribution. 
That the fourth order terms have this form was already utilized in the discussion at the 
beginning of Section 3°1. 


We are now able to give explicit expressions for the functions 9; of (3.7). 
In fact, from (3.12), (3.35) and (3.38) we obtain readily: 


1 1 [xt + 2mxt + mx} 
t x — WwW 


| ete (27 
+ ) i |A |? — aa” 
ni age \- (10/6)((t/4) — pe) (az |? — x) I 
| 


SUO 


1 fae mae 
La. : : th Ka 
02 (wt) == = —- : da 


AC Di 0 | 
1 Nd BW \ ery, A) ae 


+ 
# 
a Gi). 
gelo 


(2m)? 


SEZ 207) 


y il UE TT HD NS 0 
f= hf Pate BV MA 
2 (2x2) G— w 12 4 DE i 


~ (2m)? 


This means, by looking at (3.38), (3.39)-(3.41) that we have obtained for 
all o expressions that depend on o, and o, on one side and on functions 
Ay (t), A, (6), 4, (t) on the other (aj, «, and B; are given by (3.28), (3.29) as fune- 
tions of % and 6,). The former are given by -.N scattering; if for instance 
we take into account only the (3, 3) resonance in x-.N° rescattering corrections 
as done by CHEW, GOLDBERGER, Low and NAMBU (8) one finds 


gi (st) = — 207 (st) = = pa is D (i = ! Vile | sin? 033(8) 
e ale tL 2h? E—m k 3 
Sin? 0, (5) 
‘LR 


(3.44) 


| 6; (st) = — 2o; (st) = — 


(7?) S. MANDELSTAM: Phys. Rev., 115, 1741, 1752 (1959) 
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where 
Su: +(s — m? — w?)? — u?m? 
eee oe, gt 


E == = 
24/8 IS 


and 033 is the usual (3, 3) r-N phaseshift. ‘ 

The functions Z;(t), À, (t), 7, (t) come instead from the NN—=>2x process; 
unhappily we have not much information about it. The fact that they repre- 
sent the whole 2x7 s and p-wave contribution to such a reaction would allow, 
however, the use of unitarity in order to specify their properties. For its 
discussion we refer to the work of Bowcock, CorrIngHAM and LURIÉ (!*) in 
which (under the assumption of a sharp resonance in p-Wave pion-pion scat- 
tering), A, (t), and 7; (t) (A; = 0) are obtained in terms of certain constants. 
In a future paper, consequences of the x-7 resonance in .N°-.N° scattering, will 
be investigated. 


4. — Discussion and conclusions. 


We shall briefly sum up and discuss the results obtained in the preceding 
paragraph. We shall use the information that we have get on the structure 
of the ¢,’§ in order to discuss the possibility of obtaining integral equations 
for the partial wave amplitudes in VN scattering. 

A detailed derivation and discussion of these partial wave integral equations 
will be given in a later paper. We shall here shortly note the main features 
of these equations. 

Tt is clear that our ¢,’s, satisfying a Mandelstam representation, present 
Singularities only on the real k? axis. Going to the 7, by means of (2.32) 
does not introduce any extra singularity, with the exception of a kinematical 
singularity (due to the Vw present in 7) which can be, however, singled out in 
the singlet case. Partial wave integral equations will be satisfied by amplitudes 
proportional to the a7, b7 and ©, (as defined, for instance, in G.N.O., for- 
mula (B, 11) in Appendix IT) and presenting the right analytic and convergence 
properties. It is easy to see from the above mentioned formula in G.N.O. 
and from (2.32) that these functions will be obtained by multiplying the am- 


plitudes a, b and 0 by the kinematical factor: V (k?-+m2)/k?. The resulting 
amplitudes, which can be expressed in terms of the projections of the cs over 
partial waves, will present the following singularities in the k? plane: 

a) a cut for positive values of k?, starting from 0; 

b) a cut going from —p?/4 to — co, due to the one meson pole; 

c) a cut going from —p? to —oo, due to the two pion contribution; 


d) farther cuts on the negative real axis, due to many-pion states. 
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The jump of the amplitude at the cut «) 1 is just given by the imaginary 
part of the amplitude considered, for positive values of the energy; the jumps 
at b) and c) can be calculated starting from the formulae given in the previous 
paragraph, and the two pion contribution can be given explicitly in terms of 
pion-nucleon and pion-pion interaction; the higher cuts will contribute with 
a weak dependence on the energy, which will be approximate by a low order 
polynomial in k?. This means that the different amplitudes will contain some 
constants, which will enter as parameters into our theory. The number of them 
will clearly depend on the region of the variables in which we wish to reach 
definite conclusions, i.e. the energy up to which one is interested in. The 
theory is however meaningful only if the number of parameters is rather 
restricted. This philosophy has as background the assumption that the weight 
functions (of the Mandelstam representations) in all the phenomena we are 
interested in have no «strange » behaviour, so that the main features are still 
understandable just from the positions of the singularities. This is a theoretical 
assumption whose consequences are reflected in our results in which the whole 
aspect of low energy .N-N scattering is given in terms of only a few para- 
meters. 


We want to discuss here briefly the difference between our approach and 
that of GOLDBERGER, OEHME and NAMBU (**). They used dispersion relations 
at fixed momentum transfer ¢ of the form 


ai 1 j 
(A) Ee) = — ale ae =i, 
f OP = Ye TT EST 
(2m)* (CIO 


Re Tiss 1 [= F(a, t) im F(z, 1) 
IT 


GOLDBERGER and OEHME have analysed the last integral by considering the 
two pion contributions to it. 

Let us analyse the relation (4.1) in some detail. First we note that it is 
surely most appropriate to use it in the case of forward scattering (f= 0) in 
which, by virtue of the optical theorem, the imaginary part of the amplitudes 
in the first and second integrals in (4.1) is connected to N-N and N-N total 
scattering cross-sections respectively. Im this case (4.1) is transformed in a 
sum rule similar to that analysed by Puppr and STANGHELLINI for 7:-.N° scat- 
bering (18) (*). 

But if we want to obtain more information from the theory (i.e. phase 
shifts) we are forced to use the dispersion relations also for ¢ # 0, for instance, 
by doing a development in partial waves. This means to integrate { over the 


(9) G. Pur and A. STANGHELLINI: Nuovo Cimento, 5, 1305 (1957). 


(*) Added in proof. - An analysis in this direction has been recenly attempted by 
Y. Hara and H. MIYAZAWA (Tokyo U niversity preprint). 
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appropriate Legendre functions. It is clear, however, that in regions in which 
there are singularities in ¢ (poles or cuts) such developments have no meaning. 

This means that in some way one musti dentify in (4.1) the terms which, 
having singularities in t, would give contribution to many waves in .N-N 
Scattering. For instance, neglect of the high energy contribution in the first 
integral in (4.1) can be misleading if one does not single out the strong # de- 
pendence that such contributions have; i.e. the high energy contribution would 
generate a function of t (not developable) that would prevent a development 
in partial waves. This is also true for the low energy contribution, we have 
seen in Section 3 that the fourth order perturbative contribution, of dia- 
gram I of Fig. 2 is just that term which is strongly dependent on both t and 
w (low lying singularities in both variables). There is a term analogous to 
this one that has instead a strong dependence in w and ¢; this is given by the 
perturbative graph III of Fig. 2. This contribution is considered by GoLp- 
BERGER and OEHME and comes from the two pion contribution to the last 
integral in (4.1). The reason why in such an approach the graph ITI is taken 
into account while I is not included, comes from the asymmetry of eq. (4.1) 
in the treatment of t and f¢. 


We wish to thank warmly Professor S. FUBINI for his continuous interest 
in this investigation and many enlightening discussions. We thank also Pro- 
fessor M. Frerz for useful discussions; one of us (E.L.) wishes to thank him 
for his hospitality at the Theoretical Study Division of CERN. 


APPENDIX I 


The invariant operators in the nucleon-nucleon and nucleon-antinucleon elastic 
scattering. 


In order to build up all the possible independent invariant operators for 
NN and NN elastic scattering, we have at our disposal two y matrices 
[y" and y?, to be saturated between spinors as in (2.16) ] and the three inde- 
pendent vectors (2.1): P, N and A. It is clearly possible to construct from 
these five 4-vectors quite a few scalar quantities (matrices in the spin space 
of the two fermions) but many of them can be further reduced, by using the 
Dirac equation (2.17), to simpler ones. At the end of this process of reduction, 
we are left with a set of 16 possible invariants [this can be most easily seen 
by using the separation (as given by MICHEL ("‘)), of all the forms that can 


x 


(14) L. MicHEL: Proc. Phys. Soc., A 63, 514 (1950); Compt. Rend., 232, 391 (1951); 
Journ. Phys. et Rad., 12, 793 (1951). 
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be obtained by using two y matrices, into five « classes » characterised by 
«class» number Q (Q = 0, 1, 2, 3, 4 indicating that the two y matrices combine 
to give, respectively, a scalar, vector, tensor, axial, pseudoscalar matrix) |. 
These 16 invariants belong to three different classes [as we have no inde- 
pendent pseudovector at our disposal, no combination of y matrices present 
in the classes with Q —3 and Q — 4 can be saturated in order to give even- 
tually a scalar]; five of them are the Fermi invariants (Q — 0) and the others are: 


| MQ) = CE iy? È N) : (iy 5 P 1?) 3 (y" aye dy? 5 N) : (y” vy” È P y?) ; 


(AI.1) | (ytyt prypiy? N); (yrygiy™ Py?y2); (vers? W); (yey Py). 


Q=2: (yt Piy? N); (tip Pyriy?: N); yin Py? N). 
The requirement of time reversal invariance (as expressed, for instance, in 
G.N.O.) allows us to rule out (yt y2iy?-N) and (ytiy"- Py?); time reversal inva- 
riance plus use of the Dirac equation helps us to reduce (y" y?iy?-N), (y"iy" Py”), 
(yeyeyrysiy? N) and (y"yziy"-Py?ys). 

The invariant (y"iy"-P yviy?- N) can be expressed in terms of the remaining 
ones by deyeloping y" and y? in terms of the four orthogonal 4-vectors defined 
by (2.4): M, K, 4 and L. We are therefore left with the five Fermi invariants 
and four other invariants: 


() = ils (ae wy? N) : ay" -P 1”) È 


DI 
er?) Q=2: (yr Pip": N); (ygiy™ Pygiy?: N). 

If we assume now charge symmetry, only the symmetric combination 
(1° iy?-N + iy"-P1?) can be present in the development of M into invariant 
operators. Thus we now have the five Fermi invariants and three other inva- 
riants, which latter coincide with three of the invariants introduced and used 
by G.N.O. (the other two being S and P). These invariants cannot be reduced 
further by the use of the Dirac equation; they are however not linearly inde- 
pendent and can be expressed always in terms of only five linearly inde- 
pendent invariant operators. This can be done in two ways: either expressing 
(1° ty? N + iy"-P1?), (iy"-Piy?-N) and (ytiy"-P y?iy?-N) in terms of Fermi 
invariants, or expressing V, 7 and A in terms of the G.N.O. invariants. The 
former way is, at least in principle, trivial, but gives rise to a rather cumber- 
some calculation; as an even number of Fierz transformations can be applied 
to an invariant without altering it, a repeated application of the Michel tables 
for Fierz transformations of invariants belonging to different Q can be made, 
followed by a suitable application of the Dirac equation. At the end, our 
three G.N.O. invariants will be reduced to linear combinations of Fermi inva- 
riants. To go the other way round, we can express each y matrix in V, T 
and À in terms of the four orthogonal 4-vectors (2.4); a repeated application 
of the Dirac equation will now allow us to express the three Fermi invariants 
we started from in terms of the G.N.O. invariants. 


It is therefore clear that the NN and NN elastic scattering can be always 
described in terms of only five independent invariant funetions, miele 
shall take, as explained in the text, as linear combinations of Fermi invariants. 


co 
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The connection between the Fermi and the G.N.0. sets is given by the fol- 
lowing table: 


St 


V =1/(wt) [m2(t — w) G, — mt + w) 4, dm) PI, 


(t — w) 4m'\ _ m (m?(w + t) o 
Tf =—4)|—_|1——__|4,+ — = — ty — 
| At | TAN OT | wi I) di 


da mei — w) a Va (w — ft) 
(AT.3) wi + wi * A nd 


ae m>(t ra: w)? 
wtt 


NA UP SNCS ay PA?) Co (ty PEN) and, = (yty® Py Ni 

From Table (AI.3) and the definition (2.19) of the perturbative invariants 
P; it is easy to find out the transformation table (2.21) between the pertur- 
bative and the Fermi set of invariant operators. 


APPEND Ix It 


Separation of the S and P-wave two pion contribution to the unitarity condition. 


The absorption matrix (2.13) can be expressed in terms of two independent 
scalar functions A and B: 


(AIT.1) t=—Atwy OB. 


We want now to separate, out of A and B, those parts that contribute 
to the nucleon-antinucleon absorption into S and P two pion states, from 
the rest, which gives absorption into two pion states with angular momentum 
higher than 1. In order to do so, let us write A and B in the form 


(AIT.2) di (LZ — x) + AY, 5 le (4 — B) SE Db 


where ./ and 4 are defined as in (3.25) and A,, and B,, contain all of the 
S and P-wave two pion contribution. The definition of x and f has therefore 
to be such that the S and P-waves coming from them have to cancel exactly 
those coming from / and 4. 
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The § and P-wave contributions of (— + iy:QZ) to 7 will be propor- 
tional to: 


(AIT.3) Je A+ iy Q Bao, je A+ iy'Q BZ) 0,49. 


Expanding ./ and Z into Legendre polynomials, the coefficients in the 
development being given by: 


+1 +1 
(AIT.4) %(t) = |v P,(cos g) d cos gp, B,(t) SL P, (cos g)d cos , 
N =i 


a 


it is easy to see from (ATT.3) that only that part of ./ which does not depend 
on the nucleon-pion angle [proportional to (s— s)] will contribute to the S-wave, 
and only that part of ./ linearly dependent on — and that of # indepen: 
dent of — (s — s) will contribute to the P-wave. We shall therefore define x(t), 
x,(t) and B,(i) by means of the relations: 


[a dQ = [auto dD = ma,(t), 
(AIT.5) "È ; % 
|(— 1 + iy-0 20,40 = | i a +0 à Qu49, 


C4 


[we have, of course, % = a(t) + a, (t)((s —8)/4), B = B(t)]. 
From this definition it is easy to get a, «, and 6, in terms of oY and Z: 


(AIT.6) 


This result could also have been derived starting from the expressions for A 
and B in terms of partial waves that can be obtained from FRAZER and FULCO (5). 
À me shall have to deal, in using the unitarity condition, with the integral 
it over two pion intermediate states; if now we write down explicitly 
the whole of the interference terms present in this integral (by which we 
mean those terms containing products of ~~, 4, « and B with AG, OG Bap). 
it will be quite straightforward to see that it is zero, as it is indeed required 


on physical grounds. It is enough for this to note that relation (AIT.5b) can 
also be written as: 


a ; SES fh 
(AIT.7) I A ) Qua =|@-Miy:00,02. 


(°°) W. R. Frazer and J. R. Futco: Phys. Rev., 117, 1603 (1960). 
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We have in so far disregarded isotopic spin. It is easy to see that the Pauli 
principle for the two pion states implies: 


(AIT.8) AD BY = 0, Ls] BO}. 


RIASSUNTO 


Si studia l’interazione nucleone-nucleone ad energie inferiori alla soglia per process- 
anelastici, usando un metodo basato sulla rappresentazione bidimensionale dell’ami 
piezza di diffusione, nella forma recentemente proposta da Cini e Fubini. L'idea essen- 
ziale di questo metodo è quella di tener conto interamente delle singolarità presenti 
nell’ampiezza di diffusione, per valori delle variabili vicini a quelli corrispondenti alla 
regione fisica (regione la cui ampiezza dipende naturalmente dall’energia dei nucleoni 
incidenti); per ottenere questo, si sfrutta completamente la simmetria della rappre- 
sentazione di Mandelstam. Le funzioni spettrali vengono calcolate facendo uso della 
condizione di unitarietà sia nel canale nucleone-nucleone che in quello nucleone-anti- 
nucleone. In quest’ultimo canale si calcola il contributo dovuto a due pioni intermedi, 
in termini della diffusione pione-nucleone e di alcune ampiezze corrispondenti a bassi 
momenti orbitali nel processo di annichilamento in due pioni. Si discute poi sulla 
costruzione di equazioni integrali per le ampiezze di diffusione corrispondenti alle varie 
onde parziali. L'ampiezza totale di diffusione nucleone-nucleone è separata in operatori, 
ciascuno dei quali è rappresentato da una combinazione lineare di operatori di Fermi, 
combinazione scelta in modo tale che i corrispondenti coefficienti hanno semplici pro- 
prietà di scambio e soddisfano una rappresentazione di Mandelstam al quarto ordine 


perturbativo. 


d 
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The Angular Distribution of Secondary 7-Mesons 
from 4.5GeV =z and 6GeV Proton Interactions in Emulsion. 


H. H. Ary and C. M. FISHER 
H. H. Wills Physical Laboratory, University of Bristol - Bristol 


(ricevuto il 30 Marzo 1960) 


Summary. — The angular distributions of x-mesons produced in two 
hundred 4.5 GeV 7- and two hundred 6 GeV proton interactions with 
emulsion nuclei have been investigated. It is found that the pions from 
the 6 GeV proton induced interactions are consistent with an isotropic 
distribution in the O-system whilst the pions from the 4.5 GeV 77 induced 
interactions have an anisotropic distribution peaked in the forward 
direction. The degree of anisotropy is larger with smaller N, values. 


In a previous paper, ALY et al. (1) hereafter referred to as paper I, an 
analysis of the interactions of 4.5 GeV 7-mesons with emulsion nuclei was 
described. A similar experiment has been undertaken to investigate the nature 
of 6 GeV proton interactions. As a result of these experiments an interesting 
property of the angular distribution of the secondary x-mesons in the centre- 
of-momentum systems, hereafter referred to as the C-systems, of the inter- 
actions has been revealed. It is the purpose of the present paper to describe 
these results. 

The observations were made in two stacks of stripped Ilford G-5 emulsions, 
one exposed to the external 4.5 GeV z- beam and the other to the internal 
6 GeV proton beam of the Berkeley Bevatron. In each experiment 200 inter- 
actions were found by «along the track » scanning. The interactions, recog- 
nized by the production of heavily ionizing tracks or at least two shower tracks, 
were analysed in two groups having N Ae 7 where N, is the number of 


(*) H. H. Any, J. G. DurHIE and C. M. FisHeR: Phil. Mag., 45, 993 (1959). 
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heavily ionizing prongs, i.e. tracks having b* > 1.4 where b* is the ratio of the 
blob density of the track to that of a beam track in the vicinity. 

The procedure adopted in both experiments was the same as was described 
in paper I. The angles of emission of all, shower particles, having b* < 1.4, 
were determined and a sample of shower particles from each experiment iden- 
tified using the method described by DANIEL et al. (2) depending on the va- 
riation of the ionization with velocity in nuclear emulsion. The results for 
the 4.5 GeV x- interactions are given in «I». In the 6 GeV proton induced 
interactions 48 tracks were investigated of which 14 were identified as protons, 
31 as pions and 3 could not be identified. Of the tracks identified as protons 
all but one had angles of emission of less than 10° with the primary direction. 
The laboratory distributions were thus readily corrected for the « proton con- 
tamination» in the manner de- 


scribed in « I ». 

Two methods were adopted to 
investigate the C-system angular 
distributions which do not require a 
knowledge of the Lorentz factor y, of 
the C-system. However it is assumed «| 
that the particles are created in a 
single collision system. 


De) der 


i) The Duller-Walker (3) dis- 
tributions in terms of the varia- 
bles log F(0)/(1—F(0)) and log tg 6, 
where (0) is the fraction of the 
total number of particles having 
angles of emission less than 0, are 
shown in Fig. 1. For an isotropic 
distribution in the (-system the 


Fig. 1. — Log(Fj(1—F)) plots for: & 

a) 6GeV proton interactions having 

N,<7; bh) 6 GeV proton interactions 

having N, > 7; c) 4.5 GeV x interac- 

tions having N, < 7 and a, SEEN Ta 01 sO a 
interactions having N, > 7. o 


(2) R. R. DANIEL, J. H. DAVIES, J. H. MuLvery and D. H. PERKINS: Phil. Mag., 


43, 753 (1952). | 
(3) N. M. DuLLeR and W. D. WALKER: Phys. Rev., 93, 215 (1954). 
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relation between these variables will be 


F(0) x lo | 
1— FO) 2 to (1 2(1 — F(6))) ’ 


(1) 2 log y, tg 0 = log 


where 8,/B,=1-—-ò and B, is the velocity of the C-system in the laboratory 
whilst B, is the velocity of the shower particles in the O-system. Values of B., 
were obtained from the scattering measurements as described in paper I. No 
dependence on the angle of emission 0 was observed and the mean value 
obtained for shower particles identified as m-mesons were f* (4.5 GeV 7-) = 
= 0.865 + 0.024 and B* (6 GeV protons) = 0.88 + .03. Using the value of £. 
for a single nucleon collision partner we obtain 6 ~ 0.05 in both experiments. 
If the relation (1) is plotted then a curve is obtained similar to curve A2 of 
Duller and Walker (5). This is well represented by a straight line of slope 
~1.9. The results of the experimental distributions are given in the Table I. 
The effect of missing due to pions emitted backwards in the C-system and 
appearing as «grey » tracks in the laboratory is taken into account in these 
distributions by increasing the total number of tracks by 5% of the number 
of grey prongs. This is a small correction based on the results of KAMAL (1). 


TABLE I: 
| o S 12 
4.5 GeV r- No AT 4.02 1.48 +.2 = 
N, <7 .52 +.03 1.37 2.2 LOL 
6 GeV protons N,>7 Chil tee 1.87 +.2 — 
NT .44+.02 1:85+.2 Da 
o = dispersion of differential distribution; 
S= slope of Duller-Walker distribution; 
P= Pearson probability of agreement with isotropy. 


The distribution of pions from the 6 GeV proton induced interactions is 
consistent with isotropic emission in the C-system. The distribution of pions 
from the 4.5 GeV x- induced interactions shows a deviation from isotropy. 
It is interesting to note that the deviation trom isotropy is larger for smaller 
N, values. This may indicate that in high N, interactions secondary processes 


are involved which tend to make the distribution more isotropic. In the 6 GeV 
proton case this effect is absent as might be expected. 


(4) A. Kamat: Thesis (University of Bristol, 1957). 
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ii) A further method of investigating the C-system distribution is pro- 
vided by the differential angular distribution in the variable x = log tg 0 
(LINDERN (5)). If the distribution in the C-system is isotropic this plot will 
have dispersion o = 0.394 about its mean value provided 6 = (8:/B.) —1 —0. 
If 6 #0 the dispersion will be changed by a factor (1— 6) Thus for 


6 = 0.05, o = .384 for an isotropic distribution. The experimental distribu- 
tions are shown in Figs. 2 and 3. 


| ol À 
A ov A 
SU 
| 
> 
= HD 
‘810 È 
=| è 
rs 
ers d= 
-1.0 0 +1.0 
X= Log Tg 6, 


e 
20} 
É 2 
5 È 
2 S 
510 x 
© (S) 
= © 
= 
-1.0 0 +1.0 
x= Log Tg 6, X=Log 194 
Fig. 2. — Differential distribution in Fig. 3. — Differential distribution in 
log tg 0, for 6 GeV proton interactions log tg 0, for 4.5 GeV pion interactions 
having: a) N<"7 and 6) N,>7. The having: c) N,<7 and d) N, <7. The 
mooth curves are Gaussians having the smooth curves are Gaussians having the 
same dispersions as the histograms. same dispersions as the histograms. 


It is apparent that the distribution from the proton interactions is con- 
sistent with isotropic emission in the C-system for events having N,> 7 
and deviates only slightly from isotropy for events having N, <7, whilst 


(5) L. von LINDERN: Nuovo Cimento, 5, 491 (1957). 
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the x—-meson distributions are anisotropic the deviations from isotropy being 
largest for small N, events. The corrections for «grey » pions would tend to 
increase the dispersions because they are emitted backwards in the O-system, 
however this correction is small and is not included in these results. 

It was shown in paper «I» that for low N, events at least, the median 
angle of the laboratory distribution is consistent with the assumptions of 
a) single nucleon collision partner and 6) that equal numbers of particles are 
emitted, on average, backwards and forwards in the C-system of a 4.5 GeV 
2- nucleon interaction. A similar result has been obtained for 6 GeV proton 
interactions. If we assume that, for events having N,<7 the interaction 
is with a single nucleon then the angles may be transformed to the C-system 
using the relation 
(2) vate 0}= Di Ü a 
B./Bs + cos 6 
where 9 is the laboratory angle of emission and 0 is the corresponding angle 
in the C-system. Using the values of BY deduced from the scattering meas- 

urements as described in «I», the 
ee pi laboratory distributions, corrected 


ss le for the proton component were trans- 
formed to the O-system. The distri- 
E butions obtained are shown in 

ie 20f 7, Fig. 4. 
= | Clearly there will be a loss of 
| particles emitted backwards in the 
10} O-system and appearing as grey 
tracks. From considerations of the 
0 ee es ee nets o velocity transformation it is evident 
HO 40/69 0,2 0! 20 2060 that all particles emitted with cos 6 
oe greater than — 0.6 appear as shower 
n _ particles in the laboratory. Thus any 


correction for «grey » m-mesons is 
30 30 restricted to the region — .6 > 
> cos? >—- 1.0. The correction for 


A 20 20 
= 
10 | | 10 
Fig. 4. — Angular distribution of secon- 
dary pions in the C-system of an inter- 
0 MINT lp action with a single nucleon: a) 4.5 GeV 
AO SUO O2 002 446. =O m interactions having N, <7; b) 6.0 GeV 
ees proton interactions having N, <7. 
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grey tracks estimated from the results of Kamal is shown in Fig. 4 as a 
shaded region. We may perform a y? test on the part of these distributions 
having — .6< + cos ÿ < +1.0 to examine their departure from isotropy. The 
result is that the proton distribution gives agreement with isotropy on the 
30% level whilst the pion distributions yields a probability of less than 1%. 
This confirms the previous indications that the pion distribution is anisotropic 
whilst the proton distribution is consistent with isotropy. 

In Section 10 of paper «I» the angular distribution of shower particles 
in the effective C-system of a composite 4.5 GeV pion nucleus interaction was 
given which includes events of all N,. The transformation was chosen such 
that the distribution has equal numbers of tracks in the backwards and forward 
directions. This distribution therefore does not represent the distribution in’ 
the O-system of a pion-nucleon interaction. 

It is clear from Fig. 4 that the anisotopy of the distribution of secondary 
pions from the pion induced interactions is in the form of a peak in the forward 
direction, a result which may explain why the collision partner for this group 
of events (7< N,) was found to be less than one proton mass unit (see pa- 
per «I »). This distribution is also in agreement with that obtained by the 
Berkeley group working with a hydrogen filled diffusion chamber (MAENCHEN 
et al. (°), WALKER (")). 

To summarize the results of this investigation we therefore have that the 
C-system angular distribution of secondary pions produced in 6 GeV proton 
interactions with emulsion nuclei deviates only slightly from isotropy whilst 
the distribution of pions from 4.5 GeV #--interactions is anisotropic with a 
peak in the forward direction. 


Discussion. 


KALBACH et al. (8) have investigated the angular distribution of secondary 
pions from proton-proton interactions in emulsion at 6 GeV. These authors 
obtain an anisotropic distribution peaking in the backwards and forwards 
directions of the C-system, however the anisotropy was found to be smallest 
where the statistics are best. It is evident that if a small anisotropy exists 
in the elementary nucleon-nucleon collision then it may be reduced by se- 
condary scattering interactions within the nucleus to give the result observed 


in the present experiment. 


(6) G. MAENCHEN, W. B. Fow Ler, W. M. PoweLL and R. W. WRIGHT: Phys. 
Rev., 108, 850 (1957). 

(7) W. D. WALKER: Phys. Rev., 108, 872 (19517): | 

(8) R. M. KarBac®, J. J. Lorp and C. H. Tsao: Phys. Rev., 113, 330 (1959). 
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Results obtained by FRIEDLANDER (°) in an investigation of 9 GeV proton 
interactions in emulsion are consistent with isotropic emission in the O-system. 

The results may easily be understood in terms of the isobar model of Stern- 
heimer and Lindenbaum (1). If the result of a proton nucleon interaction is to 
produce, on average, two isobars which decay isotropically in their rest systems 
then the resulting C-system distribution will be approximately isotropic if 
the velocities of the isobars in the 0-system are small compared with the velo- 
cities of the pions into which they decay. If the velocities of the isobars in- 
crease with the available energy quicker than the velocities of their decay 
products, then at higher energies the distribution will become anisotropic. 
Observations by the Polish Group (C10K et al. (1)) have provided strong evi- 
dence in favour of a two centres model for some high energy interactions. The 
application of such a model has been discussed in some detail by Coccont (©), 
KRAUSHAAR and MARKS (1%) and others. 

The observed anisotropy in the 4.5 GeV x- interactions can be explained 
if the process is one of inelastic excitation of the target nucleon with a recoil 
pion. The recoil pion will give rise to a peak in the forward direction whilst 
the isobar decays isotropically in the C-system. This model was discussed in 
detail by STERNHEIMER and LINDENBAUM (1°). 
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RIASSUNTO (*) 


Si sono studiate le distribuzioni angolari dei mesoni x prodotti in duecento inte- 
razioni di x di 4.5 GeV e in duecento interazioni di protoni di 6 GeV con nuclei della 
emulsione. Si trova che i pioni da interazioni indotte da protoni di 6 GeV sono coe- 
renti con una distribuzione isotropica nel sistema ©, mentre i pioni da interazioni 
indotte da =~ di 4.5 GeV hanno una distribuzione anisotropica con picco nella dire - 
zione anteriore. Il grado di anisotropia aumenta al diminuire dei valori N,;. 


(*) Traduzione a cura della Redazione. 
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(La responsabilità scientifica degli seritti inseriti in questa rubrica è completamente lasciata 
dalla Direzione del periodico ai singoli autori) 


Proton Magnetic Resonance Study of Hydrates by Spinning Crystals. 


R. CHIDAMBARAM 


Department of Physics, Indian Institute of Science - Bangalore 


(ricevuto l11 Maggio 1960: 


Sample spinning has been employed in nuclear magnetic resonance by earlier 
investigators with different objectives in view as, for example, to overcome par- 
tially field inhomogeneity in high-resolution work (1:2), to reduce the proton spin- 
lattice relaxation time in solid chlorobenzenes for particular proton magnetic reso- 
nance frequencies (3), and to prove the invariance of the second moment of a reso- 
nance line to motional narrowing (45). ANDREW has also indicated that the obser- 
vation of the spectrum on a rapidly rotated specimen should assist in determining 
the sources of spectral broadening (°). 

It is the purpose of this note to suggest that sample spinning can also be of 
assistance in structural investigations of crystals. An important application is to 
hydrates in which, by this technique, the resonance of the water molecule protons 
may often to a large extent be isolated from the resonances due to single protons 
or reorienting ammonium ions, if present. The basis of the method is that if a 
single crystal is spun about an axis, about which the nuclear magnetic resonance 
is not symmetric, then, in the vicinity of resonance, the r.f. voltage across the 
sample coil would be modulated at twice the frequency of the spin even in the 
absence of an external field modulation. The frequency of the spin should be much 
lower that the width of the resonance so that there may be no effects due to mo- 
tional narrowing. 

Let us consider the case of an identical interacting spin- pair (for example, 
the protons in a bound water molecule), all the pairs in the crystal being similarly 


oriented in space. The line shape function can be written (7), assuming Gaussian 
broadening, as 


(1) g(h) = VE >, exp [— {h + (—1) «(3 cos? cos? 6 — 1) 22627 
(*) F. BLocw: Phys. Rev., 94, 496 (1954). 
(@) W. A. ANDERSON and J. T. ARNOLD: Phys. Rev., 94, 497 (1954). 
@) D. E. WoESsNER and H. S. GurowsKy: Journ. Chem. Phys., 29, 804 (1958). 
(a) E. R. ANDREW, A. BRADBURY and R. G. EADps: Nature, 182, 1659 (1958). 
() I. J. Lowe: Phys. Rev. Lett., 2, 285 (1959). 
(*) E. R. ANDREW: Arch. Sci. (Geneva), 12, 103 (1959). 
(?) G. E. PAKE: Journ. Chem. Phys., 16, 227 (1948). 
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where h is the difference between the applied field H and the resonance field H* 
for an isolated nucleus: (H — H*); x=$yr-*; y is the magnetic moment of the 
nucleus and r is the pair separation; 5 is the angle between r and the plane con- 
taining the static field H and perpendicular to the axis of the r.f. coil: œ is the 
angle between the projection of r on the above plane and H; Bis the broadening 
factor which is assumed to be independent of g in the following discussion. 

If now the erystal is rotated about an axis coinciding with the axis of the r.f. 
coil at an angular frequency @, (<1/T,), g in eq. (1) may be replaced by o,f. 
If we write 


dx 
h, = (1 + a — 7 008" over. 


3% E SA 
he =|— h nor [A 2p , 


e = 3a cos? 0/24/28 , 
eq. (1) reduces to 


1 | 
(3) g(h, t) = TE D exp[—{h;— e cos 20,1}? ]. 
= iîi=1,2 


Using a phase-sensitive detector with a reference synchronised to the fourth har- 
monic of the frequency of spin, we can obtain an output signal at this frequency. 
The amplitude of the fourth harmonic can be shown from eq. (3) to be 


(4) SA) = = > exp [— hj] Die sed (— ALS 
28427 i=1,2 rat k=0 

where 

Ò ip (2h —1)(2k + 1) wos (27 -— 1) 

(5) F7 9rk-1(y ate 1)!Æ!(r— K&)! (2k —1) È 


An alternative form for s,(h) is the following: 


1 (h dl” 
Br/ 2x D ll {> È n 2° d(&/2)" (exp [— e?/2] 1612). 
i=1,2 ! /2 


(6) S4(h) = 


where J, is the modified Bessel function of the first order. It can be seen from 
eqs. (4) and (6) that s,(h) has two maxima at 


(7) h = + «(1 —# cos? 0), 


(provided we neglect the very small effect of the presence of one maximum on 
position of the other) and that the height of each peak is given by 


1 
(8) sa exp [— e°/2]1,(e°/2) . 


2315 


108 R. CHIDAMBARAN 

1.00 | The peaks correspond to re- 
| sonances at p= +45, bon 
De at constant 6 in a static 
0.90 crystal, as is to be expected 
| on physical considerations. 
0.25 The assumption of Gaussian 
broadening with a  inde- 
2 O pendent of g is unlikely to 
20.251 cause a significant error In 
x eq. (7). It is preferable to 
-050F 7 use the fourth harmonic 
re rather than the second har- 
EU monic because the latter 
1.00 | L does not appear to have 
-10 .8 _6 .4 2 0 2 4 6 8 10 maximum amplitude at ah 
h in Gauss which is a convenient fune- 

Fig. 1. — The fourth harmonic s,(h) output signal shape for a tion of « and è. 
hydrate proton pair with a = 5.36 G and B—1.5 G. The broken Fig. il gives the s,(h) 


curve is the sum of two Gaussians with B—1.5 G and centred 


Sa output curve calculated from 
at the peak positions of s,(h). 


eq. (4) for x=5.36 G (corre- 

sponding to the protons in 
a water molecule), a broadening factor f=1.5 G and §6=59° 5’ (e=1). The peaks 
are seen to be sharper than those of Gaussians with the same B, indicated in the 
figure by the broken curve. 
Fig. 2 gives the variation of 
peak separation and peak height 
(p) with 6 for the above va- 
lues of « and f. 

The method of investigation 
described above appears well- 
suited for the study of hydrate 
protons in the presence of iso- 
lated protons or reorienting 
ammonium ions since the latter 
have approximately axially 
symmetric resonances in many 
cases and their contributions lo 
to the fourth harmonie output 
signal are likely to be small. 


Peak separation in Gauss 
Relative signal peak strength 


15 30 45 60 75 90 
6 in degrees 


7 Fig. 2. — Variation of the peak separation (——) and the 
It is felt that this method sup- peak height p (-——) with 6 for a hydrate proton pair with 
plementing the normal static #— 5.36 G and 6=1.5 G: 


crystal method may be useful 


in structural investigations and we are shortly undertaking experiments in this 
direction. 


* OK % 
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Several authors have used the clas- 
sical technique of critical absorption 
for establishing bounds for the tran- 
sition energies H,, of mesic X-rays, and 
hence for the meson masses involved ('). 
The K-edge taken as an infinitely sharp 
discontinuity in the absorption coeffi- 
cient o, of the absorber, is used as a 
reference energy #7; one concludes that 
Hiy> By, if the X-ray is «strongly ab- 
sorbed », and #;,< H,, if it is « trans- 
mitted ». 

Actually, the K-edge is not a dis- 
continuous jump of o at the Z,, but 
rather a region where o varies by a 
large factor (~ 5) over a narrow range 
(e.g. ~ H,/10® for Pb). Provided that 
this variation of o is known in detail, 
one can determine an Z;; falling in this 
range very accurately (say, to +H,/10*) 


(*) Research supported by a joint program 
of the Office of Naval Research and the U.S. 
Atomic Energy Commission. 

(**) National Science Foundation predoc- 
toral fellow, 1959-60. 

(*+*) Natioaal Science Foundation 
fellow, 1959-60. Temporary address, 
Geneva, Switzerland. 

(1) For a summary, see D. WEST: Itep. on 
Progr. Phys. London, 21, 271 (1958). 


senior 
CERN, 
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by measuring o (£;;). We have applied 
this principle to the 3D +2P transition 
of mesic phosphorus, which according 
to the early work of KosLov (?) lies 
near E, for Pb. Our work was stimulated 
by the remarks of PETERMANN and 
YAMAGUCHI (3) and made possible by 
the recent work of BEARDEN (4) on the 
Pb K-edge (°). 

Fig. la shows positions of the two 
main fine structure components of the 
3D-2P transition with respect to Bear- 
den’s Pb K-edge for three assumed muon 
masses; the insert in this figure is a fine 
structure level diagram of this 3D+2P 
transition. Fig. 1b shows the absorption 
curves predicted for the three situations 
indicated in Fig. la. These curves illus- 
trate the sensitivity of the method 


(3) S. KosLov, V. FITcH and J. RAIN- 
WATER: Phys. Rev., 95, 291 (1954); S. KosLov: 
Thesis (Columbia, 1954); Nevis Report no. 19 
(unpublished). 

(8) A. PETERMANN and Y. 
Phys. Rev. Lett., 2, 359 (1959). 

(4) A. J. BEARDEN: Phys. Rev. 
240 (1960). 

(5) A preliminary account of the present 
work was given at the 1960 N.Y. Meeting of 
the A.P.S. 


YAMAGUCHI: 


Lett., 4, 
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8 
= 
N 
Lel 5 
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È 
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DI 
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87.8 88.0 88.2 88.4 
E (keV) 


Fig. la. — Positions of the two main components of the 3D — 2P transition of u-mesic phosphorus 
with respect to the Pb K-edge, assuming for My the values: 1) (206.56) mg; 2) (206.74) ma; 
3) (206.92) mg: o are the total absorptionc oefficients as measured by BEARDEN (*), Sapp are 
the «apparent » absorption coefficients defined in the text. The insert is the relevant level 
diagram (not to scale) with the relative line intensities indicated by number in [ ]. 


sa 
ES ld 
3 
S 
L 
= 
re) 
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E 
v 
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© 
a 
= 
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a 
x 
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Fig. 1b. — Curves (1) through (3) are Pb absorption curves predicted from the a 


pepsi i Japp OL Fig. la, 
rt 16 hree muon masses assumed there. The full circles with flags are our experimental points 
with their statistical standard deviations; 


ee : the dashed curve corresponds to assumption (2) but 
includes a correction for 10% Compton degradation in the target. 


N.B. - 1 mil Pb = 2,86-10=2 g/em?. 
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inasmuch as they correspond to varia- 
tions of the muon mass in steps of 0.104. 

u-mesic X-rays in the energy range 
(10--160) keV from a 2 g/em?/P target 
were analysed, with and without Pb ab- 


11I 


peaks of interest amount to about 10% 
and are due to accidentals and the 
Compton pedestals from higher transi- 
tions (mainly 2P +18). They were 
investigated in separate runs with ex- 


sorbers, with a Nal(Tl) spectrometer tended energy range. 
T Hi T T I T x 
3D—+ 2P Pb Absorber 
45 n 
| ies | 
1 2 6 4 °° 
| 6| ‘Nal (TU 
Ò ur '  P.-larget 
SS 
© 2,000 + 0 2 Zones Ov È 
S (123456)=u7 
> (56 Nal(Tl))- 
& 4F 30 Y 
LO 
N | 
fai BAC) 
4 1,000 -— Y _| 
= 0 Absorber 
= a 
Ò 4 mils Pb 
oes mils Pb 
12 mils Pb 
0 | | | all | [ ] ere = 
0 40 80 120 160 
E (KeV)—> 
Fig. 2. — Typical u-mesic X-ray spectra from P observed, without and through Pb absorbers, 


by means of the Nal(TI) spectrometer schemat 
of the small peak at —127 keV to 


(din. diameter, 0.25 in. thick). Fig. 2 
reproduces the spectra recorded with 
0, 4, 8, and 12 mils of Pb; the insert 
shows the experimental setup used. Each 
spectrum in Fig. 2 corresponds to 
1.2-106 x -stops. The energy scale was 
calibrated with 1°Tm (53 and 84 keV 
lines). The zero absorber curve in Fig. 2 
shows an essentially symmetric peak at 
88 keV, corresponding to the 3D —2P 
line sought, while the 12-mils absorber 
curve shows a skew peak near 74 keV 
which we attribute to a combination of 
the Pb-K, fluorescence radiation with 
the transmitted fraction of the 3D—2P 
line. The 4- and 8-mils curves are 
readily understood. as intermediate situ- 
ations. The backgrounds underlying the 
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ically indicated in the insert. The assignment 


a 4-2 transition is tentative. 


The apparent transmission corre- 
sponding to the geometry used was 
extracted from the data in Fig. 2 in 
two independent ways after subtracting 
the background: (A) The spectra in the 
(60—110) keV range were decomposed 
into two gaussians centered about 74 and 
88 keV, of widths corresponding to reso- 
lutions of 24 and 22% respectively, 
experimentally determined for the same 
geometry (°). The area of the 88 keV 
gaussian was taken as the transmitted 
3D—2P intensity. (B) Counts falling 
into the (88100) keV energy band 
were taken as due to transmitted 3D— 2P 


(*) Using «impregnated » targets; see text. 
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photons. A correction of 5% was 
applied in both cases to account for the 
Pb-Kg fluorescence radiation (~ 85 keV) 
which would tend to simulate trans- 
mitted 88 keV X-rays (7). Methods (4) 
and (B) gave within +5% identical 
results and were, hence, averaged to 
following trasmissions: 4 mils 
absorber, (44-+2)%; 8 mils Pb, (21+2)% 
and 12 mils Pb, (12+3)%. The relative 
transmission intensities are plotted as 
full circles in Fig. 1b. The curves in 
that figure were computed using the 
apparent absorption coefficient o,,, (r-h-s. 
scale of Fig. la) for our geometry rather 
than the true absorption coefficient 0,4 
measured by Bearden in good geometry. 
One has 0,4= Capp t+ Oeon Where « repre- 
sents the effect of oblique absorber trav- 
ersals, while o,,, is the attenuation coef- 
ficient for coherent scattering. The rela- 
tion between o;,, and o,,, was established 
empirically by measuring the transmis- 
sion of the radiation from two P targets 
(identical with the mesic target proper), 
uniformly impregnated with the 17°Tm 
and ?7Co (122 keV) activities respectively. 
Only the high energy side of the Nal 
photopeaks was used to exclude the 
effects of the target degradation. Values 
of o,,, at 53, 84, and 122 keV were so 
obtained and extrapolated to the energies 
just above and just below the K-edge, 
yielding 


obtain 


Oapp = (9.8 +0.6) cm?/g 4 


and 


Sapp = (2.60.2) cm?/g, 


respectively. Assuming o,,, to be con- 
Stant across the K-edge, we find 


Coon = (0.4+0.2) em?/g 5 


in reasonable agreement with theory (8). 
The experimental value «x-1=1.35+0.05 


() We are indebted to Prof. G. WENTZEL 
for drawing our attention to this point. For- 
tunately, the Kg fluorescence yield is small. 

(*) K. SIEGBAHN: Beta and Gamma Ray 
Spectroscopy (New York, 1955). : 


J. LATHROP, R. A. LUNDY, V. L. TELEGDI, R. WINSTON and D. D. YOVANOVITCH 


agrees well with the value 1.38 predicted 
numerically for this geometry. 

Fig. 1b shows that the experimental 
points lie closest to curve (2), i.e. that 
the 3Ds >2P; component falls near the 
center of the K-edge. The Compton 
degradation within the target, so far 
neglected, tends to simulate a smaller 
for the situation considered in 
Fig. la. The fraction of photons of 
88 keV degraded by <8 keV in our 
target, i.e. into an energy band that would 
fall under the transmitted «88 keV » 
photopeak, is <20%, while our best 
estimate of this fraction is 10%. Using 
this estimate, one predicts for situa- 
tion (2) in Fig. la the dashed curve 
in’ Fig. 10: 

One obtains the best fit for our 
points, assuming a 10% degradation 
correction, when the 3D;—2P, line is 
assumed to lie at 87992 keV on Bear- 
den’s curve. This corresponds, accept- 
ing the calculation of reference (3), to 
my, = 206.74 m,. 

Uncertainties in this value arise pri- 
marily from our limited knowledge of 
the Compton degradation in the target. 
Assuming either zero or maximum (20%) 
degradation, we obtain the limits 


Capp 


(1) = (206.74 +60) Mo, 


This is to be compared with 


(2) my = (QED) = (206.77+0.013) m, , 
the muon mass derived from Columbia’s 
corrected value (*) for the muon magnetic 
moment, assuming that quantum electro- 
dynamics (QED) predicts the muon 
g-factor correctly. 

Several minor effects have been neg- 
lected in the analysis leading to eq. (1). 
The h.f. splitting (coupling muon-nucleus) 


(*) R. L. GARWIN, G. W. HUTCHINSON, 
S. PENMAN and M. M. SHAPIRO: Phys. Fev. 
Lett., 2, 516 (1959). 
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is significantly larger (1°) (~ 20 eV) than 
the radiation width only for the 3D, + 2P, 
transition. This transition lies, however, 
in a region of the Pb absorption curve 
where the latter is assumed to be con- 
stant for the lack of detailed information. 
Small variations in cx within a few 
hundreds of eV from the K-edge (Kronig 
effect) are indeed known for other ele- 
ments (11) but would not have been 
detected for Pb with Bearden’s resolu- 
tion (2). Whether the 3D;—2P, «line » 


(19) The h.f.s. entry of Table I in ref. (3) 
is distorted by a misprint. In connection with 
the energies of the three lines, the h.f.s. effects 
were not discussed in ref. (3), 

(11) See, for instance, L. G. 
Mod. Phys., 31, 616 (1959). 

(2) A. J. BEARDEN: private communication. 


PARRAT: Fev. 


8 - Il Nuovo Cimento. 
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falls onto such a variation in 0, is an 
open question but, in any case, the h.f.s. 
may decrease the consequences of such 
a coincidence. One may also consider a 
h.f. type coupling between the muon and 
an unpaired electron in the P atom (7). 
This coupling is, however, even for 1s 
electrons about a factor (m,/m.)°(m./ M) 
smaller than the h.f.s. effect just con- 
sidered. 


CODE 


We are grateful to Dr. A. BEARDEN 
for making his results available prior to 
publication. We are grateful to Professor 
G. WEeNTZEL for much constructive crit- 
icism, and to Dr. S. PENMAN for useful 
discussions. 
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The work of West and associates (1) 
has clearly demonstrated the advantages 
of proportional counters (P.C.) over 
Nal(Tl) crystals in the spectrometry of 
low-energy mesic X-rays. For a study 
of the critical absorption of the 88 keV 
3D-2P transition of mesic phosphorus 
— such as has recently been completed 
here (?) by means of a Nal(Tl) spectro- 
meter — the use of a P.C. offers the 
following specific advantages: 


(*) Research supported by a joint program 
of the Office of Naval Research and the U.S. 
Atomic Energy Commission. 

(**) National Science Foundation predoc- 
toral fellow, 1959-60. 

{*-") National Science Foundation senior 
fellow, 1959-60. Temporary address, CERN, 
Geneva, Switzerland. 

(*) D. WEST, R. BATCHELOR and E. F. 
BRADLEY: Proc. Phys. Soc. London, A 68, 801 
(1955); D. WEST and E. F. BRADLEY: Phil. 
Mag., 1, 97 (1956); 2, 957 (1957). 

(?) J. LATHROP, R. A. LUNDY, V. L. TE- 
LEGDI, R. WINSTON and D. D. YOVANOVITCH: 
preceding work. This experiment will be referred 
to here as I. 


(a) Cleaner separation of the photo- 
peak due to the transmitted 88 keV line 
from the photopeak due to the 74 keV 
Pb-K, fluorescence excited in the ab- 
sorbers; 


(b) Smaller contamination of the 
«88 keV » photopeak by 3D—2P radi- 
ation Compton degraded in the phos- 
phorus target. 


Both (a) and (b) are consequences of 
the (3—4) times better energy resolution 
(at 88 keV) given by the P.C. The 
Xe-K, escape peaks which accompany 
the photopeaks of a Xe-filled P.C. share 
their advantages and afford in addition 
the possibility of independent checks. 

P.C. spectrometry should enable one 
to decrease appreciably some of the 
Systematic uncertainties encountered in I 
and thus to obtain a more accurate value 
of the muon mass. We describe here an 
independent measurement of the relevant 


absorption coefficient by this technique. 


This measurement is further based on 
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a more refined determination of the 
effective absorber thicknesses than I 


(see « Calibration », below). 

A P.C. was substituted for the Nal(Tl) 
assembly in a setup very similar to the 
one shown as an insert in Fig. 2 of I. 
It consisted of 2in. diameter, 11 in. 
long Al pipe of 1/32 in. wall thickness, 


provided with a 0.005in. W center 
wire, and filled with 2.5 atm of Xe 
and 2.5 cm Hg of CH,. The two end 


sections of the pipe were wrapped with 
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deteriorated. During the experiment, 
such a purification had to be performed 
every 24 h. 

The P.C. was operated at 4500V, 
and. gave pulses with less than 2 us 
jitter. A prompt coincidence of 
resolution with the signal indicating a 
muon stop in the target was made, and 
used to gate a 100 channel pulse-height 
analyzer (3). The threshold of this gate 
corresponded to an X-ray energy of 
<10 keV, and was kept constant. The 


3 LS 


Ka(Pb) 


805202. 
(74 keV) | 


(88 keV) 


\ “| 


$ ; Background 
War 
vi ttt, 


m 0.004 /n Pb a 
| EE | SE 
: 20 40 60 80 


Channel number 


Fig. 1. Typical u-mesic spectra from phosphorus target observed without and through Pb ab- 


sorber, by means of Xe-filled (2.5 atm) PIO. 


The errors indicated are statistical standard 


deviations. 


Cd foils, leaving an unshielded effective 
volume of only 3 in. length in the middle. 
The P.C. was permanently attached to 
a purifying system which enabled one 
to pass the filling mixture readily over 
hot (350 °C) Ca whenever the resolution 
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integrating time constants of the P.C. 
and of the preamplifier and amplifier 
following it were so chosen that the 


(3) One section of an RIDL 400-channel 
analyzer. 


116 


instantaneous counting rates of the 
experiment (1.5:10% P.C. pulses/s) en- 
tailed no impairment of energy resolution 
nor appreciable dead-time losses. 

Mesic X-ray spectra from a 2 g/em? 
red phosphorus target, in the energy 
range (10114) keV, were recorded with 
Pb absorbers of zero, 4 and 8 mils 
thickness (4). Each run included 3:10 
muon stops in the target, and took 
about 2h. After each run, the P.C. 
spectrometer was calibrated with an 
241Am source (60 keV); and the variation 
in gain was kept to <1% in the course 
of the entire experiment. Typical zero 
and 4 mils absorber spectra, each re- 
sulting from 1.5-10° muon stops, are 
reproduced in Fig. 1. Essentially the 
same peaks as in Fig. 2 of I are observed, 


except that the 3D—2P photopeak 
(88 keV) and the Pb-K, photopeak 
(74 keV) are accompanied by escape 


peaks lying ~ 30 keV (the Xe-K, energy) 
lower. The peak tentatively assigned 
to a 5-3 transition appears in the 
zero absorber spectrum, but overlaps 
with the Pb-K, escape peak in the 
4 mils spectrum. The peaks of interest 
to the transmission measurement are 
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Transmitted intensity —» 


thickness (mils Pb) 


Absorber 


Fig. 2. — Curves (1) though (5) are Pb absorption 
curves calculated from Bearden's values of 64,4 
across the Pb K-edge (corrected for our geom- 
etry), corresponding to the following assumed 
muon masses: 1) 206.62 ma; 2) 206.68 me; 
3) 206.74 meg; 4) 206.80 mg; and 5) 206.86 me. 
The dashed curve corresponds to my = 206.76 me 
and is the best fit to our experimental points 
as quoted in the last column of Table I. 


TABLE I. — List of experimental trasmitted intensities of 3D +2P mesonic transition 
in phosphorus for given Pb absorber thicknesses and of relevant correction factors. 


The last column gives corrected transmitted intensities as plotted on Fig. 2. 
ISIN RI Correction factors: 

Pb nes CASE Corrected 
absorber Li ISS E transmitted 
thickness | (experimental trp I compe (>) ina aio 

values) 
0 mils 1.00+0.04 1.00 1.00 1.00 1.00-+0.04 
4 mils 0.42+0.03 0.95 £0.01 0.94 0.97 +0.01 0.36 +0.03 
8 mils 0.22 10.02 0.90 +0.01 0.89 0.94 10.03 0.17 40.02 
(*) The Compton correction factors correspond to upper limit estimates of the 
a ee final mass evaluation the full magnitude of the effect is allowed as its 
uncertainty. 


(*) Pb sheets identical with those used by 
Bearden in his K-edge measurement [ref, (5)] 
were employed. 


seen to ride on a sloping background, 
the energy dependence of which appears 
to follow the efficiency curve of our P.C. 
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This background is due partly (30%) 
to non beam-coincident events (acciden- 
tals) and partly to prompt photons of 
uncertain origin (e.g. from the de-excit- 
ation of u-capture products, 7r-mesic 
background, e-bremsstrahlung, ete.). 

The data analysis is summarized in 
Table I. After background subtraction, 
the sum of counts in the 3D -2P photo 
and escape peaks was normalized to 
the number obtained with zero absorber. 
To these normalized experimental values 
(column I) we applied corrections for (i) 
the Pb-K; fluorescence radiation (the 
yield of which was determined separately 
with the help of the 120 keV gammas 
from °*Co); (ii) coherent scattering by Pb; 
and (ii) Compton scattering. The rele- 
vant factors used are listed in Table I, 
and yield the corrected intensities listed 
in its last column. These are plotted 
in Fig. 2 as full circles, with errors 
which include statistical standard devi- 
ations as well as systematic uncertainties 
in the listed corrections. 


1. — Calibration. 


To relate the known (5) Pb absorption 
coefficients o,, to our experimental re- 
sults it suffices (neglecting coherent scat- 
tering in the Pb foils and Compton 
degradation in the P target) to deter- 
mine the relative transmitted intensity 7° 
vs. Pb foil thickness x for a y-source 
for which the mass absorption coefficient 
is known, which reproduces the muon 
stop distribution in the P target, and 
for which the detection efficiency in the 
P.C. is identically to that at 88 keV 
photon energy. In fact, knowing T(ox) = 
=T(0) (o dimensionless) one can predict, 
for our experimental set-up, the trans- 
mitted intensity for photons having arbi- 
trary mass absorption coefficients. In 
practice, this calibration was performed 


(5) A. J. BEARDEN: Phys. Rev. Lett., 4, 


240 (1960). 


Gr) 
a 
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with a 60 keV point (Z4LA mM), 
with the Xe pressure reduced to 1.1 atm 
to bring about the same efficieney as 
obtains for 88 keV photons at 2.5 atm. 
The absorption measurements were car- 
ried out for a series of source positions 
throughout the target volume, weighting 
the results with the known muon stop 
distribution. Finally, the absolute ab- 
sorption coefficient o,, for the 60 keV 
radiation was determined in good geom- 
etry. The result of the calibration was 
a plot of T(o) vs. In o, giving an essen- 
tially straight line. The effect of the 
geometry of the setup is summarized 
by a factor f=1.47 (accurate to +4%) 
by which the actual thicknesses of the 
Pb foils have to be multiplied to obtain 
their effective thicknesses (6). A cor- 
rection for Compton degradation in the 
target in the calibration runs affects f 
Diya eal OF" 

Fig. 2 shows transmission curves pre- 
dicted for a set of assumed muon masses, 
using the calculations of Petermann and 
Yamaguchi (7), ref. (5), and this f. We 
emphasize that the plotted experimental 
points are corrected for Pb-K, fluorescent 
radiation, coherent scattering, and for 
Compton degradation in the target. The 
fraction of photons which are Compton 
degraded and which one would tend to 
include in the observed « 88 keV » photo 
and escape peaks is estimated as < 7%. 
The dashed curve in Fig. 2 is a best 
fit to our experimental points, and cor- 
responds to a u- mass. 


source 


(1) m, = (206.76 +453) ID, > 

the asymmetry in the assigned errors 
being due to an allowance for our im- 
perfect knowledge of the actual magni- 
tude of the Compton degradation in the 
P target. The cautionary remarks at 


(°) The work of ref. (*), done in different 
geometry, gave f = 1.35. 
(7) A. PETERMANN and Y. 


Phys. Rev. Lett.. 2, 359 (1959). 


YAMAGUCHI: 


118 J. LATHROP, R. A. LUNDY, S. PENMAN, V. L. TELEGDI, R. WINSTON, ETC. 


the end of ref. (2), concerning the pos- 
sible effects of various minor atomic 
phenomena on this type of experiment, 
fully apply here as well. 

(1) is in excellent agreement with 
the mass value obtained in I, and lies 
even closer to 


(2) m, (QED) = (206.77-+0.013) m, , 


the value derived from Columbia’s 
u+ moment measurement (8), assuming 
that quantum electrodynamies (QED) 
predicts the muon g-factor g,, correctly (°). 
Conversely, accepting (1) and the result 


of ref. (8), one finds 


; CONCA à n e 0.00012 
(3) (9y/2 = IF} = 0.00113 9.90017 ? 


while QED predicts to lowest order 


(4) (gu/2 = Don ca 


a/2% = 0.00116. 


(8) R. L. GARWIN, G. W. HUTCHINSON, 
S. PENMAN and M. M. SHAPIRO: Phys. Rev. 
Lett., 2, 516 (1959). Also Nevis Report no. 79 
(July 1959) (unpublished). 

(°) In comparing (1) and (2), we rely im- 
plicitly on the 7CP theorem which guarantees 
identical magnetic moments for anti-particles. 
Cfr. G. GRAWERT, G. LÜDERS and H. ROLLNIK: 
Fortschr. d. Phys., 7, 318 (1959). 


The best present estimate for the 
r+-ut mass difference is (1°) 


(5) my —m,+ = (66.41+0.10) m, , 
which yields with (1) 


(6) m, = (273.17+0.11) m,, 


in agreement with the accepted value 
(7) Mr = (273.27 +-0-11) me 
and with the limits (1°) 


(273.18+0.11) m.<M7< 


<(273.51+0.04) m,. 
** * 


We are greatly indebted to Dr. R. 
West (Harwell) for much helpful advice 
on proportional counter design, and to 
Prof. C. S. Wu for very useful corre- 
spondence. We wish to thank Prof. G. 
WenrzeL for several illuminating dis- 
cussions, and Dr. K. M. Crowe for the 
loan of a tank of Xe gas. 


(0) Rees 
541 (1957). 


CROWE: Nuovo Cimento, 5, 
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Cosmic Ray Intensity Increase on May 4, 1960. 


O. R. SANTOCHI, J. R. Manzano and J. G. ROEDERER 


Comision Nacional de Energia Atémica - Buenos Aires 


(ricevuto i] 17 Giugno 1960) 


On May 4, 1960, the Ellsworth neu- counting rate, in the interval (1030+ 
tron monitor (*) detected a sudden, bigin- 1045) UT. After this initial peak, 
crease of about 110% in quarter hourly intensity dropped to the previous level 


110 
100% 


790 


Section I }80 Section I 


10%o 
e Li 
s3HRae=3 tei + + + ta 
7 TL dr 122 ae “PA ie He) al 12 È if es 
Fig. 1. — Percentage intensity variation referred to preflare level, for sections I and II of the 


Ellsworth neutron monitor. 15 min interval data. 


——— after about two hours. Ellsworth was 
located in the 09 impact zone for a solar 


(*) Geomagnetic coordinates: Lat. 67° S i Le ‘ 3 
Long. 14.7° E. Operated in collaboration with event occuring at 1030 (1). 
the Instituto Antartico Argentino and the Uni- 


versity of California. 


(@) J. W. Frror: Phys. Rev., 94, 1017 (1954). 
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120 


Fig. 1 gives the 15 min counting rate 
percentage intensity variation of neutron 
monitor sections I and II, referred to the 
preflare level. Unfortunately, only 15 min 
readings are available. Section II shows a 
relative increase systematically higher 
than section I: 


SL 151 "1.100002. 


C(t):32 


500 


200 


100 


+ 6 Qe 
10.45 


10.30 11.00 


Fig. 2. 


ES ONE 15 30 


O. R. SANTOCHI, J. R. MANZANO and J. G. ROEDERER 


additional flux. Indeed, for other time 
variations, like the July 1959 Forbush de- 
creases, which have a flatter spectrum (2): 
percentage variations of sections IT andT 
of the Ellsworth neutron monitor show 
a ratio 1.03+0.02. 

In order to determine the approximate 
time of the beginning of the increase, 
we plotted the integrated additional 


LA 


U.T. 


12.00 13.00 


— Integrated additional counting rate (sum of both sections) for 15 min via The 


extr. £ in th ta ap à j j 1 
apolation of the curve to zero counting rate gives an idea of the time of the beginning of 
the increase. 


According to the geometry of this neutron 
monitor (section II lying above section I, 
Berkeley type), the coupling function 
for section I should be slightly displaced 
towards higher energies than that of 
section IT, due to the increased absorbing 
material above section I with respect to 
section IT. The fact that both sections do 
show an appreciable difference during 
this increase is an indication for a steep 
energy spectrum, or, what is the same 
for a small absorption length of the 


counting rate 
t 


to 


where I is the 15 min counting rate 
and J, the mean counting rate before 


(?) J. G. ROEDERER, O. R. SANTOCHI, 
J. C. ANDERSON, J. M. Carposo and J. R. 
M ANZANO: Preliminary report on cosmic ray 
intensity during magnetic storms in July 1959, 
preprint C.N.E.A., 1959. 
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the increase. Fig. 2 shows the result. 

If we assume tentatively a very fast 
increase, the time of rise lasting not 
more than a couple of minutes, one may 
extrapolate the curve to 0 counting rate 
in order to get some information about 
the probable starting time. We obtain 
the values 1033 and 1039 UT as esti- 
mated lower and upper time limits for 
a sudden intensity rise. Taking this 
into account, the lower and upper limits 
of maximum intensity increase would 
have been of 138% and 277%, respec- 
tively, averaged over the time interval 
between start and 1045 UT. These con- 
siderations, of course, do not exclude the 
possibility of a slower rise starting at any 
time between 1030 and 1033. 
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Finally, neutron monitors of Mina 
Aguilar (geom. Lat. 11.60 S, Long. 3.10 E) 
and Buenos Aires (geom. Lat. 23.308, 
Long. 9.4° E) did not show any abnormal 
behaviour on the 4th of May. As it 
now seems to be usual for all interesting 
cosmie ray events, Ushuaia had a power 
supply failure starting exactly at 1030 UT. 
Ushuaia was located in the 04 impact 
zone. 


* * * 


This report was possible only due to 
the extremely efficient work of our 
technicist, Mr. R. RAsTELLI, who keeps 
operating the Ellsworth neutron monitor 
in spite of innumerable difficulties. 
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Finite Value of the Wave Function Renormalization Constant 
in Quantum Electrodynamics. 


I. BIALYNICKI-BIRULA (*) 


University of Rochester - Rochester, N. Y. 


(ricevuto il 20 Giugno 1960) 


In a recent paper by JOHNSON and 
ZUMINO (1) the question has been raised 
whether there exists a special gauge 
which gives finite value of the electron 
wave function renormalization constant. 
The aim of this note is to prove that 
there exists a class of gauges which give 
a finite (as far as ultraviolet (?) diver- 
gences are concerned) value for this 
constant in every order of perturbation 
theory. It was noticed by Zumino (3) 
that one can get finite value for the 
Z constant (*) in the second order using 
Landau’s gauge defined by the trans- 
verse form of the photon propagator, 
(1) Dik) =(—9,+k,k,k2)k2. 


ev 
The second order result remains un- 
changed if Landau’s gauge is modified 


(*) On leave of absence from Warsaw Uni- 
versity, Warsaw, Poland. 

(1) K. JoHNsoN and B. Zumino: 
Rev. Lett., 3, 351 (1959). 

(*) We will not discuss here the infrared 
divergence, assuming that it is removed by a 
small photon mass À. 

(5) B. ZuMINo: Journ. Math. Phys., 1, 1 
(1960). 

() Z denotes Z, and Zz; which are equal 
due to the Ward identity. 


Phys. 


in the following way, 
(2) D¥(k) = Di,(k) —k,k,k-4*f(e,k), 


where the function f(e, k) is given by a 
power series in renormalized charge, 


(3) f(e, &) = fa(k)e®+ f.(k)et4.... 


Our procedure will be to show that the 
f;(k) coefficients can be chosen in such 
a way that Z is finite in every order of 
perturbation theory. The proof will 
proceed by induction. Let us assume 
that Z is finite up to the 2n-th order 
and let us consider the vertex part 
up to the 2(n +1)-th order. This vertex 
part is a sum of the contributions 
from all irreducible diagrams and from 
skeleton diagrams witb inserted self- 
energies and vertex parts. After mass 
renormalization all inserted electron self- 
energies and vertex parts become finite 
on the basis of our induction assump- 
tion. All inserted photon self-energies 
become finite after charge renormaliza- 
tion because the divergent factor Z, is 
absorbed into the charge at the vertices 
joined by the photon propagator under 
consideration. The arguments show that 
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the divergence is only logarithmic (5) for 
both irreducible and reducible diagrams 
since it appears at the final integration. 
This divergence can be always cancelled 
by an appropriate choice of the function 
În(k) which appears for the first time in 
the calculation of the 2(m-+1)-th order 
correction to the vertex part. The con- 
tribution from f,(k) can be obtained by 
inserting the modified photon propa- 
gator (2) into the lowest order vertex 
part. This term has the form 


(e2)ntt dtk 1 


CT) NES 


(6) p—k=—m 


1 
=e ge km kj,(k) . 


This integral is also logarithmically di- 
vergent so the divergent terms of the 


(®) By logarithmic divergence we mean any 
divergence of the form /d*kk-*(n k°)". 
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vertex part in 2(n+1)-th order can be 
cancelled if we choose the function f,(k) 
in the form of polynomial 


(5) În(k) = Ung + Any In KP + 
+ ...+ ann k?)”, 


with the properly adjusted coefficients 4,;. 
Our argument shows that these coeffi- 
cients are finite. Using this procedure 
we can make Z finite in every order, 
computing step by step the functions f,(k). 
This prescription is not unique because 
we can always change the f,(k) function 
in such a way that the divergent part of 
the integral (3) is not affected. It cor- 
responds to a gauge transformation which 
multiplies Z by a finite factor. 


* OK * 


The author is greatly indebted to 
Dip CO de (origine db (0, Wwe EAN 
E. C. G. SuDARSHAN for many helpful 
discussions and essential suggestions. 


IL NUOVO CIMENTO 


Wit, SOUL, IN Al 


1° Luglio 1960 


Equivalence Principle «Paradox» in the Motion of a Gyroscope. 


L. I. SCHIFF (*) 


Faculté des Sciences - Orsay (**) 


(ricevuto il 4 Luglio 1960) 


In this note we describe an apparent 
violation of the equivalence principle 
that arises in connection with the pre- 
cession of the spin axis of a gyroscope 
moving in a uniform gravitational field. 
The «paradox» is then resolved by a 
realistic specification of what is meant 
by the term uniform. 

The rate of change of the spin an- 
gular momentum S, of a gyroscope, as 
measured by a co-moving observer, has 
been derived in a recent paper (!). If 
the gyroscope moves in the gravitational 
field external to a non-rotating, spher- 
ically symmetric mass m, and is acted 
on by a non-gravitational force F that 
exerts no torque on the gyroscope, the 
spin equation of motion is 


(1) (dS, /dt) =Q~x Sh : 


where 
(2) 92 = (4 /2mye?) (F xv) + 


+ (8Gm/2e?73)(r xv). 


(*) Supported in part by the United States 
Air Force through the Air Force Office of Scien - 
tific Research. 

(**) Permanent address: Institute of Theor- 
etical Physics, Department of Physics, Stanford 
University, Stanford, Cal. 

() L. I. SCHIFF: Proc. of the National Aca- 
demy of Sciences (June 1960). 


Here G is the Newtonian gravitational 
constant, ¢ is the speed of light, my is 
the rest mass of the gyroscope, r is its 
position vector with respect to the center 
of the attracting mass m, and v=(dr/dt) 
The first term on the right side of (2) 
is the special-relativistic Thomas pre- 
cession, and the second term is the gen- 
eral-relativistic de Sitter-Fokker effect. 

It might be supposed that a uniform 
gravitational field can be obtained by 
allowing m and r to become very large 
in such a way that the local acceleration 
of gravity 


(3) g = Gm, 


a definite value. Then if the 


motion is quasi-static, 


retains 


(4) EF = mr, 


Le 
where r=r/r, and 


(5) Q = (2g/c?)(r x v) . 
As a simple example of eq. (5), we 
assume that T is directed vertically 


upward along the positive ¢-axis and 
v is directed horizontally along the pos- 
itive x-axis. Then Q is along the pos- 


itive y-axis, and has the magnitude 
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{2gv/c?). If then S, is vertical and the 
gyroscope moves a distance Ar, the 
upper end of its spin axis tilts forward 
through the angle 


(6) On = AS,| /|S,| = (2g Axe?) J 


According to the equivalence prin- 
ciple, this tilt angle should agree with 
that found when the gyroscope is sub- 
jected to an upward acceleration g in 
the absence of the gravitational field. 
This may be obtained from eqs. (1) 
and (2) by setting m=0, retaining the 
expression (4) for F, and replacing v by 
vigir. It follows that the tilt angle 
measured by a co-moving observer is 
(g Ax/2c?), which is one quarter that 
given by eq. (6). 

This « paradox » may be resolved by 
noting that the gravitational field is not 
strictly uniform so long as m and r are 
finite. The essential parameter in the 
Schwarzschild metric is (@m/c?r), which 
is always less than or of order unity. 
In all practical cases, it is very much 
smaller than unity, and indeed the 
derivation of eqs. (1) and (2) is only 


valid when this is the case. But from 
eq. (3), we see that 
(7) P = (Gm/c?r) = (gr/c?) , 


so that this parameter becomes arbi- 


trarily large if g is fixed and r increases 
without limit. 

A measure of the non-uniformity of 
the field is provided by the angle of 
divergence 0,=(Ax/r) between the field 
lines at the ends of the path of the 
gyroscope. It should be noted that this 
path must cut across the field lines, 
since Q given by (5) is zero if v is 
parallel to Fr. Now 0, may be expressed 
in terms of the tilt angle (6) and the 
Schwarzschild parameter (7), as 0,= 
=(07/2P). Since P<1, we see that 
0 > 07. Thus if the precession of the 
gyroscope axis is to be significant, the 
gravitational field cannot be regarded 
as uniform, and the equivalence prin- 
ciple is not applicable. 

We conclude that the tilt angle is 
actually given by eq. (6). Of this amount, 
one quarter mav be calculated from the 
equivalence principle. The other three 
quarters may be thought of physically 
as the result of an almost completely 
unsuccessful attempt of the spin axis of 
the moving gyroscope to follow the 
changing direction of the gravitational 
field. The coupling between the gravi- 
tational field and the inertial coor- 
dinate system of the gyroscope is very 
weak, of order the Schwarzschild para- 
meter P, so that the tilt angle is 
very much smaller than the divergence 
angle. 
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International Review of Cytology, vo- Gli argomenti trattati nei tredici articoli 
lume VIIT, 1959. Editors: G. H. non sono legati da alcun nesso logico; 
BOURNE and J. F. DANIELLI. vi si possono tuttavia riconoscere alcuni 
Academic Press, New York and degli orientamenti attuali della ricerca 
London. biologica. 

L’estendersi e il perfezionarsi della 
ricerca al microscopio elettronico hanno 
riportato il metodo morfologico a un 
livello di prima importanza nell’interpre- 
tazione dei fenomeni vitali. Le immagini 


L’ottavo volume di questa Rivista 
raccoglie, come i precedenti, una serie 
di sintetici articoli di aggiornamento su 
disparati problemi di biologia cellulare. 
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fornite dal microscopio elettronico, spesso 
soggette fino a poco tempo fa all’impu- 
tazione di artefatti, vengono ora guar- 
date con maggior fiducia, per la gene- 
rale corrispondenza trovata fra queste e 
i risultati di altri mezzi più indiretti di 
indagine. Ne è conseguenza l’interpre- 
tazione funzionale, dinamica, che viene 
oggi facilmente data alle immagini elet- 
troniche. I primi quattro articoli della 
Rivista rispecchiano questa tendenza. 
The structure of cytoplasm di C. OBER- 
LING, è una rassegna delle attuali cono- 
scenze sulla morfologia, particolarmente 
submicroscopica, del citoplasma; in Wall 
organization in plant cells di R. D. PRE- 
STON, viene avanzata una interpretazione 
della struttura e della biosintesi delle 
fibrille che compongono la membrana di 
cellule vegetali, in base a dati ricavati 
dall’analisi chimica, dall’esame al micro- 
scopio elettronico, e da diagrammi di 
diffrazione di raggi X. In The cell surface 
of Paramecium di C. F. EHRET ed 
E. L. Powers, viene discussa l’interpre- 
tazione della complessa struttura rive- 
lata dal microscopio elettronico in corri- 
spondenza alla superficie di un orga- 
nismo unicellulare, mentre in Submicro- 
scopie morphology of the synapse di 
E. De RoBERTIS, è avanzata un’ipotesi 
sulla conduzione nervosa a livello delle 
sinapsi, basata in gran parte sui dati 
morfologici forniti dal microscopio elet- 
tronico. 

Due articoli sul problema della per- 
meabilità cellulare testimoniano un altro 
orientamento della ricerca biologica mo- 
derna: lo sforzo di definire con rigore i 
limiti di applicabilità delle leggi fisiche 
e chimico-fisiche ai sistemi viventi: 
Osmotic properties of living cells di 
D. A. T. Dick, è un tentativo di inter- 
pretare la permeabilità delle cellule ani- 
mali in base alle sole proprietà osmo- 
tiche dei costituenti elementari del pro- 
toplasma; in Sodium and potassium 
movements in nerve, muscle and red cells 
di I. M. GLyNN, viene riconosciuta la 
necessità di ammettere un’attivita bio- 
logica della membrana per giustificare 
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a distribuzione degli ioni di potassio e 
sodio in cellule specializzate; in Trace 
elements in cellular function di B. L. VAL- 
LEE e F. L. HocH, la distribuzione nella 
cellula di metalli, legati in complessi 
metallo-proteici per costituire composti 
dotati di attività biologica specifica, 
viene studiata a mezzo della spettro- 
grafia di emissione, metodo finora poco 
usato nella ricerca biochimica. Anche 
l’estendersi di applicazioni tecniche più 
pertinenti alla Fisica testimonia la sen- 
tita necessità di porre i risultati della 
ricerca biologica in termini esatti. 

I rimanenti articoli trattano 
menti di Citologia speciale (The mam- 
malian reticulocyte di L. M. LOWENSTEIN, 
The physiology of chromatophores di 
M. FINGERMAN, Pinocytosis di H. Hor- 
TER) o di Istologia (The fibrous com- 
ponents of connective tissue with special 
reference to the elastic fiber di D. A. HALL, 
Experimental heterotopic ossification di 
J. B. BRIDGES). A survey of metabolic 
nuclei di 


argo- 


studies on isolated mammalian 
D. B. Roopyn, è una raccolta biblio- 
grafica dei dati sulla composizione e sulle 
proprietà enzimatiche del nucleo, rica- 
vati dall’analisi di nuclei isolati per 
omogeneizzazione e centrifugazione fra- 
zionata. 

G. MARIN 


N. R. HANSON — Patterns of Discov- 
ery. Cambridge University Press, 
1958, pp. 240. 

Lo scopo di questo libro è così 
indicato dall’autore nella introduzione: 
«The approach and method of this 
essay is unusual. I have chosen not 
to isolate general philosophical issues 
- the nature of observation, the status 
of facts, the logic of causality, and the 
character of physical theory — and 
use the conclusions of such inquiries as 
lenses through which to view particle 
theory. Rather the reverse: the inade- 
quacy of philosophical discussions of 


this subjects has inclined me to give a 
different priority. Particle theory will 
be the lens through which this peren- 
nial philosophical problem will be view- 
ed». Dopo una dichiarazione program- 
matica di questo genere si dovrebbe 
richiedere la recensione del volume da 
parte di un filosofo. Tuttavia essendo 
il saggio indirizzato esplicitamente sia 
alle persone con interessi prevalente- 
mente filosofici che ai fisici, un giudizio 
da parte di un fisico appare opportuno. 

L’idea di illustrare e chiarificare con- 
cetti come osservazione, causalità © 
teoria fisica, mediante un processo 
‘costruttivo descrivendo come questi si 
determinano naturalmente nei processi 
mentali del fisico, non è nuova e non 
mancano saggi particolarmente signi- 
ficativi nella letteratura sull’argomento : 
ad esempio alcuni capitoli di Analysis 
of Matter di Bertrand Russel. Ci sembra 
invece che N. R. Hanson non riesca a 
svolgere efficacemente il suo programma 
in questo volume Patterns of Discovery. 
Il difetto principale del libro è rappre- 
sentato da una descrizione eccessiva- 
mente dettagliata di situazioni parti- 
‘colari in cui va perduto il vero scopo 
della discussione. Dato il metodo esem- 
plificativo usato nella trattazione dei 
diversi problemi il lettore sente alla 
fine di ogni capitolo la necessità di una 
sintesi delle osservazioni fatte che tut- 
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tavia manca completamente. Se si con- 
sidera la dichiarazione programmatica 
dell'autore riportata precedentemente si 
è indotti a pensare che un qualsiasi 
tentativo di sintesi in termini logici 
sia stato volutamente evitato dall’autore. 

Al contrario gli esempi usati nel- 
l’illustrare la discussione sono ben scelti 
e appropriati: progredendo nella lettura 
del saggio si nota un po’ con rammarico 
che vi è molto materiale interessante 
male organizzato. Ad esempio nell’ul- 
timo capitolo è posto il problema del 
significato della coincidenza della mec- 
canica classica e della meccanica quan- 
tistica per grandi numeri quantici. Un 
possibile modo di rispondere al problema 
consiste nel riconoscere in questa coin- 
cidenza identità strutturale di «frasi » 
(nel senso della logica simbolica) appar- 
tenenti 3 linguaggi diversi. Questa con- 
clusione è raggiunta in definitiva anche 
dal presente autore: tuttavia il lettore 
vi è portato attraverso discussioni det- 
tagliate dei concetti dinamici della mec- 
canica classica e della meccanica quan- 
tistica il cui effetto è di sviare notevol- 
mente l’attenzione dal problema posto, 

Il Volume è diviso in sei capitoli e 
due Appendici. Gli argomenti trattati 
sono nell’ordine: Observation - Facts 
- Causality - Theories - Classical Particle 
Physics - Elementary Particle Physics. 
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